Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,653,697 C→T 5.3% intergenic (+358/‑158) 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,6970CT 5.3% 584.6 / 12.0 186intergenic (+358/‑158)sseA/ryfA3‑mercaptopyruvate sulfurtransferase/novel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base C (84/92);  new base T (5/5);  total (89/97)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

Only 100 of 189 total aligned reads displayed.

GTGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGTTAGTGTAT  >  NC_000913/2653548‑2653846
                                                                                                                                                     |                                                                                                                                                     
gtgtTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATc                                                                                                                                                       >  1:336678/1‑150 (MQ=255)
gtgtTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATc                                                                                                                                                       <  2:630481/150‑1 (MQ=255)
 tgtTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCa                                                                                                                                                      <  2:223301/150‑1 (MQ=255)
  gtTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAg                                                                                                                                                     >  2:1948050/1‑150 (MQ=255)
  gtTCGGTAATAAAATAACCTGGCTTATTAGTCAGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAg                                                                                                                                                     <  1:1140294/150‑1 (MQ=255)
   ttCTGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAgg                                                                                                                                                    <  1:1998595/150‑1 (MQ=255)
     cGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCg                                                                                                                                                  <  2:527987/150‑1 (MQ=255)
      ggTAATAAAATAACCTTGCTTATTAGTCCGAATTCATACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCgg                                                                                                                                                 <  1:51167/150‑1 (MQ=255)
      ggTAATAAAATAACCTGTCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCgg                                                                                                                                                 <  2:1062093/150‑1 (MQ=255)
         aataaaataaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCTGATTTTGGCGGttt                                                                                                                                              <  2:1383845/150‑1 (MQ=255)
         aataaaataaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGttt                                                                                                                                              <  2:812897/150‑1 (MQ=255)
           taatataaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGttttt                                                                                                                                            <  2:1987915/150‑1 (MQ=255)
           taaaataaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTTCTCAAAATTAAAAATTCTAACCGGTAACAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAATCAAAAAACCGCCTGATTTTGGCGGttttt                                                                                                                                            <  1:292376/150‑1 (MQ=255)
             aaataaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTttg                                                                                                                                          <  1:1698529/150‑1 (MQ=255)
             aaataaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTttg                                                                                                                                          >  2:83770/1‑150 (MQ=255)
              aataaCCTGGATTATTAGTCCGAATTCAGACAAATATAAATAAATACTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTttgt                                                                                                                                         <  2:1997756/150‑1 (MQ=255)
               ataaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGTTATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTttgtt                                                                                                                                        <  2:717008/150‑1 (MQ=255)
                  aCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCt                                                                                                                                     <  2:542341/150‑1 (MQ=255)
                    cTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTgg                                                                                                                                   <  1:1718327/150‑1 (MQ=255)
                    cTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGACAAAATCAGGCGGTTTATTGTTGCTgg                                                                                                                                   >  2:595438/1‑150 (MQ=255)
                    cTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATAAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTgg                                                                                                                                   <  2:2015801/150‑1 (MQ=255)
                      ggCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTc                                                                                                                                 >  1:2003372/1‑150 (MQ=255)
                      ggCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAAACGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTc                                                                                                                                 <  2:1498796/150‑1 (MQ=255)
                      ggCTTATTAGTCCGAATTCAGACAAAAATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAAGATGTAACTTCACTTTCCTTTAAAAAACAAAACACCGCCTCATTTTGGCGGTTTTttgttttttaaa                                                                                                                                 <  1:224127/150‑8 (MQ=38)
                       gCTTATTATTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTcc                                                                                                                                <  1:443523/150‑1 (MQ=255)
                       gCTTATTAGTCCTAATTCAGACAAATATAAATAAATCATGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTCTTTGTTGCTGGTcc                                                                                                                                <  1:190287/150‑1 (MQ=255)
                       gCTTATTAGTCAGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTcc                                                                                                                                <  2:1072747/150‑1 (MQ=255)
                        cTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCg                                                                                                                               <  1:1651716/150‑1 (MQ=255)
                        cTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCg                                                                                                                               >  1:655870/1‑150 (MQ=255)
                        cTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCg                                                                                                                               >  2:791164/1‑150 (MQ=255)
                        cTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCgg                                                                                                                              <  2:1721842/150‑1 (MQ=255)
                              aGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTcgc                                                                                                                         >  2:169918/1‑150 (MQ=255)
                              aGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTcgc                                                                                                                         >  2:1541820/1‑150 (MQ=255)
                               gTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTcgcg                                                                                                                        >  2:1897782/1‑150 (MQ=255)
                               gTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTcgcg                                                                                                                        <  1:1593472/150‑1 (MQ=255)
                                tCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCgg                                                                                                                       <  2:939967/150‑1 (MQ=255)
                                 ccGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATTTAACTTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCTGTTTTTTGTTTCTGGTCCGGTTCGCGGc                                                                                                                      >  2:119851/1‑150 (MQ=255)
                                 ccGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCTCCAAAATCAGGCGGTTTTTTGTTGCTGGTCAGGTTCGCGGc                                                                                                                      >  1:243454/1‑150 (MQ=255)
                                  cGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGcc                                                                                                                     >  1:240880/1‑150 (MQ=255)
                                    aaTTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGACAAAATCAGGCGGTTTTTTGTTGCTCGTCCGGTTCGCGGACtt                                                                                                                   >  2:2203979/1‑150 (MQ=255)
                                      ttCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTc                                                                                                                 >  1:2113416/1‑150 (MQ=255)
                                        cAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTccc                                                                                                               <  1:1625367/150‑2 (MQ=255)
                                          gACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCcagc                                                                                                             >  2:1524228/1‑150 (MQ=255)
                                          gACAAATATAAATAAATCCTGCTCAAAATTAAAAATTATAAGCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGTGTTTTTGTTGCTGGTCCGGTTCGACGCCTTTCcagc                                                                                                             >  2:130231/1‑150 (MQ=255)
                                                    aataaatCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGAGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTAtt                                                                                                   >  1:2138105/1‑150 (MQ=255)
                                                            cTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCgtagt                                                                                           >  1:1469420/1‑150 (MQ=255)
                                                               cTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTATTAAt                                                                                        >  2:2084917/1‑150 (MQ=255)
                                                                   aaaTTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTATGTTGCTGGTCCGGTTCGCTGCCTTTCCAGCAGGTTTTATTACCGTAGTAATGCaa                                                                                    >  2:400406/1‑150 (MQ=255)
                                                                     aTTAAAAATTCTTACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAgc                                                                                  <  2:1605480/150‑1 (MQ=255)
                                                                     aTTAAAAATTCTAACCGGTAAAAGATATTACTTAAACTTGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAgc                                                                                  <  2:1504738/150‑1 (MQ=255)
                                                                       tAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAgcgc                                                                                >  1:796622/1‑150 (MQ=255)
                                                                       tAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAgcgc                                                                                >  2:635781/1‑150 (MQ=255)
                                                                        aaaaaTTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGcgcg                                                                               <  2:375102/150‑1 (MQ=255)
                                                                            aTTCTAACCGGTAAATGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCTCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGtctc                                                                           <  2:346421/150‑1 (MQ=255)
                                                                                tAAACGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCg                                                                       <  2:32500/150‑1 (MQ=255)
                                                                                 aaCCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTAAGCgg                                                                      >  2:1358053/1‑150 (MQ=255)
                                                                                   ccGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGgag                                                                    <  2:1751717/150‑1 (MQ=255)
                                                                                    cGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGgaga                                                                   <  1:1963167/150‑1 (MQ=255)
                                                                                     ggTAAAAGTTATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGAc                                                                  <  1:576750/150‑1 (MQ=255)
                                                                                     ggTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGAc                                                                  >  2:11289/1‑150 (MQ=255)
                                                                                     ggTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGAc                                                                  <  1:1897782/150‑1 (MQ=255)
                                                                                        aaaaGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAAt                                                               >  1:1456966/1‑150 (MQ=255)
                                                                                         aaaGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATa                                                              >  1:1202/1‑150 (MQ=255)
                                                                                          aaGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTATCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATAc                                                             >  1:522031/1‑150 (MQ=255)
                                                                                             aaattaattaAACTTTTAAATTCAATTTCCTTTAAAATACAAAAAAACGCCATAATTTGGCGGTTTTTTTTTGCTTGTCTGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCg                                                          <  2:683142/148‑1 (MQ=25)
                                                                                               atTACTTAAACATTTATATTCACTTTCCTTTAAAAAACAAACAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGcc                                                        <  1:1679025/150‑1 (MQ=255)
                                                                                               atTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGcc                                                        <  1:83770/150‑1 (MQ=255)
                                                                                                  aCTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGt                                                     >  2:1159807/1‑150 (MQ=255)
                                                                                                   cTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTa                                                    >  2:559933/1‑150 (MQ=255)
                                                                                                    ttAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTaa                                                   >  1:1080623/1‑150 (MQ=255)
                                                                                                     tAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCTCCAGTaaa                                                  >  1:1116819/1‑149 (MQ=255)
                                                                                                        aCATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAActct                                               >  2:888812/1‑150 (MQ=255)
                                                                                                           tGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCtt                                            >  2:1840491/1‑150 (MQ=255)
                                                                                                            gTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGAGGTGTTTTGTTGCTGGTCAGTTTCGAGGCCTTTCAAGAAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGAGGATACAATACTCGCCAGAAACTCTCttt                                           >  1:406543/1‑150 (MQ=255)
                                                                                                              aaaTTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTACCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCttttt                                         >  2:1451088/1‑150 (MQ=255)
                                                                                                                       ttCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCaa                                <  1:2203979/150‑1 (MQ=255)
                                                                                                                        tCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCaaa                               <  1:1643459/150‑1 (MQ=255)
                                                                                                                         ccTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCaaaa                              >  2:1934366/1‑150 (MQ=255)
                                                                                                                         ccTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTAAAGCaaaa                              >  2:1793525/1‑150 (MQ=255)
                                                                                                                          cTTTAAAAAACAAAATACCGCCTACATCAGGCTGTTTTTTGTTGCTTGTCCTGTTCGCGGCCTTTCCAGCATGTTGTATTACCGTAGTAATGCAAGGGCGTCTCATCGTATACAATACTCTCCAGTAACTCTCTTTTTGTCAAGCAAAag                             <  2:1890307/150‑1 (MQ=37)
                                                                                                                          cTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAag                             >  2:613216/1‑150 (MQ=255)
                                                                                                                          cTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAag                             >  2:1682117/1‑150 (MQ=255)
                                                                                                                              aaaaaacaaaaaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCACCAGGTTGTATTACAGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAagagat                         >  2:569196/1‑149 (MQ=255)
                                                                                                                                aaatcaaaaaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGtt                       <  2:252871/150‑1 (MQ=255)
                                                                                                                                aaaacaaaaaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGtt                       >  1:945966/1‑150 (MQ=255)
                                                                                                                                 ccaaaaaaaacCGCCATAATCTGGCTGTTTTTTTTTGCTTGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTATTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGtta                      <  1:777155/146‑1 (MQ=37)
                                                                                                                                      aaaaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTattatt                 >  1:1153299/1‑150 (MQ=255)
                                                                                                                                      aaaaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTattatt                 >  1:2101101/1‑150 (MQ=255)
                                                                                                                                       aaaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTg                <  1:1541820/150‑1 (MQ=255)
                                                                                                                                        aaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGt               <  1:989526/150‑1 (MQ=255)
                                                                                                                                          acCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGTAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTc             <  2:2178980/150‑1 (MQ=255)
                                                                                                                                          acCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTc             <  1:567251/150‑1 (MQ=255)
                                                                                                                                             gCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGt          >  1:2042700/1‑150 (MQ=255)
                                                                                                                                              ccttcaTCAGGGGGTTTTTTGTTTCTTGTTCGGTTCGCGGCCTTTCCTGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGtt         <  1:1764846/145‑1 (MQ=38)
                                                                                                                                              ccAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGtt         <  1:1229298/150‑1 (MQ=255)
                                                                                                                                              ccAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGtt         >  2:1629624/1‑150 (MQ=255)
                                                                                                                                              ccAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGtt         <  1:717790/150‑1 (MQ=255)
                                                                                                                                               cAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGTTa        <  2:662564/150‑1 (MQ=255)
                                                                                                                                               cAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGTTa        >  1:1402623/1‑150 (MQ=255)
                                                                                                                                                     cAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGAAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGTTAGTGtat  >  2:1114140/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTCTGTTAGTGTAT  >  NC_000913/2653548‑2653846

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: