Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,653,695 A→T 6.1% intergenic (+356/‑160) 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,6950AT 6.1% 1178.1 / 41.3 367intergenic (+356/‑160)sseA/ryfA3‑mercaptopyruvate sulfurtransferase/novel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base A (181/160);  new base T (11/13);  total (192/175)
Fisher's exact test for biased strand distribution p-value = 5.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 373 total aligned reads displayed.

TGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGTTAGT  >  NC_000913/2653549‑2653842
                                                                                                                                                  |                                                                                                                                                   
  ttCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAATCAAAAAACCGCCAAAATCAgg                                                                                                                                                >  1:4742061/1‑150 (MQ=255)
  ttCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAgg                                                                                                                                                >  1:3381089/1‑150 (MQ=255)
  ttCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAgg                                                                                                                                                >  1:4198504/1‑150 (MQ=255)
  ttCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCACAAAATTAAAAATAATAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAATAAAACGCCAAAATCAgc                                                                                                                                                >  1:4926688/1‑149 (MQ=255)
    cGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCg                                                                                                                                              <  2:4270952/150‑1 (MQ=255)
    cGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCg                                                                                                                                              >  1:4831891/1‑150 (MQ=255)
       taataaAATAACCTGGCTTATTAGACCGAATTCAGACAAATATAAATAAATCCTGCTCAAAAATAAAAATTCTAACCGGTAAACGCTATTACTTAAACATGGAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGtt                                                                                                                                           <  2:4454533/150‑1 (MQ=255)
         ataaaataaCCTGTCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATTTAAATTCACTTTCCTTTAAAAAACAAAAAACAGCCAAAATCAGGCGGtttt                                                                                                                                         >  1:3737118/1‑150 (MQ=255)
         ataaaataaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGtttt                                                                                                                                         <  1:5494600/150‑1 (MQ=255)
          taaaataaCCTGGATTATTAGTCCGAATTCAGACAAATATAAAGAACTCCTGCTCAAATTTAATAATTCTATCCGGTAATATATATTACTTAATCATGTAAATTCTCTTTCCTTTATAAATCATACTACCGCCTGATTTTGGCGGttttt                                                                                                                                        <  2:4116265/150‑1 (MQ=255)
           aaaataaCCTGGATTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGtttttt                                                                                                                                       <  2:6070984/150‑1 (MQ=255)
             aataaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTGTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGttgtttgt                                                                                                                                     >  1:6697132/1‑150 (MQ=255)
              ataaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTttgtt                                                                                                                                    <  2:2033783/150‑1 (MQ=255)
              ataaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTttgtt                                                                                                                                    <  1:4042249/150‑1 (MQ=255)
               taaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTttgttg                                                                                                                                   >  1:4586290/1‑150 (MQ=255)
                caCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGc                                                                                                                                  <  2:2168256/149‑1 (MQ=255)
                 aCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGCAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCt                                                                                                                                 <  2:2494879/150‑1 (MQ=255)
                    tGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATACTGCTCAAAATTAAAAAATCTAACCGGTAAAATATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAAACGCCAAAATCAGGCGGTTTTTTGTTGATGGt                                                                                                                              >  1:5998865/1‑150 (MQ=255)
                     ggCTTATTAGTCCTAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTc                                                                                                                             <  2:6470734/150‑1 (MQ=255)
                      gCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGc                                                                                                                                  <  2:3556054/144‑1 (MQ=255)
                      gCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGc                                                                                                                                  >  1:3556054/1‑144 (MQ=255)
                        ttattaGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCCAAATCAGGCGGTTTTTTGTTGCTGGTCCgg                                                                                                                          >  2:3915046/1‑150 (MQ=255)
                        ttattaGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAAAATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCgg                                                                                                                          >  1:3163855/1‑150 (MQ=255)
                          attaGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAACACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGtt                                                                                                                        <  1:194669/150‑1 (MQ=255)
                               tCCGAATTCAGACAGATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCgg                                                                                                                   >  2:3124583/1‑150 (MQ=255)
                                ccGAACTCAGACAAAGATAAATAAATCCTGCACAAAATTAAAAATTCTAACCGGTAACAGATATTACTTAAACACCGAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGATTTTTGTTGCTGGTCCGGTTCGCGGc                                                                                                                  <  1:6393660/150‑1 (MQ=255)
                                      tCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTcc                                                                                                            <  1:1285777/150‑1 (MQ=255)
                                        aGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCcag                                                                                                          <  2:5091368/150‑1 (MQ=255)
                                        aGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCcag                                                                                                          >  1:2008547/1‑150 (MQ=255)
                                         gACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCcagc                                                                                                         >  2:4073917/1‑150 (MQ=255)
                                          aCAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCcagca                                                                                                        >  2:4824673/1‑150 (MQ=255)
                                            aaaTATCAATAAAACCAGCTCAAAATTAAAAATTCTAACAGGTAAAAGATATTACTGCAACATCTAAATTCACATTCCATTAAAAAACAAAAAAACGCCTAATTCAGGCGGTATGTAGTTTCTGGTCCGGTTCGCGTACTTTCcagcatg                                                                                                      >  2:159980/1‑148 (MQ=255)
                                              atatAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTGCGCGGCCTTTCCAGCAGGtt                                                                                                    >  2:4486563/1‑150 (MQ=255)
                                                ataaataaATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGt                                                                                                  <  2:2117897/150‑1 (MQ=255)
                                                ataaataaATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGt                                                                                                  >  1:5182633/1‑150 (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)
                                                          ccTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCgtag                                                                                        >  1:1618956/1‑150 (MQ=255)
                                                            tGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCgtagta                                                                                      >  2:5525977/1‑150 (MQ=255)
                                                              cTCAAAATTAAAAATTCTAACCGGGAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAAt                                                                                    <  2:4694583/150‑1 (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)
                                                                        aaaaTTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGt                                                                          >  1:3060066/1‑150 (MQ=255)
                                                                        aaaaTTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGt                                                                          >  2:3653970/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)
                                                                               tAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCg                                                                   >  2:1899376/1‑150 (MQ=255)
                                                                               tAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCg                                                                   >  2:1995493/1‑150 (MQ=255)
                                                                                 cccGGGAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGga                                                                 <  2:3029716/149‑1 (MQ=255)
                                                                                 aCCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGga                                                                 >  1:4168086/1‑150 (MQ=255)
                                                                                   cGGTCACAGATCTTACTGACACATGTAACTTCACTTTCCTTTAAACAACACACCACCGCCAACCTCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGgaga                                                               <  2:6697132/150‑1 (MQ=37)
                                                                                    ggTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGAc                                                              >  1:5649113/1‑150 (MQ=255)
                                                                                        aaaGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATa                                                          <  1:3486896/150‑1 (MQ=255)
                                                                                        aaaGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATa                                                          >  1:2105044/1‑150 (MQ=255)
                                                                                                  cTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTa                                                >  1:1968049/1‑150 (MQ=255)
                                                                                                  attAAACATGTAAATTCACTTTCCTTTAAAAAACAAAACACCGCCAAAATCTGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTa                                                <  1:6613293/149‑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:1287660/150‑1 (MQ=255)
                                                                                                                  cACTTTCATTTAAAAAACAAAAAACCGCCAAAATCAGGCGTTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCaa                                <  2:294059/150‑1 (MQ=255)
                                                                                                                        ccTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCaaaa                          <  2:178802/150‑1 (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)
                                                                                                                                  acaaaaaacCGCCTGATTTTGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGttatt                >  1:2467574/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)
                                                                                                                                      aaaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTg            >  2:798947/1‑150 (MQ=255)
                                                                                                                                       aaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGt           <  1:2954219/150‑1 (MQ=255)
                                                                                                                                       aaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGt           >  1:4497168/1‑150 (MQ=255)
                                                                                                                                        aacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGtt          <  1:2768322/150‑1 (MQ=255)
                                                                                                                                          ccGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCt        <  1:5508474/150‑1 (MQ=255)
                                                                                                                                           cGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTAGCGGCCTTTACAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCACAGCGGAGACAATACTCGCCAGTAAATCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTATg       >  1:1830053/1‑150 (MQ=37)
                                                                                                                                            gCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGt      >  2:2823411/1‑150 (MQ=255)
                                                                                                                                             ccAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGtt     >  2:397754/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)
                                                                                                                                                aaaTCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGTTAgt  <  1:1497707/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGTTAGT  >  NC_000913/2653549‑2653842

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: