Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,270,143 C→T 9.7% intergenic (+541/+73) yhaC → / ← rnpB pentapetide repeats‑related protein/RNase P, M1 RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,270,1430CT 9.7% 457.7 / 30.3 176intergenic (+541/+73)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base C (90/66);  new base T (4/13);  total (95/81)
Fisher's exact test for biased strand distribution p-value = 9.51e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 184 total aligned reads displayed.

GAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTC  >  NC_000913/3269994‑3270289
                                                                                                                                                     |                                                                                                                                                  
 aaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAATAATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTCTAGcg                                                                                                                                                   >  2:3783892/1‑150 (MQ=255)
 aaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACACTCATTAAAATTTATGCCCCACCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAAGTAAAACAGGTGAAACAGACCGATAAGCCGCGTTCTTTTGGGGATAGcg                                                                                                                                                   >  1:3836/1‑150 (MQ=11)
   aCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtc                                                                                                                                                 >  2:2746140/1‑150 (MQ=37)
     tGACCGATCAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAACAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtcgt                                                                                                                                               <  2:1923277/150‑1 (MQ=255)
        ccGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGa                                                                                                                                            <  2:2815920/150‑1 (MQ=37)
         cgctAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGACCTTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGCAGCGGTTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAc                                                                                                                                           <  2:2376805/147‑1 (MQ=11)
              aGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtca                                                                                                                                      <  1:748641/150‑1 (MQ=37)
                ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACAGAACGATAAGCCGCTTTCTTTTGGGTAAAGCGTAGTGGACACtcaaa                                                                                                                                    >  2:1048061/1‑148 (MQ=11)
                 cGCTTTCTTTTGGGTATAGTGTCGTGGACACTAAATAAAATTTATGCCCCACCAAAATCAAAAAACCCCCCAATCCGGTTTTTAAGTAAAACACGTGAAAAAGACCGATAAAACCCGTTCTTTTTTGTATATCGTCGTGGAACCtcattc                                                                                                                                   >  1:3646420/1‑150 (MQ=2)
                     ttctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCt                                                                                                                               <  1:2803761/150‑1 (MQ=37)
                               tataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAGAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCt                                                                                                                     <  1:876778/150‑1 (MQ=37)
                               tataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCt                                                                                                                     <  1:3537689/150‑1 (MQ=37)
                                tttgtgtCGTGGACAGTCATTCATCTGTCTGGCCCTCCTACAGCAAAATCGCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGCTAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCGCCTc                                                                                                                    <  1:3503570/147‑1 (MQ=11)
                                ataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTc                                                                                                                    <  2:1362640/150‑1 (MQ=37)
                                 taGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGGATAGCGTCGTGGACAGGCCATCAACTTTCAGCCCCGcc                                                                                                                   >  2:3557066/1‑150 (MQ=11)
                                     gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                                                                               <  1:1746948/150‑1 (MQ=37)
                                          ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaaa                                                                                                          <  2:2442099/150‑1 (MQ=37)
                                              aGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCCATCAACATTCTTCCCCTCCACAAGCAAAAAcc                                                                                                      >  2:1168213/1‑150 (MQ=11)
                                              aGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTc                                                                                                                    <  1:3329629/136‑1 (MQ=37)
                                              aGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTc                                                                                                                    >  2:3329629/1‑136 (MQ=37)
                                              aGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCAATCAACTTTATGCCCCTCCAAAAGCaaaaaac                                                                                                      >  2:1387600/1‑148 (MQ=255)
                                              aGTCATTCATCTTTCTCCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTGGGGTATAGCGTCGAGGAAAGTCATTCATAATGCTGCCCCGCCAAAAGCAAAAAcc                                                                                                      >  2:2329240/1‑150 (MQ=11)
                                                 cattcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTTTAGCGTCGTGGACAGTAATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgc                                                                                                   <  1:3945854/150‑1 (MQ=255)
                                                  attcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATGTTTATGCCCCTCCAAAAGAAAAAACccgcc                                                                                                  >  1:1112024/1‑150 (MQ=255)
                                                  attcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgcc                                                                                                  >  1:510481/1‑150 (MQ=37)
                                                  attcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCGTTCTGCCCCTCCAAAAGCAAAAACccgcc                                                                                                  >  2:1230398/1‑150 (MQ=37)
                                                  attcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGAGGACAGTCATTCATCGTTCTGCCCCGCCAAAAGCAAAAACccgcc                                                                                                  >  2:1968783/1‑150 (MQ=255)
                                                  attcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCGTTCTTTTGGGTATAGCGTCGGGGACAGTCATTCATAATTCTGCCCCAACAAAAGCACAAACccgcc                                                                                                  >  2:3042834/1‑150 (MQ=11)
                                                   ttcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgccg                                                                                                 >  1:2597377/1‑150 (MQ=37)
                                                   ttcatCTTTCTGCCCCTCCAAAAGCAAAAACACGCCGAAGCGGGTTTATACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTGTGGGTATAGCGTCGTGGACAGTCATTCATATTTCAGCCCCTCCAAAAGCAAAAACacgccg                                                                                                 >  1:2148799/1‑150 (MQ=11)
                                                      atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAg                                                                                              <  1:2457432/150‑1 (MQ=17)
                                                         tttCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGCCAGTCATTCATCGTTCTGCCCCTCCCAAAGCAAAAACCCGCCGAAGCgg                                                                                           >  1:1132445/1‑150 (MQ=255)
                                                           tCTGCCCCTCAAAAAGCAAAAACCCGCCGAAGCGGTTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGt                                                                                         <  1:1578190/150‑1 (MQ=255)
                                                             tGCCCCTCCAACAGCACCAACCCGCCGCAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                                                       <  2:3721856/150‑1 (MQ=255)
                                                             tGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCGTTCTGCCCCTCCAAAAGCAATAACCCGCCGAAGCGGGttt                                                                                       >  1:633212/1‑150 (MQ=255)
                                                              gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGGCAGTCATTCATCTTTCTGCCCCGCCAAAAGCAAAAACCCGCCGCAGCGGGtttt                                                                                      >  1:2987209/1‑150 (MQ=255)
                                                              gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGtttt                                                                                      <  1:2746140/150‑1 (MQ=255)
                                                              gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtttt                                                                                      >  1:738844/1‑150 (MQ=37)
                                                                cccTCCAACAGCAAAAACCCGCCGAAGCGGGTTTTTACGTACAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTa                                                                                    <  2:3632784/150‑1 (MQ=12)
                                                                cccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTa                                                                                    <  2:2466041/150‑1 (MQ=37)
                                                                 ccTCCAACAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                   <  1:2157186/150‑1 (MQ=37)
                                                                     cAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa                                                                               >  2:1415500/1‑150 (MQ=37)
                                                                      aaaaGCAAAACCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                                              <  2:3926175/150‑1 (MQ=37)
                                                                        aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGGTTTTACGTAAATc                                                                            >  2:1019271/1‑150 (MQ=38)
                                                                           cAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTGTTTACGTAAATCAgg                                                                         >  2:1320699/1‑150 (MQ=38)
                                                                            aaataCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGt                                                                        <  2:2165976/150‑1 (MQ=17)
                                                                            aaaaaCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGGATAGCGTCGTGGACAGTCCTTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTGTTTACGTAAATCAGGt                                                                        >  1:3152632/1‑150 (MQ=255)
                                                                             aaaaCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTg                                                                       >  1:1435881/1‑150 (MQ=37)
                                                                                  ccgccgAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACt                                                                  >  2:2399658/1‑150 (MQ=37)
                                                                                  ccaccgAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACt                                                                  >  1:3483710/4‑150 (MQ=37)
                                                                                     ccgAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGAc                                                               >  2:2343947/1‑150 (MQ=37)
                                                                                     ccgAAGCGGGTTTTTACGTAAACCAGCTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAACTCAGGTGAAACTGAc                                                               <  2:3647087/150‑1 (MQ=255)
                                                                                      ttcggCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTCTAGCGTCGTGGCCAGTCATTCATCTTTCTGCCCCTCCCAAAGCCAAATCCCGCCGAAGCGGGGTGTTACGTAAATCAGGTGAAACTGAcc                                                              >  2:431725/5‑150 (MQ=11)
                                                                                        aaGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGa                                                            >  2:324186/1‑150 (MQ=37)
                                                                                           cGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATaa                                                         >  2:1682019/1‑150 (MQ=37)
                                                                                           cGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTTCGTACATCAGGTGAAACTGACCGATaa                                                         >  1:1971421/1‑150 (MQ=255)
                                                                                           cGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAAGTAAATCAGGTGAAAAAGACCGATaa                                                         >  2:2341007/1‑150 (MQ=255)
                                                                                           cGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGGTTTTACGTAAATCAGGTGAAACTGACCCATaa                                                         >  1:3347454/1‑150 (MQ=255)
                                                                                            gggTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACAGACCGATAAg                                                        >  1:1011164/1‑150 (MQ=37)
                                                                                            gggTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCCAAGCGGGTTTTTAAGTAAATAAGGTGAAACTGACCGATAAg                                                        >  1:1054548/1‑150 (MQ=255)
                                                                                            gggTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCACCAAAAGCAAAAACCCGCCGACGCGGGTTTTTACGTAAATCAGGTGAAACAGACCGATAAg                                                        >  1:1011187/1‑150 (MQ=255)
                                                                                              gTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGcc                                                      <  1:3553834/150‑1 (MQ=37)
                                                                                                ttttACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCGCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCgg                                                    >  1:3686665/1‑150 (MQ=38)
                                                                                                ttttACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTAATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAAGTAAATCAGGTGAAACTGACCGATAAGCCgg                                                    >  1:3187785/1‑150 (MQ=255)
                                                                                                ctttACGTCTATCAGGTGAAACTGTCAGATAAGCCGCTTTCTTTTGGGTATAGAGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCgg                                                    <  2:2201235/149‑1 (MQ=11)
                                                                                                  ttACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGt                                                  >  1:2244626/1‑150 (MQ=39)
                                                                                                  ttACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGt                                                  >  1:418962/1‑150 (MQ=39)
                                                                                                     cGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCt                                               >  2:2922861/1‑150 (MQ=255)
                                                                                                      gTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCCGCCCCTCCAAAAGCAAAAACCCGCCGAACCGGGTTTTTACGTAAAACAGGTGAAACTGACCGTTAAGCCGGGGTCTg                                              >  2:336854/1‑150 (MQ=11)
                                                                                                         catcAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGGAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCg                                           >  2:3359001/4‑150 (MQ=255)
                                                                                                         aaCCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTGATCTTTGTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTTAGTAAATCAGGTGAAAATGACCGATAATCCGGGTTCTGTCg                                           >  1:2403299/1‑150 (MQ=14)
                                                                                                           tcAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGGAACTGACCGATAAGCCGGGTTCTGTCGTg                                         <  2:239875/149‑1 (MQ=255)
                                                                                                           tcAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTg                                         >  2:3855430/2‑150 (MQ=255)
                                                                                                           tcAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTg                                         <  2:2124546/149‑1 (MQ=255)
                                                                                                           ccAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTg                                         >  2:124120/1‑150 (MQ=255)
                                                                                                           ccAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTg                                         <  1:733616/150‑1 (MQ=255)
                                                                                                               gTGAAACTGCCCGCTAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACTGTCATTCCTCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGGGTTCTGTCGTGGACa                                     <  2:2912599/150‑1 (MQ=14)
                                                                                                               gTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACa                                     <  2:787194/150‑1 (MQ=255)
                                                                                                                 gAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAgg                                   >  2:2654120/1‑149 (MQ=255)
                                                                                                                     cTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCAACTTTCTGCCCCTCCAAAATCAAAAAACCGCCCAAGCGGGTTTTTAAGTAAATAAGGTGAAAAATACCGATAAGCCGGGTTCTGTCGTGGACAGtcagt                               >  1:1396140/1‑148 (MQ=17)
                                                                                                                        aCCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGGCCGTCAATCAACTTTCAGCCCCACCAAAAACAAAAAACCCCCCAACCGGGTTTTTAAGTAAATAAGGTGAAAAAGAACGAAAAGCCCGGTTCTGTAGTGGACAGtcaataaa                            >  1:1774678/1‑145 (MQ=9)
                                                                                                                          cGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCt                          <  2:1435881/150‑1 (MQ=255)
                                                                                                                           gATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTa                         >  2:3284778/1‑150 (MQ=255)
                                                                                                                            aTAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAg                        <  2:2597377/150‑1 (MQ=255)
                                                                                                                              aaGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGc                      <  1:2943934/150‑1 (MQ=255)
                                                                                                                                gCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCa                    <  2:1734661/150‑1 (MQ=255)
                                                                                                                                 ccGCCTTCTTTTGGGTATAGAGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtcattcat                            >  2:2208947/1‑141 (MQ=255)
                                                                                                                                  cGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCAATCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTAATCTAGGCCAGc                  >  2:3429732/1‑150 (MQ=255)
                                                                                                                                   gCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAAGTAAATCAGGTGAAACAGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTATGCCAGCa                 >  1:3019786/1‑150 (MQ=255)
                                                                                                                                    ctttctttTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCaa                >  2:3741919/1‑150 (MQ=255)
                                                                                                                                     tttctttTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCCCCGAAGCGGGTTTTTACGTAAATAATGTGAAACAGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAAt               >  1:1657038/1‑150 (MQ=255)
                                                                                                                                     tttctttTGGGTATAGCGTCGTGGACAGTCATTCATATTTATGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACAGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCaaa               >  1:1445002/1‑149 (MQ=255)
                                                                                                                                        ctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGc            <  1:3855430/150‑1 (MQ=255)
                                                                                                                                        cttcTGGGTATAGTGTCGTGGAGAGTGATTCATCTGTCTGCCCCTCCAAAAGCAAAAACCGGCCGAAGCGGGTTTTTAGGTACATCAGGTGAATCGGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGc            <  2:2244992/150‑1 (MQ=18)
                                                                                                                                          tttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGt                                          <  1:2851770/118‑1 (MQ=255)
                                                                                                                                          tttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGt                                          >  2:2851770/1‑118 (MQ=255)
                                                                                                                                            tGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCAc        <  2:2707060/150‑1 (MQ=255)
                                                                                                                                               gTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTgg     <  1:1449017/150‑1 (MQ=255)
                                                                                                                                               gTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTgg     >  2:2275186/1‑150 (MQ=255)
                                                                                                                                               gTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGATTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTgg     <  1:1630277/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTC  >  NC_000913/3269994‑3270289

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: