Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,448,277 T→A 5.9% intergenic (‑128/+37) 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,2770TA 5.9% 427.1 / 11.5 145intergenic (‑128/+37)rplN/rpsQ50S ribosomal subunit protein L14/30S ribosomal subunit protein S17
Reads supporting (aligned to +/- strand):  ref base T (73/63);  new base A (4/5);  total (77/68)
Fisher's exact test for biased strand distribution p-value = 7.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 147 total aligned reads displayed.

GCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTTCTCGT  >  NC_000913/3448128‑3448426
                                                                                                                                                     |                                                                                                                                                     
gCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCATCAAGGATATGGGTAGAAAAAATAAACGGCTCAttt                                                                                                                                                       >  1:2201900/1‑150 (MQ=255)
 cATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTc                                                                                                                                                      <  2:941905/150‑1 (MQ=255)
  aTAGTCTGTTATTGGATCATTTTAGTGCTCCGCTAATTTAAACTAGTACTTTGTCCCGAAAAGCAGGTCGTTTAAAATCCCCATCTCGAGGGCGCGGCATTATAAGACCGCTTCAAGGTTATGGGTAGAAAAACTAAACGGATCAATTCt                                                                                                                                                     <  2:1058026/150‑1 (MQ=255)
      tCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGc                                                                                                                                                 <  2:3607/150‑1 (MQ=255)
       cTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGcc                                                                                                                                                >  2:862178/1‑150 (MQ=255)
        tGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCg                                                                                                                                               <  2:712031/150‑1 (MQ=255)
         gTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGt                                                                                                                                              >  1:1201597/1‑150 (MQ=255)
         gTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAACCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGt                                                                                                                                              >  1:813519/1‑150 (MQ=255)
          ttCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCAGAAATGAGCCGtt                                                                                                                                             <  2:1284693/150‑1 (MQ=255)
          tacTTGGATCATTTTAGTGCTCCGCTAATGTCAACTCCTACTGACACCCGAAAATCAGGCCGTTAAAAATCCCCATATCGAGGGCGCGGCATGATAACCCCGCGTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGtt                                                                                                                                             <  1:292885/148‑1 (MQ=255)
           tCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATTACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGttt                                                                                                                                            >  1:1504470/1‑150 (MQ=255)
           tCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCAGAAATGAGCCGttt                                                                                                                                            <  2:2099026/150‑1 (MQ=255)
           tCTTGGATCATTTTAGTGCTACGCTATTGTCAACTAATACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTGCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGttt                                                                                                                                            <  1:115246/150‑1 (MQ=255)
              tGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTAtt                                                                                                                                         >  1:1395806/1‑150 (MQ=255)
               ggATCATTGTAGAGCTCAGCTAATGCCAACTTCTACTGAGACACGAAAATCAGGTCGTTAAACATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCAGAAATGAGCCGTTTAttt                                                                                                                                        <  1:369860/150‑2 (MQ=255)
                gATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCg                                                                                                                                       >  1:415776/1‑150 (MQ=255)
                   cATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCTTAt                                                                                                                                    >  1:163820/1‑150 (MQ=255)
                   cATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTAt                                                                                                                                    >  1:453076/1‑150 (MQ=255)
                    catttAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTAtt                                                                                                                                   <  2:1407898/148‑1 (MQ=255)
                        tAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTgaga                                                                                                                               >  2:661619/1‑150 (MQ=255)
                           tGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGagagag                                                                                                                            >  2:2172251/1‑150 (MQ=255)
                             cTCCGCTAATGTCAACTACTACTGAGACACGAAAATCAGGTCGTTAAAAATCCCCATCTCGAGGGCGCGGCATTATAACCCCGCTTCACGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAgtg                                                                                                                          <  1:2047356/150‑1 (MQ=255)
                              tCCTCTAATGTCAACTACTACTTATACCCTAAAATCAGGTGGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAgtgt                                                                                                                         <  1:210240/150‑1 (MQ=255)
                              tCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAgtgt                                                                                                                         >  2:1308142/1‑150 (MQ=255)
                                 gCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACt                                                                                                                      <  1:30031/150‑1 (MQ=255)
                                   tAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGt                                                                                                                    >  2:1225018/1‑150 (MQ=255)
                                   tAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGt                                                                                                                    <  1:1149138/150‑1 (MQ=255)
                                   tAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGt                                                                                                                    <  2:608003/150‑1 (MQ=255)
                                   tAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGCGCGCTGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGTCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGt                                                                                                                    >  1:1564036/1‑150 (MQ=255)
                                   tAATGTCAACTACTACTGAGAACCGAAAATCATGTCGTTAAAAATCCCCATATAGAGGTCGCGGCATTATAACACCGATTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGt                                                                                                                    <  1:59051/150‑1 (MQ=255)
                                   aaaTTTCAACTACTACTCAGACCCCAAAATCAGGTCGTTAAATATCCCCATATCGATGGCGCGGCATTATAACACCGCTTAAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGt                                                                                                                    <  2:2208956/149‑1 (MQ=255)
                                       gTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTAGTATTGAGAGAGTGTACTGTATTa                                                                                                                >  2:2022452/1‑150 (MQ=255)
                                        tCAAATACTACTGATATCCGAAAATCAGGTCGATAAAAATGCTCATATCGAGGGCGCGGAATTATACCACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTAc                                                                                                               <  1:1927468/150‑1 (MQ=255)
                                         cAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGTGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTAATTTCTGAGCCGTTTATTCGTATTTAGAGATTGTACTGTATTACa                                                                                                              >  1:376032/1‑150 (MQ=255)
                                         cAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACa                                                                                                              >  2:1825993/1‑150 (MQ=255)
                                         cAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACa                                                                                                              <  2:2171135/150‑1 (MQ=255)
                                          tactactactGAGTCCAGAAAATCAGGTCGTTAAAAATCCCCATATCGATTGCGCGGCATTTTAACAGCGCTTCAAGGTTATGGGTAGAAAAAATAAACGGCTCATTTGTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAg                                                                                                             <  1:91228/149‑1 (MQ=255)
                                           actactactGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTATAAAAAATAAACGCCTCATTTCTGAGCCGTTTATTCGTATTGAGAGATTGTACTTTATTACAGa                                                                                                            >  2:1211867/1‑150 (MQ=255)
                                               ctactGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCg                                                                                                        <  1:817744/150‑1 (MQ=255)
                                                  ctGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACCCCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTTTTGAGAGAGTGTACTGTATTACAGAACCGCtt                                                                                                     >  2:897422/1‑150 (MQ=255)
                                                   tGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAATGTACTGTATTACAGAACCGCttt                                                                                                    >  2:1863460/1‑150 (MQ=255)
                                                     agaCCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTtct                                                                                                  <  1:1294556/150‑1 (MQ=255)
                                                        cccGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTa                                                                                               >  2:17506/1‑150 (MQ=255)
                                                        cccGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTa                                                                                               >  2:402417/1‑150 (MQ=255)
                                                        cccGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTa                                                                                               >  2:1070871/1‑150 (MQ=255)
                                                         ccGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTAc                                                                                              >  1:2187169/1‑150 (MQ=255)
                                                         ccGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTAc                                                                                              >  2:958507/1‑150 (MQ=255)
                                                          cGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGAATTACAGAACCGCTTTCTCTACa                                                                                             >  2:1059077/1‑150 (MQ=255)
                                                           gAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACaa                                                                                            >  1:1912184/1‑150 (MQ=255)
                                                              aaTCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAAcgc                                                                                         <  2:415776/150‑1 (MQ=255)
                                                                   ggTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACACCGCGAAcc                                                                                    >  1:222910/1‑150 (MQ=255)
                                                                   ggTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAAcc                                                                                    <  2:793321/150‑1 (MQ=255)
                                                                     tCGTTAAAAATCCCCATATCGAGGGCGCCGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTGTGAGACGTTTATTCGTATTGAGATAGTGTACTGTATTACAGCACCGCTTTCTCTACAACGCGAACCAg                                                                                  >  1:897999/1‑150 (MQ=255)
                                                                      cGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGc                                                                                 >  2:668035/1‑150 (MQ=255)
                                                                       gTTAAAAATCCCCATATCGCGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCTTTTATTCGTATTGAGAGAGTGTACTGTATTACAGCACCGCTTTCTCTACAACGCGAACCAGCg                                                                                >  2:1308893/1‑150 (MQ=255)
                                                                        ttAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGt                                                                               >  2:1037488/1‑150 (MQ=255)
                                                                          aaaaaTCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTcc                                                                             <  1:1277364/150‑1 (MQ=255)
                                                                            aaaTCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAg                                                                           <  1:1239972/150‑1 (MQ=255)
                                                                                   cATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAATGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAg                                                                    <  1:132723/150‑1 (MQ=255)
                                                                                   cATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAg                                                                    <  1:973141/150‑1 (MQ=255)
                                                                                   cATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAg                                                                    >  1:2043326/1‑150 (MQ=255)
                                                                                   cATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGCGTAGAAAAAATAAACGGCTCATTTCTGCGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAg                                                                    <  2:1435814/150‑1 (MQ=255)
                                                                                       tCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCtt                                                                <  2:76348/150‑1 (MQ=255)
                                                                                         gAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGGTTATTCGTATTGAGAGAGTGTACTGTTTTACAGAACGGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTgg                                                              >  2:178304/1‑150 (MQ=255)
                                                                                         gAGGGCGCGGCATTAAAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTgg                                                              <  2:666545/150‑1 (MQ=255)
                                                                                                 ggCATTATAACACCGCTTCAAGTATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAAGACTAAATCCTGGAcgctggt                                                      <  2:1556829/150‑7 (MQ=255)
                                                                                                  gCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTGTCTCTACAAAGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGAc                                                     >  2:1251953/1‑150 (MQ=255)
                                                                                                    aTTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACgg                                                   <  1:1092747/150‑1 (MQ=255)
                                                                                                     ttATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGc                                                  <  2:1757994/150‑1 (MQ=255)
                                                                                                        taACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCAtt                                               <  1:72874/150‑1 (MQ=255)
                                                                                                         aaCACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTc                                              >  2:1547469/1‑150 (MQ=255)
                                                                                                          acacCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCTTAATGAGAAAGTGAAATGTAATACATAACAGCATTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTcg                                             >  2:498325/1‑150 (MQ=255)
                                                                                                          acacCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTcg                                             >  2:684481/1‑150 (MQ=255)
                                                                                                            acCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGCGCGAGTGTCCTGTATTACAGAACCGCGTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTcgcg                                           <  1:1439515/150‑1 (MQ=255)
                                                                                                               gCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCATGGACAGCGGACGGCATTCGCGGAt                                        >  2:148122/1‑150 (MQ=255)
                                                                                                                cTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGAGAACAAGCGTCCAGGATTTAGTAATGAACAGCGGACGGCATTCGCGGAtt                                       >  1:1495736/1‑150 (MQ=255)
                                                                                                                 ttCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGAttt                                      <  1:209597/150‑1 (MQ=255)
                                                                                                                  tCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGAttta                                     >  1:1680847/1‑149 (MQ=255)
                                                                                                                  tCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTc                                     >  1:1285038/1‑150 (MQ=255)
                                                                                                                   aaaGGATATTGTTAGAAAAAATAAACGGCTCATTTCTGATCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGATTTCTCTACAACTCTAAACAGCGTCCAGTATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCa                                    <  2:376032/149‑1 (MQ=255)
                                                                                                                           tGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACAGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTc                            >  2:427967/1‑150 (MQ=255)
                                                                                                                            gggTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCa                           <  1:2172251/150‑1 (MQ=255)
                                                                                                                            gggTAGAAAAAATAAACGGCTAATTTCTGAGCCTTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAAAGCGAACCAGCGTCCAGGATTTAGACTAGGACAGCGGTAGGAATTCGCGGATTTCAACCACGTCa                           >  2:1742870/1‑150 (MQ=255)
                                                                                                                             ggTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCAc                          >  2:2020328/1‑150 (MQ=255)
                                                                                                                              gTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCAcc                         >  2:2135474/1‑150 (MQ=255)
                                                                                                                               caatacgaaTAAACGGCTCAGAAATGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCACCCACGTCACCg                        >  2:47556/8‑150 (MQ=38)
                                                                                                                                 gAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGAt                      >  2:1682743/1‑150 (MQ=255)
                                                                                                                                  aaaaaaTAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTAATGTATTACAGATCCGCTGTCTCTACAACGCGAAACAGCGTCCAGGATTTAGTCTTTGACAGCGGACGGCATTAGCGGATTTAAACCACGTCACCGATa                     >  1:375129/1‑150 (MQ=255)
                                                                                                                                   aaaaaTAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGATAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTCGACAGCGGACTGCATTCGCGGATTTCAACCACGTCACCGATAc                    <  1:582475/150‑1 (MQ=255)
                                                                                                                                   aaaaaTAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCCGAGTCCAGGATTTAGTATTGGACAGCGGACGGCATTCGCGGATTTCACCCACGGCACCGATAc                    >  2:1326447/1‑150 (MQ=255)
                                                                                                                                       aTAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCa                >  1:1136911/1‑150 (MQ=255)
                                                                                                                                        tAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCACCCACGTCACCGATACCGCAt               >  1:88005/1‑150 (MQ=255)
                                                                                                                                           aCGGGTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGACCCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTGTTGTACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCAttag            <  2:148588/150‑3 (MQ=255)
                                                                                                                                             ggCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCgtt          >  1:1075933/1‑150 (MQ=255)
                                                                                                                                               cTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCgttgt        <  2:243824/150‑1 (MQ=255)
                                                                                                                                               cTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCgttgt        <  1:862178/150‑1 (MQ=255)
                                                                                                                                               atCATTTCTGAGCCGTTTATTCGTATTGAGAGATTGTACTGTATTACAGAACCGCTATCTCTACAACTCGAAGCAGCGTCCAGGATTTAGTCTTGGAAAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCgtttt        <  2:1363522/149‑3 (MQ=255)
                                                                                                                                                 cATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTtc      <  1:1341928/150‑1 (MQ=255)
                                                                                                                                                 cATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTtc      <  1:684481/150‑1 (MQ=255)
                                                                                                                                                     tCTGAGCCGTTTATTCGTATTGAGAGAGTGGACTGTATTACTGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGTTTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTTCTCGt  <  1:1211867/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTTCTCGT  >  NC_000913/3448128‑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: