Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,448,275 T→G 5.3% 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.3% 490.1 / 12.1 167intergenic (‑126/+39)rplN/rpsQ50S ribosomal subunit protein L14/30S ribosomal subunit protein S17
Reads supporting (aligned to +/- strand):  ref base T (83/75);  new base G (4/5);  total (87/80)
Fisher's exact test for biased strand distribution p-value = 7.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

Only 100 of 168 total aligned reads displayed.

CATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTTCTC  >  NC_000913/3448129‑3448424
                                                                                                                                                  |                                                                                                                                                     
cATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTc                                                                                                                                                    >  2:1530252/1‑150 (MQ=255)
cATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTc                                                                                                                                                    <  1:673992/150‑1 (MQ=255)
 aTAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCt                                                                                                                                                   >  1:1799322/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                                                                                                                                          >  1:1125915/1‑150 (MQ=255)
          tCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGttt                                                                                                                                          >  2:1027980/1‑150 (MQ=255)
           cTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTa                                                                                                                                         <  2:1041059/150‑1 (MQ=255)
           cTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTa                                                                                                                                         <  2:1538543/150‑1 (MQ=255)
           cTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGGCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTa                                                                                                                                         <  2:327585/150‑1 (MQ=255)
            ttGGATCATTTTAGTGCTCCGCTAATTTCAACTACTACTGAGACCCGAAAATCAGGTCGTTATAAATCCCCATATCGAGGTCGCGTCATTATAACACCGCTTCAAGGATATGTGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTAt                                                                                                                                        <  1:1079772/150‑1 (MQ=255)
            ttGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGATACCCGAAAATCAGGTCTTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTGCACGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTAt                                                                                                                                        <  1:1514112/150‑1 (MQ=255)
            ttGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCTAAAATCCGGTCGTTAAAAATCCCCATATCGCGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATCTCTGAGCCGTTTAt                                                                                                                                        <  1:1950280/150‑1 (MQ=255)
            ttGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTAt                                                                                                                                        <  1:1341483/150‑1 (MQ=255)
             tGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTAtt                                                                                                                                       <  2:1425794/150‑1 (MQ=255)
              ggATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTc                                                                                                                                      >  2:1593742/1‑150 (MQ=255)
               gATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGGATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGGCGTTTATTCg                                                                                                                                     >  2:1270894/1‑150 (MQ=255)
               gATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCg                                                                                                                                     >  2:1544641/1‑150 (MQ=255)
               gATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCg                                                                                                                                     <  2:1271654/150‑1 (MQ=255)
               gATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCg                                                                                                                                     >  2:745862/1‑150 (MQ=255)
               gATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCAGAAATGAGCCGTTTAtttt                                                                                                                                     <  2:617169/150‑3 (MQ=255)
               gATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGCGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACCGCAGATTTCTGCACCGTTTAttag                                                                                                                                     >  2:1737742/1‑148 (MQ=255)
                  cATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCAGAAATGAGCCGTTTAttttttc                                                                                                                                  <  2:1042634/150‑6 (MQ=255)
                    ttttAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTg                                                                                                                                <  2:1837530/150‑1 (MQ=255)
                    ttttAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTg                                                                                                                                >  2:1420810/1‑150 (MQ=255)
                         gTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTgagaga                                                                                                                           >  2:1443611/1‑150 (MQ=255)
                          tGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGagagag                                                                                                                          <  1:400430/150‑1 (MQ=255)
                             gccGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTAAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTACAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAgtgt                                                                                                                       <  1:493562/149‑1 (MQ=255)
                              ccGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAATTCCCCATATCGAGGGCGAGGCATTATAACAACGCTGCAAGGATATGTTTAGAAAAAATAAACGGCTCATATCTTACCCGATTATTCGAATTGATAGAGTGTa                                                                                                                      >  1:342158/1‑150 (MQ=255)
                                 cTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGATAGTGTACTg                                                                                                                   >  1:1007123/1‑150 (MQ=255)
                                   aaTGTCTACTACTACTGATACCACAAAATCTGGTCTTTAAAAATCCCCATAACGAGGGCGCGGCATTATAACACCGCTTCATGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTa                                                                                                                 <  1:1698935/150‑1 (MQ=255)
                                      ttCAACTACTATTCAGACCCTAAAATCAGGTGGTTAAAAATCCCCATATCCCGGGCGCGGCATTTTAACACGGCTTCAAGTAATTGGGAATAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTAc                                                                                                             <  2:1395379/149‑1 (MQ=255)
                                      gTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTa                                                                                                              <  1:564727/150‑1 (MQ=255)
                                      gTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTa                                                                                                              >  1:1613156/1‑150 (MQ=255)
                                       tCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTCAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTAc                                                                                                             <  1:1159987/150‑1 (MQ=255)
                                        cAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACa                                                                                                            <  1:75783/150‑1 (MQ=255)
                                            ttctaatGAGACACTAAATTCAGGTCGTTATAAATCCACATATCGAGGGCTCGGCATTATATTACCGCTTCAATGATATTGGTATAAAAAATAATCTGCTCATTTCTGAGCCGTTTATTCTTATTGAGAGAGTGTACTGTATTTCAGAAc                                                                                                        <  1:1250294/148‑1 (MQ=255)
                                            tactactGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAAc                                                                                                        <  1:1129256/150‑1 (MQ=255)
                                              ctactGAGACCCGAAAATCAGGGGGTTAAAAATCCCCATTTCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCg                                                                                                      <  2:208111/150‑1 (MQ=255)
                                                 ctGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCTAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCtt                                                                                                   <  2:500004/150‑1 (MQ=255)
                                                 ctGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCtt                                                                                                   <  2:1549722/150‑1 (MQ=255)
                                                 ctGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCtt                                                                                                   >  2:1214296/1‑150 (MQ=255)
                                                 atgataaccGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCTTTATAACCCCGCTTCAAGGATATGGGTAGAAAAATTAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCtt                                                                                                   <  2:1520959/143‑1 (MQ=255)
                                                   gagaCCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCttta                                                                                                 >  2:147624/1‑149 (MQ=255)
                                                   gagaCCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTCTATTACAGCACCGCtttg                                                                                                 >  2:1282414/1‑149 (MQ=255)
                                                    agaCCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTtct                                                                                                <  1:1727179/150‑1 (MQ=255)
                                                      aCCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGTATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTctct                                                                                              <  2:906655/150‑1 (MQ=255)
                                                      aCCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTctct                                                                                              >  1:1578130/1‑150 (MQ=255)
                                                       cccGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGGAGAAAAAATAAACGGCTCATTTCTGAGACGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTa                                                                                             >  2:1214268/1‑150 (MQ=255)
                                                        ccGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTAc                                                                                            >  2:1429219/1‑150 (MQ=255)
                                                           aaaaTCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAAc                                                                                         <  1:1632750/150‑1 (MQ=255)
                                                           aaaaTCAGGTCGTTAAAAATCCCCAAATCGAGGGCGCGGCATTATAACACCGCTTCAAGTATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAAc                                                                                         <  1:945653/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)
                                                                         aaaacaccccATATCGCGCGCGCGGCACTATAACACCGCTCCAAGGATAGGGGTACAAAAAATAAACGGCTCATTTCTGAGCCGATTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTcc                                                                           <  2:1867792/144‑1 (MQ=255)
                                                                            aaTCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATACACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAgg                                                                        <  2:1643203/150‑1 (MQ=255)
                                                                                cccATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGAttt                                                                    >  2:548651/1‑150 (MQ=255)
                                                                                 ccATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTa                                                                   >  1:1583955/1‑150 (MQ=255)
                                                                                   atatCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGt                                                                 >  2:1225937/1‑150 (MQ=255)
                                                                                        gAGGGAGAGGCTTTATAAGACCGCTTCAAGTATATTGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGATCCAGCGTCCAGGATTTAGTCTTgg                                                            <  2:222747/150‑1 (MQ=255)
                                                                                           ggCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAAGACTAAATCCTGGAcg                                                         <  2:1796907/150‑2 (MQ=255)
                                                                                            gcgcGGAATTATAACACCGCTTCAAGTATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAAGACTAAATCCTGGAcgc                                                        <  1:792343/150‑3 (MQ=255)
                                                                                             cgcgGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGc                                                       >  1:814212/1‑150 (MQ=255)
                                                                                               cgGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCgg                                                     >  2:1135720/1‑150 (MQ=255)
                                                                                               cgGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAAGACTAAATCCTGGAcgctgg                                                     <  1:1544641/150‑6 (MQ=255)
                                                                                                       taACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCAtt                                             >  2:1403858/1‑150 (MQ=255)
                                                                                                          cacCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTcgc                                          >  1:1516351/1‑150 (MQ=255)
                                                                                                            ccGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCgg                                        >  2:1066719/1‑150 (MQ=255)
                                                                                                               cTTCAAGGATATTGGTAGAATAAATAATCTGCTCATTTATGAGTCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGCACCAGCGTCCAGGATTTAGTCTTTGACAGCGGACGGCATTCGCGGAtt                                     <  2:376690/150‑1 (MQ=255)
                                                                                                                  cAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCa                                  >  2:1878935/1‑150 (MQ=255)
                                                                                                                   aaGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCaa                                 <  2:1444821/150‑1 (MQ=255)
                                                                                                                   aaGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCaa                                 >  1:266080/1‑150 (MQ=255)
                                                                                                                     ggATATGGGTAGAAAAAATAAACGGCTCATTTCGGAGCCGTTGATTCGGATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAAAGCGAACCAGCGTCCAGGAGTTAGTCTTGGACAGCGGACGGCATTCGCTGATTTCAAca                               >  2:352887/1‑149 (MQ=255)
                                                                                                                       atatGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCAc                             >  1:824927/1‑150 (MQ=255)
                                                                                                                        tatGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACg                            <  1:1846461/150‑1 (MQ=255)
                                                                                                                          tGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTc                          >  1:626486/1‑150 (MQ=255)
                                                                                                                            ggTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCAc                        >  2:2008151/1‑150 (MQ=255)
                                                                                                                            ggTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTAAACCACGTCAc                        >  1:844764/1‑150 (MQ=255)
                                                                                                                              tAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCg                      <  1:1821553/150‑1 (MQ=255)
                                                                                                                               aGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGa                     >  2:801553/1‑150 (MQ=255)
                                                                                                                               aGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGa                     >  1:1017498/1‑150 (MQ=255)
                                                                                                                                gAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGAt                    >  2:1946019/1‑150 (MQ=255)
                                                                                                                                gAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAAAGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTAAACCACGTCCCCGAt                    >  2:1360727/1‑150 (MQ=255)
                                                                                                                                   aaaaTAAACGGCTCATTTCTGAGCCGTTGATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGGCCAGGAGTTAGTCTTGGACAGCGGCAGGCATTCGCGGATTTCAACCCCGTCACCGagccc                 >  1:1275740/1‑146 (MQ=255)
                                                                                                                                    aaaTAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCg                <  1:235905/150‑1 (MQ=255)
                                                                                                                                       tAAACGGCTCATTTCTGAGCCGTGTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCAt             >  1:1310526/1‑150 (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)
                                                                                                                                           cGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCgt         <  1:56237/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)
                                                                                                                                               tCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCTAACCAGCGTACAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCgttgtt     >  1:1641548/1‑150 (MQ=255)
                                                                                                                                                cATTTCTGAGCCGTTTTTTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTTCAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGAGGATTTCAACCACGTCACCGATACCGCATTCGTTGTtc    <  1:869423/150‑1 (MQ=255)
                                                                                                                                                 aTTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTtct   <  2:1799322/150‑1 (MQ=255)
                                                                                                                                                  tttCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTtctc  <  1:147624/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTTCTC  >  NC_000913/3448129‑3448424

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: