Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 560,670 C→T 5.2% P325L (CCG→CTG)  sfmD → putative outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913560,6700CT 5.2% 248.4 / 2.6 78P325L (CCG→CTG) sfmDputative outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base C (36/36);  new base T (2/2);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGACGGCGGGCGATTTTCGTAGCGG  >  NC_000913/560524‑560817
                                                                                                                                                  |                                                                                                                                                   
tatGTATCCTTATAGCCTGCAAGGGTTTGCGCCATCTGTTCGTGGGATTGCCAGTACGGCGGCCCAGCTAAGGATTCGGGAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCg                                                                                                                                                    <  6:225995/148‑1 (MQ=255)
   gTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCTGACa                                                                                                                                                 >  5:81974/1‑148 (MQ=255)
       ccTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCtt                                                                                                                                             >  7:214620/1‑148 (MQ=255)
       ccTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCACGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCACTGCGCTTTTTAAATTACAGATTTTCACCAGACATCtt                                                                                                                                             >  4:102306/1‑148 (MQ=255)
        cTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCtt                                                                                                                                             >  7:164847/1‑147 (MQ=255)
          gATAGCCAGAAATTGTTTTCCCCAAAGGTACTTGGGATTGACAGTACGTCGGCCAAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAATCTATGTTTACCAAGGCGCTTTTGAAATTACAGATTTGCCCACGACATCTTCa                                                                                                                                           >  8:324606/1‑147 (MQ=255)
           aTAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCCAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATATGCACCCGACATCTTCaa                                                                                                                                          <  1:234898/147‑1 (MQ=255)
           aTAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCaa                                                                                                                                          <  8:164847/147‑1 (MQ=255)
            tAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCTGACATCTTCAAAt                                                                                                                                        >  3:238434/1‑148 (MQ=255)
              gCCAGCAAGGTTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCTCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATg                                                                                                                                       <  1:116004/147‑1 (MQ=255)
                  gCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAACGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGAt                                                                                                                                  >  5:35477/1‑148 (MQ=255)
                     aGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTg                                                                                                                               <  8:214620/148‑1 (MQ=255)
                      gggTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTgg                                                                                                                              <  5:232126/148‑1 (MQ=255)
                       ggTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGa                                                                                                                             <  1:264389/148‑1 (MQ=255)
                           tGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCACGGCGCTTTTGAAATTACAGATTTGCACCTGACATCTTCAAATGGCGATCTTGACGt                                                                                                                          <  6:81974/147‑1 (MQ=255)
                            gCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTc                                                                                                                         <  6:388732/147‑1 (MQ=255)
                            gCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCa                                                                                                                        >  3:305540/1‑148 (MQ=255)
                             ccccAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCAc                                                                                                                       >  5:89455/1‑148 (MQ=255)
                             ccccAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACAGCTTCAAATGGCGATCTGGACGTCa                                                                                                                        <  4:322145/147‑1 (MQ=255)
                             ccccAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCACCCGGCGCTTTTGAAATTACAGCTTTGCACCCGCCATCTGCAAATGGCGATCTGGAAGTCAc                                                                                                                       >  7:358357/1‑148 (MQ=255)
                             ccccAACGGTACGGGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCAc                                                                                                                       <  8:381590/148‑1 (MQ=255)
                               ccAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCa                                                                                                                     >  5:143422/1‑148 (MQ=255)
                                  aCGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTTACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATcg                                                                                                                  >  8:133406/1‑148 (MQ=255)
                                            ggATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGcgcg                                                                                                         <  2:27947/147‑1 (MQ=255)
                                            ggATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGcgcg                                                                                                         <  2:283335/147‑1 (MQ=255)
                                            ggATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGa                                                                                                        <  7:410424/148‑1 (MQ=255)
                                             gATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGAt                                                                                                       <  4:26262/148‑1 (MQ=255)
                                               ttGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCaa                                                                                                                                          <  5:359730/111‑1 (MQ=255)
                                               ttGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCaa                                                                                                                                          >  6:359730/1‑111 (MQ=255)
                                                 gCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCtt                                                                                                                                             <  2:189623/106‑1 (MQ=255)
                                                 gCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCtt                                                                                                                                             >  1:189623/1‑106 (MQ=255)
                                                  cccGTACGTCTGCCCAGCTAACGATTCGTCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAACTGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCaa                                                                                                  <  3:102306/148‑1 (MQ=255)
                                                  cccGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCaa                                                                                                  <  6:143422/148‑1 (MQ=255)
                                                     gTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCaa                                                                                                                                          >  8:224803/1‑105 (MQ=255)
                                                     gTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCaa                                                                                                                                          <  7:224803/105‑1 (MQ=255)
                                                      tACGGCGGCCCAGCTTACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATcag                                                                                              <  7:133406/148‑1 (MQ=255)
                                                        cggcggCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCTAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACCGCTTTGCACCCGACATCTTCAAATGGCTATCTGGACGTCACCATCGACGAGCGCGATGGCAATcagc                                                                                             <  6:78219/147‑1 (MQ=255)
                                                        cggcggCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATcag                                                                                              >  5:350133/1‑146 (MQ=255)
                                                                 aGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTaca                                                                                    >  2:337989/1‑147 (MQ=255)
                                                                 aGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTaca                                                                                    >  3:198149/1‑147 (MQ=255)
                                                                  gCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATc                                                                                                                                               <  6:365990/87‑1 (MQ=255)
                                                                  gCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATc                                                                                                                                               >  5:365990/1‑87 (MQ=255)
                                                                  gCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTacac                                                                                   >  5:9282/1‑147 (MQ=255)
                                                                   cTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACaa                                                                                 >  2:340670/1‑148 (MQ=255)
                                                                    tAACGATACGGCAAAATCGTATTATTAACTATCAAAGCTATGTTTTCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCCAATGGCGCTCTGGACGTCACCATCGACCAGCGCGATCGCAATCAGCAGAATTAcacact                                                                                >  2:293736/1‑146 (MQ=255)
                                                                       cGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTcc                                                                             >  2:847/1‑148 (MQ=255)
                                                                         aTTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACACTACATCTTCAAATGGCGATATGGACGTAACCATCGACGAGCGCGATGGCAATCATCAGAATTACACAATTCCg                                                                            >  3:71221/1‑147 (MQ=255)
                                                                                 aaaTGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACa                                                                   <  4:305540/148‑1 (MQ=255)
                                                                                 aaaTGGTTTTATTAACTATCAAAGCGATGATTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAATTGGCGAACTGGACGTCACCCTCGACGAGCGCGATGGCAATCCGCAGAATTACACAATTCCGTATTCAACa                                                                   >  2:291209/1‑148 (MQ=255)
                                                                                   aTGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGt                                                                 >  8:215107/1‑148 (MQ=255)
                                                                                           attatCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCATTCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAAttt                                                         >  5:11227/1‑148 (MQ=255)
                                                                                               tctatcAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGCCGTCACCATCGACGAGCCCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTAc                                                      >  6:257209/1‑147 (MQ=255)
                                                                                                   tcAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCTGACATCTTCAAATGGCGATCTGGACGTGACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAAcg                                                  <  4:238434/147‑1 (MQ=255)
                                                                                                   tcAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAAcgc                                                 >  1:177140/1‑148 (MQ=255)
                                                                                                    cAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAAcgc                                                 >  1:318922/1‑147 (MQ=255)
                                                                                                     aaaGCTATGTTTTCCCCGGCGATTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGTCGTGACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGa                                               <  1:293736/148‑1 (MQ=255)
                                                                                                     aaaGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGa                                               <  5:215345/148‑1 (MQ=255)
                                                                                                     aaaGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGa                                               <  6:35477/148‑1 (MQ=255)
                                                                                                     aaaGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATCCAACAGTGCCAATTTTACAAcgcg                                                <  1:337989/147‑1 (MQ=255)
                                                                                                       aGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAg                                             >  5:381013/1‑148 (MQ=255)
                                                                                                       aGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAg                                             >  4:115129/1‑148 (MQ=255)
                                                                                                             gTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGtt                                       >  4:209921/1‑148 (MQ=255)
                                                                                                                tCCCCCTGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTc                                     <  4:71221/147‑1 (MQ=255)
                                                                                                                tCCCCCTGAGCTTTTGAAATTACAGTATTGCACCCGTCATCTTCAAATGGCGATCTTTACGTCACCATCGACGAGCGCGATTGCAATGAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTc                                     <  6:381013/147‑1 (MQ=255)
                                                                                                                  ccccGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCaa                                   <  3:115129/147‑1 (MQ=255)
                                                                                                                  ccccGGCGCTTTTGAAATTACAGATTTGCACAAGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAAATCCGTATTCAACAGTGCAAATTTTACAACGCGAAGTGCGTTTCaaa                                  <  3:67271/148‑1 (MQ=255)
                                                                                                                    ccGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGAGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAAtt                                >  1:292774/1‑148 (MQ=255)
                                                                                                                     cGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTTACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAAc                                                                    <  1:265139/111‑1 (MQ=255)
                                                                                                                     cGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTTACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAAc                                                                    >  2:265139/1‑111 (MQ=255)
                                                                                                                       gcgcTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGa                             >  7:4651/1‑148 (MQ=255)
                                                                                                                             ttGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGAc                       >  4:229329/1‑148 (MQ=255)
                                                                                                                                aaaTTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGAcggc                    <  4:198149/148‑1 (MQ=255)
                                                                                                                                 aaTTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGAcggcg                   <  6:350133/148‑1 (MQ=255)
                                                                                                                                     aCAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAggg                                           <  2:380934/120‑1 (MQ=255)
                                                                                                                                     aCAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAggg                                           >  1:380934/1‑120 (MQ=255)
                                                                                                                                        gATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGACGGCGGGCGAtt            >  7:15030/1‑148 (MQ=255)
                                                                                                                                         aTTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGACGGCGGGCGAttt           <  1:340670/148‑1 (MQ=255)
                                                                                                                                              cACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGACGGCGGGCGATTTTCGTa      <  3:209921/148‑1 (MQ=255)
                                                                                                                                                cccGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGACGGCGGGCGATTTTCGTAGc    <  6:89455/148‑1 (MQ=255)
                                                                                                                                                 ccGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGACTGCGGGCGATTTTCGTAGCg   >  3:200290/1‑148 (MQ=255)
                                                                                                                                                 ccGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGACGGCGGGCGATTTTCGTAGCg   >  7:228328/1‑148 (MQ=255)
                                                                                                                                                  cGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGACGGCGGGCGATTTTCGTAGCgg  >  6:182641/1‑148 (MQ=255)
                                                                                                                                                  cGACATCTTCAAATGGCGATCTGGACGGCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGACGGCGGGCGATTTTCGTAGCgg  <  6:11227/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGACGGCGGGCGATTTTCGTAGCGG  >  NC_000913/560524‑560817

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: