Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,448,278 C→A 6.7% intergenic (‑129/+36) rplN ← / ← rpsQ 50S ribosomal subunit protein L14/30S ribosomal subunit protein S17

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,448,2780CA 6.7% 486.8 / 16.0 168intergenic (‑129/+36)rplN/rpsQ50S ribosomal subunit protein L14/30S ribosomal subunit protein S17
Reads supporting (aligned to +/- strand):  ref base C (81/75);  new base A (6/6);  total (87/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 170 total aligned reads displayed.

CATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTTCTCGT  >  NC_000913/3448129‑3448426
                                                                                                                                                     |                                                                                                                                                    
cATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTc                                                                                                                                                      <  1:673992/150‑1 (MQ=255)
cATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTc                                                                                                                                                      >  2:1530252/1‑150 (MQ=255)
 aTAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCt                                                                                                                                                     >  1:1799322/1‑150 (MQ=255)
    gTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAg                                                                                                                                                  >  1:624024/1‑150 (MQ=255)
       tGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCg                                                                                                                                               >  1:1840754/1‑150 (MQ=255)
        gTTCTTGGATCATTTTAGTGCTACGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGt                                                                                                                                              >  1:613489/1‑150 (MQ=255)
         ttCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGtt                                                                                                                                             <  2:867078/150‑1 (MQ=255)
          tCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACGCGAAAATCTGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACTCCGCTTCAAGGATATGTGTAGAAAAAATAAACGGCTCATTTCTGAGCCGttt                                                                                                                                            <  2:1584428/150‑1 (MQ=255)
          tCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGttt                                                                                                                                            >  2:1027980/1‑150 (MQ=255)
          tCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGttt                                                                                                                                            >  1:1125915/1‑150 (MQ=255)
          tCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCAGAAATGAGCCGttt                                                                                                                                            <  2:1139029/150‑1 (MQ=255)
           cTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTa                                                                                                                                           <  2:1041059/150‑1 (MQ=255)
            ttGGATCATTTTAGTGCTCCGCTAATTTCAACTACTACTGAGACCCGAAAATCAGGTCGTTATAAATCCCCATATCGAGGTCGCGTCATTATAACACCGCTTCAAGGATATGTGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTAt                                                                                                                                          <  1:1079772/150‑1 (MQ=255)
             tGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTAtt                                                                                                                                         <  2:1425794/150‑1 (MQ=255)
               gATCATTTTAGTGCTCCGCTAATGTCAACTAGTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTTGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCg                                                                                                                                       >  1:1811615/1‑150 (MQ=255)
               gATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGGATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGGCGTTTATTCg                                                                                                                                       >  2:1270894/1‑150 (MQ=255)
               gATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCg                                                                                                                                       <  2:1264793/150‑1 (MQ=255)
               gATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCg                                                                                                                                       >  2:745862/1‑150 (MQ=255)
               gATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCg                                                                                                                                       >  1:1565396/1‑150 (MQ=255)
               gATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCAGAAATGAGCCGTTTAtttt                                                                                                                                       <  2:617169/150‑3 (MQ=255)
                  cATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCAGAAATGAGCCGTTTAttttttc                                                                                                                                    <  2:1042634/150‑6 (MQ=255)
                   aTTTTAGTGCTCCTCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTAtt                                                                                                                                   <  2:1612116/150‑1 (MQ=255)
                    ttttAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTg                                                                                                                                  <  2:1837530/150‑1 (MQ=255)
                             gccGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTAAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTACAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAgtgt                                                                                                                         <  1:493562/149‑1 (MQ=255)
                              ccGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAATTCCCCATATCGAGGGCGAGGCATTATAACAACGCTGCAAGGATATGTTTAGAAAAAATAAACGGCTCATATCTTACCCGATTATTCGAATTGATAGAGTGTa                                                                                                                        >  1:342158/1‑150 (MQ=255)
                               cGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTAc                                                                                                                       <  1:1292560/150‑1 (MQ=255)
                                 cTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGATAGTGTACTg                                                                                                                     >  1:1007123/1‑150 (MQ=255)
                                     tGTCAACTACTACTGAGACACGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTAtt                                                                                                                 <  2:1762734/150‑1 (MQ=255)
                                      ttCAACTACTATTCAGACCCTAAAATCAGGTGGTTAAAAATCCCCATATCCCGGGCGCGGCATTTTAACACGGCTTCAAGTAATTGGGAATAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTAc                                                                                                               <  2:1395379/149‑1 (MQ=255)
                                      gTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTa                                                                                                                <  1:564727/150‑1 (MQ=255)
                                       tCAACTACTACTGAGACCCGAAAATCATGTAGTTAAAAATCCCCATCACGAGGGCGCGGCATTATAACACCGCTTCAAGGCTATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTAc                                                                                                               <  2:1487873/150‑1 (MQ=255)
                                       tCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTCAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTAc                                                                                                               <  1:1159987/150‑1 (MQ=255)
                                        cAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACa                                                                                                              <  1:75783/150‑1 (MQ=255)
                                           ctactactGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGACATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGaa                                                                                                           <  1:1622379/150‑1 (MQ=255)
                                            ttctaatGAGACACTAAATTCAGGTCGTTATAAATCCACATATCGAGGGCTCGGCATTATATTACCGCTTCAATGATATTGGTATAAAAAATAATCTGCTCATTTCTGAGCCGTTTATTCTTATTGAGAGAGTGTACTGTATTTCAGAAc                                                                                                          <  1:1250294/148‑1 (MQ=255)
                                            tactactGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAAc                                                                                                          <  2:1215059/150‑1 (MQ=255)
                                            tactactGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAAc                                                                                                          <  1:1129256/150‑1 (MQ=255)
                                              ctactGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCg                                                                                                        >  1:1590356/1‑150 (MQ=255)
                                              ctactGAGACCCGAAAATCAGGGGGTTAAAAATCCCCATTTCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCg                                                                                                        <  2:208111/150‑1 (MQ=255)
                                                 ctGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCTAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCtt                                                                                                     <  2:500004/150‑1 (MQ=255)
                                                 ctGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCtt                                                                                                     <  2:1549722/150‑1 (MQ=255)
                                                    agaCCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTtct                                                                                                  <  1:1727179/150‑1 (MQ=255)
                                                     gaCCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGTTATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTTTATTACAGAACCGCTTtctc                                                                                                 >  2:2004556/1‑150 (MQ=255)
                                                      aCCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGTATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTctct                                                                                                <  2:906655/150‑1 (MQ=255)
                                                       cccGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGCGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTa                                                                                               >  1:1991277/1‑150 (MQ=255)
                                                       cccGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGGAGAAAAAATAAACGGCTCATTTCTGAGACGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTa                                                                                               >  2:1214268/1‑150 (MQ=255)
                                                       cccGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTAAAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGCGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTAAAGAACCGCTTTCtata                                                                                               >  2:1427023/1‑150 (MQ=255)
                                                           aaaaTCAGGTCGTTAAAAATCCCCAAATCGAGGGCGCGGCATTATAACACCGCTTCAAGTATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAAc                                                                                           <  1:945653/150‑1 (MQ=255)
                                                                  ggTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAAcc                                                                                    <  2:250263/150‑1 (MQ=255)
                                                                     cGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGc                                                                                 >  1:1398273/1‑150 (MQ=255)
                                                                      gTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTAGTATTTAGAGAGTGTACTGTATTACAGAAACGCTTACTCTACAACGCGCACCAGCg                                                                                >  1:654299/1‑150 (MQ=255)
                                                                       ttAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAAACGCTTTATCTACAACGCGAACCAGCGt                                                                               >  2:731685/1‑150 (MQ=255)
                                                                                cccATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGAttt                                                                      <  2:1535905/150‑1 (MQ=255)
                                                                                cccATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGAttt                                                                      >  1:1687144/1‑150 (MQ=255)
                                                                                cccATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGAttt                                                                      >  2:548651/1‑150 (MQ=255)
                                                                                 ccATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTa                                                                     >  1:1583955/1‑150 (MQ=255)
                                                                                     atCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCt                                                                 >  1:549984/1‑150 (MQ=255)
                                                                                        gAGGGAGAGGCTTTATAAGACCGCTTCAAGTATATTGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGATCCAGCGTCCAGGATTTAGTCTTgg                                                              <  2:222747/150‑1 (MQ=255)
                                                                                            gcgcGGAATTATAACACCGCTTCAAGTATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAAGACTAAATCCTGGAcgc                                                          <  1:792343/150‑3 (MQ=255)
                                                                                             cgcgGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGc                                                         >  1:1529396/1‑150 (MQ=255)
                                                                                             cgcgGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGc                                                         >  1:814212/1‑150 (MQ=255)
                                                                                                 gCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTTGACAGCGGAc                                                     >  2:133744/1‑150 (MQ=255)
                                                                                                       taACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCAtt                                               >  2:1410788/1‑150 (MQ=255)
                                                                                                       taACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCACATTTATGATCCCTGAATTGGTATAGAGAGAGTGTACTGTATTAAACAGCCGCTTTCTATACAACGCGAACAACAGACCAGGCTTTAGTCTTGGTCTGCGGACCGCaat                                               >  1:344909/1‑148 (MQ=255)
                                                                                                          cacCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTcgc                                            >  1:1516351/1‑150 (MQ=255)
                                                                                                            ccGCTTCAAGGATATGGGTATAAAAAATAAACCGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCgg                                          <  2:1274695/150‑1 (MQ=255)
                                                                                                            ccGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCgg                                          >  2:1066719/1‑150 (MQ=255)
                                                                                                               cTTCAAGGATATTGGTAGAATAAATAATCTGCTCATTTATGAGTCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGCACCAGCGTCCAGGATTTAGTCTTTGACAGCGGACGGCATTCGCGGAtt                                       <  2:376690/150‑1 (MQ=255)
                                                                                                                  cAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCa                                    >  2:1878935/1‑150 (MQ=255)
                                                                                                                  cAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCa                                    <  2:1840754/150‑1 (MQ=255)
                                                                                                                   aaGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCaa                                   >  1:266080/1‑150 (MQ=255)
                                                                                                                   aaGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCaa                                   <  2:1444821/150‑1 (MQ=255)
                                                                                                                     ggATATGGGTAGAAAAAATAAACGGCTCATTTCGGAGCCGTTGATTCGGATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAAAGCGAACCAGCGTCCAGGAGTTAGTCTTGGACAGCGGACGGCATTCGCTGATTTCAAca                                 >  2:352887/1‑149 (MQ=255)
                                                                                                                       atatGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCAc                               >  1:824927/1‑150 (MQ=255)
                                                                                                                         atGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGt                             <  2:1590356/150‑1 (MQ=255)
                                                                                                                          tGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTc                            >  2:1682778/1‑150 (MQ=255)
                                                                                                                            ggTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCAc                          <  2:1457064/150‑1 (MQ=255)
                                                                                                                            ggTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTAAACCACGTCAc                          >  1:844764/1‑150 (MQ=255)
                                                                                                                             gTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCAcc                         <  2:1440024/150‑1 (MQ=255)
                                                                                                                               aGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGa                       >  2:801553/1‑150 (MQ=255)
                                                                                                                               aGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGa                       >  1:1017498/1‑150 (MQ=255)
                                                                                                                               aGAAAAAATAAACGGCTCATTTCTGAGACGTTTATTCGTAGTGAGAGAGTGTACTGTATTACAGAACCGCTTGCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTAGGACAGCGGACGGCAGTCGCGGATTTCACCCACATCACCGa                       >  1:1907666/1‑150 (MQ=255)
                                                                                                                                gAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGAt                      >  2:1946019/1‑150 (MQ=255)
                                                                                                                                 aaaaaaTAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATa                     >  1:1313673/1‑150 (MQ=255)
                                                                                                                                   aaaaTAAACGGCTCATTTCTGAGCCGTTGATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGGCCAGGAGTTAGTCTTGGACAGCGGCAGGCATTCGCGGATTTCAACCCCGTCACCGagccc                   >  1:1275740/1‑146 (MQ=255)
                                                                                                                                    aaaTAAACGGCTCATTTCTGAGCCGTTTATTCGTTTTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCg                  <  1:1479008/150‑1 (MQ=255)
                                                                                                                                        aaaCGGCTCAGAAATGAGCCGTTTATTAGTATTGAGAGAGTGTACTGTATTACAGAACCGCTATCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATGTCAACCACGTAACCGATACCGCAtt              >  1:1901745/1‑150 (MQ=255)
                                                                                                                                         aaCGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTc             >  1:432240/1‑150 (MQ=255)
                                                                                                                                         aaCGGCTCAGAAATGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTc             >  2:924691/1‑150 (MQ=255)
                                                                                                                                         aaCGGCTCAGAAATGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTc             >  2:793871/1‑150 (MQ=255)
                                                                                                                                          aCGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCGTTCTCTACAACGCGAACCAGCGTCCAGGAGTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCg            <  2:1539529/150‑1 (MQ=255)
                                                                                                                                           cGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCgt           >  1:1664694/1‑150 (MQ=255)
                                                                                                                                            ggCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCgtt          >  2:1666669/1‑150 (MQ=255)
                                                                                                                                             gCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCgttg         >  1:764115/1‑150 (MQ=255)
                                                                                                                                             gCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCgttg         <  2:1398414/150‑1 (MQ=255)
                                                                                                                                               tCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCTAACCAGCGTACAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCgttgtt       >  1:1641548/1‑150 (MQ=255)
                                                                                                                                                cATTTCTGAGCCGTTTTTTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTTCAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGAGGATTTCAACCACGTCACCGATACCGCATTCGTTGTtc      <  1:869423/150‑1 (MQ=255)
                                                                                                                                                 aTTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTtct     <  1:1420810/150‑1 (MQ=255)
                                                                                                                                                 aTTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTtct     <  2:1799322/150‑1 (MQ=255)
                                                                                                                                                    tCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTTCTCGt  >  2:1134304/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTTCTCGT  >  NC_000913/3448129‑3448426

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: