Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,653,698 A→T 5.9% intergenic (+359/‑157) sseA → / → ryfA 3‑mercaptopyruvate sulfurtransferase/novel sRNA, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,653,6980AT 5.9% 1254.0 / 44.7 371intergenic (+359/‑157)sseA/ryfA3‑mercaptopyruvate sulfurtransferase/novel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base A (181/162);  new base T (10/13);  total (192/179)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 375 total aligned reads displayed.

TGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGTTAGTGTATT  >  NC_000913/2653549‑2653847
                                                                                                                                                     |                                                                                                                                                     
  ttCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAgg                                                                                                                                                     >  1:4198504/1‑150 (MQ=255)
  ttCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAgg                                                                                                                                                     >  1:3381089/1‑150 (MQ=255)
  ttCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCACAAAATTAAAAATAATAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAATAAAACGCCAAAATCAgc                                                                                                                                                     >  1:4926688/1‑149 (MQ=255)
    cGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCg                                                                                                                                                   >  1:4831891/1‑150 (MQ=255)
    cGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCg                                                                                                                                                   <  2:4270952/150‑1 (MQ=255)
       taataaAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGtt                                                                                                                                                >  1:2054533/1‑150 (MQ=255)
         ataaaataaCCTGTCTTATTAGTCCGAATTCAGACACATATACATAAATCCTGCTCAAAATTCAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCTGATTTTGGCGGtttt                                                                                                                                              <  2:3928817/150‑1 (MQ=255)
         ataaaataaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGtttt                                                                                                                                              <  1:5494600/150‑1 (MQ=255)
          taaaataaCCTGGATTATTAGTCCGAATTCAGACAAATATAAAGAACTCCTGCTCAAATTTAATAATTCTATCCGGTAATATATATTACTTAATCATGTAAATTCTCTTTCCTTTATAAATCATACTACCGCCTGATTTTGGCGGttttt                                                                                                                                             <  2:4116265/150‑1 (MQ=255)
           aaaataaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCTGATTTTGGCGGtttttt                                                                                                                                            <  2:1692691/150‑1 (MQ=255)
           aaaataaCCTGGATTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGtttttt                                                                                                                                            <  2:6070984/150‑1 (MQ=255)
             aataaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTttgt                                                                                                                                          >  2:2051394/1‑150 (MQ=255)
             aataaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTGTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGttgtttgt                                                                                                                                          >  1:6697132/1‑150 (MQ=255)
              ataaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTttgtt                                                                                                                                         <  1:4042249/150‑1 (MQ=255)
               taaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTttgttg                                                                                                                                        >  1:4586290/1‑150 (MQ=255)
                caCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGc                                                                                                                                       <  2:2168256/149‑1 (MQ=255)
                aaCCTGGCTTATTAGTCCTAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGc                                                                                                                                       <  1:5504943/150‑1 (MQ=255)
                    tGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATACTGCTCAAAATTAAAAAATCTAACCGGTAAAATATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAAACGCCAAAATCAGGCGGTTTTTTGTTGATGGt                                                                                                                                   >  1:5998865/1‑150 (MQ=255)
                     ggCTTATTAGTCCTAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTc                                                                                                                                  <  2:6470734/150‑1 (MQ=255)
                        ttattaGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCCAAATCAGGCGGTTTTTTGTTGCTGGTCCgg                                                                                                                               >  2:3915046/1‑150 (MQ=255)
                             aGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTcgc                                                                                                                          <  2:3504072/150‑1 (MQ=255)
                               tCCGAATTCAGACAGATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCgg                                                                                                                        >  2:3124583/1‑150 (MQ=255)
                                ccGAACTCAGACAAAGATAAATAAATCCTGCACAAAATTAAAAATTCTAACCGGTAACAGATATTACTTAAACACCGAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGATTTTTGTTGCTGGTCCGGTTCGCGGc                                                                                                                       <  1:6393660/150‑1 (MQ=255)
                                  gAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCt                                                                                                                     >  1:5180898/1‑150 (MQ=255)
                                      tCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTcc                                                                                                                 <  1:1285777/150‑1 (MQ=255)
                                        aGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCcag                                                                                                               >  1:2008547/1‑150 (MQ=255)
                                              atatAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTGCGCGGCCTTTCCAGCAGGtt                                                                                                         >  2:4486563/1‑150 (MQ=255)
                                                ataaataaATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGt                                                                                                       <  2:2117897/150‑1 (MQ=255)
                                                 taaataaaTCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTa                                                                                                      <  2:900738/150‑1 (MQ=255)
                                                    ataaatCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAATTTCACTTTCGTTTAAAAAACAAAAATCCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTGCGGTTGGCGGCCTTTCCAGCAGGTTGTATTa                                                                                                   >  1:620291/1‑150 (MQ=255)
                                                    ataaatCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTa                                                                                                   <  2:4908198/150‑1 (MQ=255)
                                                     taaatCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTAc                                                                                                  >  2:1524488/1‑150 (MQ=255)
                                                      aaatCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTAcc                                                                                                 <  1:2678428/150‑1 (MQ=255)
                                                          ccTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCgtag                                                                                             >  1:1618956/1‑150 (MQ=255)
                                                          ccTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCgtag                                                                                             >  1:1954671/1‑150 (MQ=255)
                                                          ccTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCgtag                                                                                             >  2:4798545/1‑150 (MQ=255)
                                                            tGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCgtagta                                                                                           >  2:5525977/1‑150 (MQ=255)
                                                              cTCAAAATTAAAAATTCTAACCGGGAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAAt                                                                                         <  2:4694583/150‑1 (MQ=255)
                                                               tCAAAATTAAAAATTCTCACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTGTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATg                                                                                        >  2:5047020/1‑150 (MQ=255)
                                                               tCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTcc                                                                                                                                 >  1:863973/1‑109 (MQ=255)
                                                               tCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATg                                                                                        >  1:5127145/1‑150 (MQ=255)
                                                               tCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATg                                                                                        >  1:1681786/1‑150 (MQ=255)
                                                                    aTTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAgc                                                                                   <  2:4792005/150‑1 (MQ=255)
                                                                      tAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAgcgc                                                                                 <  2:1076814/150‑1 (MQ=255)
                                                                       aaaaaTTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTAACGTAGTAATGCAAGcgcg                                                                                >  1:5889369/1‑150 (MQ=255)
                                                                          aaTTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGtct                                                                             >  2:5628579/1‑150 (MQ=255)
                                                                            ttCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCa                                                                           >  1:5350399/1‑150 (MQ=255)
                                                                              cTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGc                                                                         <  1:1041445/150‑1 (MQ=255)
                                                                               tAACCGGTAAAAGATATTACTTAAACCTGTAAATTCAATTTCATTTAATAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTTTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCg                                                                        <  1:4152779/150‑1 (MQ=39)
                                                                               tAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCg                                                                        >  2:1995493/1‑150 (MQ=255)
                                                                               tAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCg                                                                        >  2:1899376/1‑150 (MQ=255)
                                                                                 aCCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGga                                                                      <  1:1524488/150‑1 (MQ=255)
                                                                                   cGGTCACAGATCTTACTGACACATGTAACTTCACTTTCCTTTAAACAACACACCACCGCCAACCTCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGgaga                                                                    <  2:6697132/150‑1 (MQ=37)
                                                                                    ggTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGAc                                                                   >  1:5649113/1‑150 (MQ=255)
                                                                                    ggTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGAc                                                                   >  2:2639480/1‑150 (MQ=255)
                                                                                      tAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACaa                                                                 >  1:2401246/1‑150 (MQ=255)
                                                                                        aaaGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATa                                                               >  1:4710751/1‑150 (MQ=255)
                                                                                        aaaGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATa                                                               <  1:3486896/150‑1 (MQ=255)
                                                                                          aGATATTACTTATTCTTGTAAATTCACTTTCATTTATAAAACAACATACCGGCAAAATCAGGCGGTTTTTTGTTTCTGGTCGGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACt                                                             <  2:4378914/150‑1 (MQ=255)
                                                                                              atTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTCGTAATGCAAGCGCGTCTCTGCGGAGACAATACTCGcc                                                         >  1:3596443/1‑150 (MQ=255)
                                                                                               ttACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCa                                                        <  1:3031641/150‑1 (MQ=255)
                                                                                                  cTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTa                                                     >  1:1968049/1‑150 (MQ=255)
                                                                                                  attAAACATGTAAATTCACTTTCCTTTAAAAAACAAAACACCGCCAAAATCTGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTa                                                     <  1:6613293/149‑1 (MQ=255)
                                                                                                   ttAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTaa                                                    <  1:3521459/150‑1 (MQ=255)
                                                                                                         aTGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAActctct                                              >  2:3704695/1‑150 (MQ=255)
                                                                                                          tGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAAt                                                                >  1:6482572/1‑131 (MQ=255)
                                                                                                          tGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCtt                                             >  1:1431651/1‑150 (MQ=255)
                                                                                                             aaaTTCACTTTCCTATAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCttttt                                          <  2:5166044/150‑1 (MQ=255)
                                                                                                              aaTTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTg                                         <  2:3466354/150‑1 (MQ=255)
                                                                                                               aTTCACTTTCCTTTAAAATTCAATATACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTGGGGTTCGCGGGCTTTGGAGGAGGTTGTATTTGGTTTGTTATGGTAGCGCGTCTCTGCGGAGAGAATAGTCGCGAGTAAGTCTCTTTTTGt                                        >  2:1385052/1‑150 (MQ=255)
                                                                                                                ttCACTTTCCTTTAAAAAACAATAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTGCGCGGCCTTTCCAGCAGGTTGTATTACCGGAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTc                                       >  1:3818729/1‑150 (MQ=255)
                                                                                                                ttCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGGAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTc                                       >  2:3249512/1‑150 (MQ=255)
                                                                                                                  gaCTTTCCTTAATAATTCAACACACCGCAAAAATCAGGCGGTTTTTTTTTGCTTGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCaa                                     <  1:5408893/149‑1 (MQ=37)
                                                                                                                  cACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATTCAAGCGCGTCTCAGCGGAGACAAAACACGCCAGTAACTCTCTTTTTGTCaa                                     >  1:11442/1‑150 (MQ=255)
                                                                                                                  cACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCaa                                     <  2:2948084/150‑1 (MQ=255)
                                                                                                                  cACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCaa                                     <  2:1287660/150‑1 (MQ=255)
                                                                                                                        ccTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCaaaa                               <  2:178802/150‑1 (MQ=255)
                                                                                                                          tttAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAaga                             >  2:3600052/1‑150 (MQ=255)
                                                                                                                           ttAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAagag                            >  2:5012508/1‑150 (MQ=255)
                                                                                                                            tAAAAAACAAACCACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGGAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAagaga                           <  2:3319923/150‑1 (MQ=255)
                                                                                                                            tAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAagaga                           <  2:5759997/150‑1 (MQ=255)
                                                                                                                              aaaaacaaaaaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGt                         >  2:6648837/1‑150 (MQ=255)
                                                                                                                               aaaacaaaaaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTcc                                                                                                                 <  2:6309072/61‑1 (MQ=255)
                                                                                                                                aaacaaaaaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGtta                       <  2:2401246/150‑1 (MQ=255)
                                                                                                                                  acaaaaaacCGCCTGATTTTGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGttatt                     >  1:6149036/1‑150 (MQ=37)
                                                                                                                                    aaaacacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTtattat                   <  2:5896846/150‑1 (MQ=255)
                                                                                                                                    aaaaaacCGCCTGATTTTGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGtct                                                                             >  1:2828309/1‑92 (MQ=21)
                                                                                                                                    aaaaaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCATTCCAGCCGGTTTTATTACCGTAGTAATTCAAGCGCGTCTCAGCGGAGACAATACTAGCCAGTAACTATCTTTTTTTCAAGACAAAGAGAGGtattat                   >  1:2304334/1‑150 (MQ=255)
                                                                                                                                     aaaaacCGCCTGATTTTGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTattatt                  >  2:2754267/1‑150 (MQ=25)
                                                                                                                                     aaaaacCGCCAAATTCGGCGGTTTTTTATTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTa                                                                                                   <  1:3163214/68‑1 (MQ=14)
                                                                                                                                      aaaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTg                 >  2:798947/1‑150 (MQ=255)
                                                                                                                                       aaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGt                <  1:2954219/150‑1 (MQ=255)
                                                                                                                                       aaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGt                >  1:4497168/1‑150 (MQ=255)
                                                                                                                                           cGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTg            >  2:3065294/1‑150 (MQ=255)
                                                                                                                                           cGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTAGCGGCCTTTACAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCACAGCGGAGACAATACTCGCCAGTAAATCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTATg            >  1:1830053/1‑150 (MQ=37)
                                                                                                                                            gCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGt           >  2:4045561/1‑150 (MQ=255)
                                                                                                                                            gCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGt           >  2:2823411/1‑150 (MQ=255)
                                                                                                                                             ccAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGtt          <  2:5998865/150‑1 (MQ=255)
                                                                                                                                                aaaTCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAGTGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGTTAgt       <  2:4850988/150‑1 (MQ=255)
                                                                                                                                                aaaTCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGTTAgt       <  2:6538546/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGTTAGTGTATT  >  NC_000913/2653549‑2653847

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