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

Only 100 of 369 total aligned reads displayed.

CGCCGTGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGTTA  >  NC_000913/2653544‑2653840
                                                                                                                                                    |                                                                                                                                                    
cgccgTGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCaa                                                                                                                                                     >  1:3310757/1‑150 (MQ=255)
       ttCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAATCAAAAAACCGCCAAAATCAgg                                                                                                                                              >  1:4742061/1‑150 (MQ=255)
       ttCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAgg                                                                                                                                              >  1:4198504/1‑150 (MQ=255)
       ttCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCACAAAATTAAAAATAATAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAATAAAACGCCAAAATCAgc                                                                                                                                              >  1:4926688/1‑149 (MQ=255)
         cGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCg                                                                                                                                            >  1:4831891/1‑150 (MQ=255)
         atttaataaAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGCTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCg                                                                                                                                            <  2:5702037/147‑1 (MQ=255)
             aataaaataaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGttt                                                                                                                                        >  2:2786928/1‑150 (MQ=255)
              ataaaataaCCTGTCTTATTAGTCCGAATTCAGACACATATACATAAATCCTGCTCAAAATTCAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCTGATTTTGGCGGtttt                                                                                                                                       <  2:3928817/150‑1 (MQ=255)
              ataaaataaCCTGTCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATTTAAATTCACTTTCCTTTAAAAAACAAAAAACAGCCAAAATCAGGCGGtttt                                                                                                                                       >  1:3737118/1‑150 (MQ=255)
               taaaataaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCTGATTTTGGCGGttttt                                                                                                                                      <  2:4116103/150‑1 (MQ=255)
               taaaataaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCTGATTTTGGCGGttttt                                                                                                                                      <  2:4718591/150‑1 (MQ=255)
               taaaataaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGttttt                                                                                                                                      <  2:490683/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)
                     aaCCTGGCTTATTAGTCCTAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGc                                                                                                                                <  1:5504943/150‑1 (MQ=255)
                      aCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGCAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCt                                                                                                                               <  2:2494879/150‑1 (MQ=255)
                       ccTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTg                                                                                                                              <  2:5097813/150‑1 (MQ=255)
                           gCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGc                                                                                                                                >  1:3556054/1‑144 (MQ=255)
                           gCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTcc                                                                                                                          <  1:1982705/150‑1 (MQ=255)
                           gCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCGACCCGGTACAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTcc                                                                                                                          <  1:5780663/150‑1 (MQ=255)
                                ttaGGCCGAATTCAGACAAATATAAATACCTCCTGCTCAAAATTACAAAGTCTAACCGGTAAACGAGAGTACTGAAACCTGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGGTc                                                                                                                     <  1:27730/150‑1 (MQ=255)
                                  aGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTcgc                                                                                                                   <  2:3504072/150‑1 (MQ=255)
                                    tCCGAATTCAGACAGATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCgg                                                                                                                 >  2:3116122/1‑150 (MQ=255)
                                    tCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCgg                                                                                                                 >  2:3868796/1‑150 (MQ=255)
                                       gAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCt                                                                                                              >  1:5180898/1‑150 (MQ=255)
                                          ttCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTc                                                                                                           <  1:6323195/150‑1 (MQ=255)
                                           tCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTcc                                                                                                          <  1:3010692/150‑1 (MQ=255)
                                           tCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTcc                                                                                                          >  2:5594484/1‑150 (MQ=255)
                                           tCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTcc                                                                                                          <  1:1285777/150‑1 (MQ=255)
                                                 aaaTATCAATAAAACCAGCTCAAAATTAAAAATTCTAACAGGTAAAAGATATTACTGCAACATCTAAATTCACATTCCATTAAAAAACAAAAAAACGCCTAATTCAGGCGGTATGTAGTTTCTGGTCCGGTTCGCGTACTTTCcagcatg                                                                                                    >  2:159980/1‑148 (MQ=255)
                                                 aaaTATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCAGTTTGCTTTAAAATACATAAAACCGCCAAAATCTGGCGGTTTTTTGTTGCTGGTCGGGTTCGCGGCCTTTCCAGCAgg                                                                                                    >  1:3116540/1‑150 (MQ=39)
                                                   atatAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGtt                                                                                                  >  2:6038486/1‑150 (MQ=255)
                                                     ataaataaATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGt                                                                                                <  2:2117897/150‑1 (MQ=255)
                                                      taaataaaTCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTa                                                                                               <  2:900738/150‑1 (MQ=255)
                                                         ataaatCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTa                                                                                            >  2:2327892/1‑150 (MQ=255)
                                                         ataaatCCGGCTCAAAATTAAAAAGTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTa                                                                                            <  1:2364880/150‑1 (MQ=255)
                                                          taaatCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTAc                                                                                           >  2:1524488/1‑150 (MQ=255)
                                                              tCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACGTGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCgta                                                                                       <  2:4937105/150‑1 (MQ=255)
                                                              tCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCgta                                                                                       <  2:5996265/150‑1 (MQ=255)
                                                               ccTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCgtag                                                                                      >  1:5487564/1‑150 (MQ=255)
                                                               ccTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCgtag                                                                                      >  2:4798353/1‑150 (MQ=255)
                                                                cTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCgtagt                                                                                     >  2:4825215/1‑150 (MQ=255)
                                                                   cTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAAt                                                                                  <  2:3238822/150‑1 (MQ=255)
                                                                    tCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTcc                                                                                                                          >  1:863973/1‑109 (MQ=255)
                                                                    tCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATg                                                                                 >  1:5127145/1‑150 (MQ=255)
                                                                      aaaaTTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCa                                                                               >  1:6201143/1‑150 (MQ=255)
                                                                           tAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAgcgc                                                                          <  2:1076814/150‑1 (MQ=255)
                                                                           tAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAgcgc                                                                          <  1:1650840/150‑1 (MQ=255)
                                                                           tAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAgcgc                                                                          <  1:3439541/150‑1 (MQ=255)
                                                                             aaaaTTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGt                                                                        <  1:2549243/150‑1 (MQ=255)
                                                                             aaaaTTCTAACCCCAAAAACATATTACTTAAACATCTAAATACACTTTCCTTTAAAAAACAAAAAACCCCCAAAATCACGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGt                                                                        <  2:620291/150‑1 (MQ=38)
                                                                               aaTTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGtct                                                                      >  2:6465581/1‑150 (MQ=255)
                                                                                 ttCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCa                                                                    >  1:5350399/1‑150 (MQ=255)
                                                                                   cTAACCGGTAACAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGTGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGc                                                                  <  2:5024633/150‑1 (MQ=255)
                                                                                   cTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGc                                                                  <  1:1041445/150‑1 (MQ=255)
                                                                                    tAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCg                                                                 >  2:1899376/1‑150 (MQ=255)
                                                                                    tAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCg                                                                 >  2:1994913/1‑150 (MQ=255)
                                                                                      cccGGGAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGga                                                               <  2:3029716/149‑1 (MQ=255)
                                                                                      aCCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGga                                                               <  2:4260872/150‑1 (MQ=255)
                                                                                      aCCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGga                                                               >  1:4168086/1‑150 (MQ=255)
                                                                                      aCCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGga                                                               <  1:1524488/150‑1 (MQ=255)
                                                                                      aCCGGTAAAAGATATTACTTAAACATGTAAATACACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGga                                                               <  2:5439919/150‑1 (MQ=255)
                                                                                        cgcgaacaGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAACACCGCCAAATTCTGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGgaga                                                             <  1:5979647/143‑1 (MQ=38)
                                                                                        cGGTCACAGATCTTACTGACACATGTAACTTCACTTTCCTTTAAACAACACACCACCGCCAACCTCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGgaga                                                             <  2:6697132/150‑1 (MQ=37)
                                                                                         ggTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGAc                                                            >  1:5649113/1‑150 (MQ=255)
                                                                                             aaaGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATa                                                        <  1:3486896/150‑1 (MQ=255)
                                                                                               aGATATTACTTATTCTTGTAAATTCACTTTCATTTATAAAACAACATACCGGCAAAATCAGGCGGTTTTTTGTTTCTGGTCGGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACt                                                      <  2:4378914/150‑1 (MQ=255)
                                                                                               aGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTCGTAATGCAAGCGCGTCTCAGCGGAGACAATAcg                                                      >  2:5832814/1‑149 (MQ=255)
                                                                                                  tatTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGc                                                   <  1:2870729/150‑1 (MQ=255)
                                                                                                       cTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTa                                              <  2:4479056/150‑1 (MQ=255)
                                                                                                              aTGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCAAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTAAGCGGAGAAAATACTCGCCAGTAActctct                                       >  1:2141128/1‑150 (MQ=255)
                                                                                                               tGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAAt                                                         <  2:6482572/131‑1 (MQ=255)
                                                                                                               tGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAAt                                                         >  1:6482572/1‑131 (MQ=255)
                                                                                                               tGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCtt                                      >  1:1431651/1‑150 (MQ=255)
                                                                                                                    aTTCACTTTCCTTTAAAATTCAATATACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTGGGGTTCGCGGGCTTTGGAGGAGGTTGTATTTGGTTTGTTATGGTAGCGCGTCTCTGCGGAGAGAATAGTCGCGAGTAAGTCTCTTTTTGt                                 >  2:1385052/1‑150 (MQ=255)
                                                                                                                     ttCACTTTCCTTTAAAAAACAATAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTGCGCGGCCTTTCCAGCAGGTTGTATTACCGGAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTc                                >  1:3818729/1‑150 (MQ=255)
                                                                                                                       gaCTTTCCTTAATAATTCAACACACCGCAAAAATCAGGCGGTTTTTTTTTGCTTGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCaa                              <  1:5408893/149‑1 (MQ=37)
                                                                                                                       cACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATTCAAGCGCGTCTCAGCGGAGACAAAACACGCCAGTAACTCTCTTTTTGTCaa                              >  1:11442/1‑150 (MQ=255)
                                                                                                                       cACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCaa                              <  1:1869153/150‑1 (MQ=255)
                                                                                                                       cACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCaa                              <  2:1287660/150‑1 (MQ=255)
                                                                                                                       cACTTTCATTTAAAAAACAAAAAACCGCCAAAATCAGGCGTTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCaa                              <  2:294059/150‑1 (MQ=255)
                                                                                                                        aCTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAg                             >  1:1768057/1‑150 (MQ=255)
                                                                                                                         cTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGATGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGc                            >  2:3466133/1‑150 (MQ=255)
                                                                                                                            tCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCaaa                         <  1:6654886/150‑1 (MQ=255)
                                                                                                                             ccTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCaaaa                        <  2:178802/150‑1 (MQ=255)
                                                                                                                               tttAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTTAAGCAAAaga                      >  1:1958215/1‑150 (MQ=255)
                                                                                                                               tttAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAaga                      >  2:3676246/1‑150 (MQ=255)
                                                                                                                               tttAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAaga                      >  2:3600052/1‑150 (MQ=255)
                                                                                                                                ttAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAagag                     >  2:4610583/1‑150 (MQ=255)
                                                                                                                                 tAAAAAACAAACCACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGGAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAagaga                    <  2:3319923/150‑1 (MQ=255)
                                                                                                                                   aaaaacaaaaaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGt                  >  2:6648837/1‑150 (MQ=255)
                                                                                                                                       acaaaaaacCGCCTGATTTTGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGttatt              >  1:6149036/1‑150 (MQ=37)
                                                                                                                                          aaaaacCGCCTGATTTTGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTattatt           >  2:2754267/1‑150 (MQ=25)
                                                                                                                                          aaaaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCa                                                 >  1:201260/1‑112 (MQ=255)
                                                                                                                                           aaaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTg          >  2:798947/1‑150 (MQ=255)
                                                                                                                                            aaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGt         >  1:4497168/1‑150 (MQ=255)
                                                                                                                                                 gCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGt    >  2:4045561/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGCCGTGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGTTA  >  NC_000913/2653544‑2653840

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: