Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,448,275 T→G 5.5% intergenic (‑126/+39) 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,2750TG 5.5% 444.1 / 11.2 145intergenic (‑126/+39)rplN/rpsQ50S ribosomal subunit protein L14/30S ribosomal subunit protein S17
Reads supporting (aligned to +/- strand):  ref base T (74/62);  new base G (3/5);  total (77/68)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

Only 100 of 147 total aligned reads displayed.

AGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTTC  >  NC_000913/3448127‑3448422
                                                                                                                                                    |                                                                                                                                                   
aGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCAtt                                                                                                                                                    >  2:1341928/1‑150 (MQ=255)
  cATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTc                                                                                                                                                  <  2:941905/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)
           ttCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCAGAAATGAGCCGtt                                                                                                                                         <  2:1284693/150‑1 (MQ=255)
            tCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATTACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGttt                                                                                                                                        >  1:1504470/1‑150 (MQ=255)
            tCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCAGAAATGAGCCGttt                                                                                                                                        <  2:2099026/150‑1 (MQ=255)
            tCTTGGATCATTTTAGTGCTACGCTATTGTCAACTAATACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTGCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGttt                                                                                                                                        <  1:115246/150‑1 (MQ=255)
             cTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTa                                                                                                                                       <  2:1039212/150‑1 (MQ=255)
              ttGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCTAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCAGAAATGAGCCGTTTAt                                                                                                                                      <  1:776615/150‑1 (MQ=255)
                ggATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTc                                                                                                                                    >  1:59094/1‑150 (MQ=255)
                ggATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATACCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCAGAAATGAGCCGTTTAttt                                                                                                                                    <  1:197827/150‑2 (MQ=255)
                 gATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCg                                                                                                                                   >  1:1492226/1‑150 (MQ=255)
                    cATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTAt                                                                                                                                >  1:453076/1‑150 (MQ=255)
                    cATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGGTTATTCGTAt                                                                                                                                >  1:1387001/1‑150 (MQ=255)
                     catttAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTAtt                                                                                                                               <  2:1407898/148‑1 (MQ=255)
                     aTTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTAtt                                                                                                                               <  1:511272/150‑1 (MQ=255)
                     aTTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTAtt                                                                                                                               <  1:1219090/150‑1 (MQ=255)
                         tAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTgaga                                                                                                                           >  2:661619/1‑150 (MQ=255)
                            tGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGagagag                                                                                                                        >  2:2172251/1‑150 (MQ=255)
                              cTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAgtg                                                                                                                      >  2:1320251/1‑150 (MQ=255)
                              cTCCGCTAATGTCAACTACTACTGAGACACGAAAATCAGGTCGTTAAAAATCCCCATCTCGAGGGCGCGGCATTATAACCCCGCTTCACGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAgtg                                                                                                                      <  1:2047356/150‑1 (MQ=255)
                               tCCTCTAATGTCAACTACTACTTATACCCTAAAATCAGGTGGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAgtgt                                                                                                                     <  1:210240/150‑1 (MQ=255)
                                  gCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACt                                                                                                                  <  1:30031/150‑1 (MQ=255)
                                    tAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGt                                                                                                                <  2:608003/150‑1 (MQ=255)
                                    tAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGt                                                                                                                >  2:1225018/1‑150 (MQ=255)
                                    tAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGt                                                                                                                <  1:1149138/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)
                                       tGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCGTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTAtt                                                                                                             >  1:1703137/1‑150 (MQ=255)
                                        gTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTAGTATTGAGAGAGTGTACTGTATTa                                                                                                            >  2:1296348/1‑150 (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)
                                                ctactGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCg                                                                                                    <  1:817744/150‑1 (MQ=255)
                                                   ctGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACCCCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTTTTGAGAGAGTGTACTGTATTACAGAACCGCtt                                                                                                 >  2:897422/1‑150 (MQ=255)
                                                    tGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCttt                                                                                                >  1:243824/1‑150 (MQ=255)
                                                    tGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAATGTACTGTATTACAGAACCGCttt                                                                                                >  2:1863460/1‑150 (MQ=255)
                                                      agaCCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTtct                                                                                              <  1:1294556/150‑1 (MQ=255)
                                                         cccGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTa                                                                                           >  2:1070871/1‑150 (MQ=255)
                                                         cccGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTa                                                                                           >  2:402417/1‑150 (MQ=255)
                                                         cccGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTa                                                                                           >  2:17506/1‑150 (MQ=255)
                                                          ccGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACCCCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTAc                                                                                          >  1:12790/1‑150 (MQ=255)
                                                          ccGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTAc                                                                                          >  2:958507/1‑150 (MQ=255)
                                                          ccGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTAc                                                                                          >  1:2187169/1‑150 (MQ=255)
                                                          ccGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTAc                                                                                          >  1:2195346/1‑150 (MQ=255)
                                                           cGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACa                                                                                         >  1:32520/1‑150 (MQ=255)
                                                           cGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGAATTACAGAACCGCTTTCTCTACa                                                                                         >  2:1059077/1‑150 (MQ=255)
                                                               aaTCAGGTCTTAAAAAATCCCAATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATACACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAAcgc                                                                                     <  1:2066353/150‑1 (MQ=255)
                                                               aaTCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAAcgc                                                                                     <  2:415776/150‑1 (MQ=255)
                                                                    ggTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACACCGCGAAcc                                                                                >  1:222910/1‑150 (MQ=255)
                                                                    ggTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAAcc                                                                                <  2:793321/150‑1 (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)
                                                                                   ccATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTa                                                                 >  1:1591806/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)
                                                                                                     aTTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACgg                                               <  1:1092747/150‑1 (MQ=255)
                                                                                                     aTTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTAATGAGAGAGTGTACTGTATTACAGAAACGCTTTCTCTACAACGCCAACCAGCGTCCAGGAATTAGTCTTGTACAGCGTACgg                                               >  1:841868/1‑150 (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)
                                                                                                               cGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGGCAGCGGACGGCATTCGCGGa                                     <  1:395608/150‑1 (MQ=255)
                                                                                                                gCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGAt                                    >  1:21081/1‑150 (MQ=255)
                                                                                                                gCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCATGGACAGCGGACGGCATTCGCGGAt                                    >  2:148122/1‑150 (MQ=255)
                                                                                                                   tCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGAttta                                 >  1:1680847/1‑149 (MQ=255)
                                                                                                                    aaaGGATATTGTTAGAAAAAATAAACGGCTCATTTCTGATCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGATTTCTCTACAACTCTAAACAGCGTCCAGTATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCa                                <  2:376032/149‑1 (MQ=255)
                                                                                                                            tGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACAGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTc                        >  2:427967/1‑150 (MQ=255)
                                                                                                                             gggTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCa                       <  2:1246936/150‑1 (MQ=255)
                                                                                                                             gggTAGAAAAAATAAACGGCTAATTTCTGAGCCTTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAAAGCGAACCAGCGTCCAGGATTTAGACTAGGACAGCGGTAGGAATTCGCGGATTTCAACCACGTCa                       >  2:1742870/1‑150 (MQ=255)
                                                                                                                              ggTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCATGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCAc                      <  1:1320251/150‑1 (MQ=255)
                                                                                                                              ggTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCAc                      >  2:2020328/1‑150 (MQ=255)
                                                                                                                               gTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCAcc                     >  2:2135474/1‑150 (MQ=255)
                                                                                                                                tataaaaaaTAAACGTCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCg                    <  1:888783/147‑1 (MQ=255)
                                                                                                                                caatacgaaTAAACGGCTCAGAAATGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCACCCACGTCACCg                    >  2:47556/8‑150 (MQ=38)
                                                                                                                                 aGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTAGTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGa                   >  1:956789/1‑150 (MQ=255)
                                                                                                                                  gAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGAt                  >  2:1682743/1‑150 (MQ=255)
                                                                                                                                   aaaaaaTAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTAATGTATTACAGATCCGCTGTCTCTACAACGCGAAACAGCGTCCAGGATTTAGTCTTTGACAGCGGACGGCATTAGCGGATTTAAACCACGTCACCGATa                 >  1:375129/1‑150 (MQ=255)
                                                                                                                                        aTAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCa            >  1:1136911/1‑150 (MQ=255)
                                                                                                                                         tAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCACCCACGTCACCGATACCGCAt           >  1:88005/1‑150 (MQ=255)
                                                                                                                                           aaCGGCTCAGAAATGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTc         >  2:1079493/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)
                                                                                                                                                atCATTTCTGAGCCGTTTATTCGTATTGAGAGATTGTACTGTATTACAGAACCGCTATCTCTACAACTCGAAGCAGCGTCCAGGATTTAGTCTTGGAAAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCgtttt    <  2:1363522/149‑3 (MQ=255)
                                                                                                                                                  cATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTtc  <  1:1341928/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTTC  >  NC_000913/3448127‑3448422

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: