Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,628,6490AC25.7% 162.6 / 13.7 74intergenic (+102/+206)nadR/ettADNA‑binding transcriptional repressor/NMN adenylyltransferase NadR/energy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base A (47/8);  new base C (0/19);  total (47/27)
Fisher's exact test for biased strand distribution p-value = 1.04e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACA  >  NC_000913/4628511‑4628788
                                                                                                                                          |                                                                                                                                              
ggAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAAc                                                                                                                                              <  2:369834/140‑1 (MQ=255)
       gTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGGAAAAACCTGCACg                                                                                                                                       >  1:1070141/1‑140 (MQ=255)
       gTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACg                                                                                                                                       >  2:1567197/1‑140 (MQ=255)
       gTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACg                                                                                                                                       >  1:231007/1‑140 (MQ=255)
       gTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACg                                                                                                                                       >  1:1322740/1‑140 (MQ=255)
            ggAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGa                                                                                                                                >  1:1108199/1‑140 (MQ=255)
              agaTGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAACCCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGAcc                                                                                                                              <  1:517566/140‑1 (MQ=255)
              agaTGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGAcc                                                                                                                              >  1:1635062/1‑140 (MQ=255)
              agaTGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGAcc                                                                                                                              >  1:1774171/1‑140 (MQ=255)
                 tgatgGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGc                                                                                                                           >  1:475702/1‑140 (MQ=255)
                  gatgGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGcc                                                                                                                          >  1:1680/1‑140 (MQ=255)
                   atgGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAACCCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAACCCTGCACGGGTGACCCGCCa                                                                                                                         <  1:1499520/140‑1 (MQ=255)
                        gaaGCAGAGATAACCGTGATAAAACTCCTAAAAGGCGAGGTAAAAAATGAGTTTTTTTGATGAGTTAAAACCCTCTCGGGAAGAGGCTGTCGAGATTAAACAGGGTTTAAAAAACCCTGCAGGGGTGCCCCGCCACGaaa                                                                                                                    <  2:1391076/138‑1 (MQ=255)
                           gCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTg                                                                                                                 >  2:136521/1‑140 (MQ=255)
                            cAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGa                                                                                                                >  2:1395249/1‑140 (MQ=255)
                            cAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGa                                                                                                                >  2:838864/1‑140 (MQ=255)
                             agagATACCCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAACCCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCAGGGGTGACCCGCCACGAAATTGAg                                                                                                               <  1:523372/140‑1 (MQ=255)
                             agagATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAACCCTGCACGGGTGACCCGCCACGAAATTGAg                                                                                                               <  1:88211/140‑1 (MQ=255)
                                 aTAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATg                                                                                                           >  2:357051/1‑140 (MQ=255)
                                   accggTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAACCCTCTCTGGAAGAGGCTGTCGAGATAAAACAAGGTTTAAAAAACCCTGCACGGGTGCCCCGCCACGAAATTGAGGATGCt                                                                                                         <  1:1356977/136‑1 (MQ=255)
                                   aaCCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCt                                                                                                         <  1:1567197/140‑1 (MQ=255)
                                   aaCCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTAAAAAAACCTGCAGGGGTGACCCGCCACGAAATTGAGGATGCt                                                                                                         <  1:310999/140‑1 (MQ=255)
                                     ccGTGATGAAACTACTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTaa                                                                                                       >  2:725813/1‑140 (MQ=255)
                                         gatgaAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCt                                                                                                   >  1:266462/1‑140 (MQ=255)
                                          atgaAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCtg                                                                                                  >  1:520553/1‑140 (MQ=255)
                                           tgaAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCtgt                                                                                                 >  2:1074102/1‑140 (MQ=255)
                                           tgaAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCtgt                                                                                                 >  1:612622/1‑140 (MQ=255)
                                               aCTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTa                                                                                             >  1:229783/1‑140 (MQ=255)
                                                  gCTCAAGGGGGAGGTAAAAAATGATTTTTTTTGATGAGTTAAAACCCTCTCTGGAAGAGGCTGTGGAGATTAAACAGGGTTTAAAAAACCCTGCGCGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGAc                                                                                          <  2:877274/140‑1 (MQ=255)
                                                   cTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGAcc                                                                                         >  2:215101/1‑140 (MQ=255)
                                                    tCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAACCCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAACCCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCg                                                                                        <  1:1504806/140‑1 (MQ=255)
                                                    aaaaaGGGGAGGTAAAAAATGATTTTTTTTGAGGATTTAAAACCCTCTCTGAAGGAGGCTGTGGAAATAAAACAAGTTTTAAAAAAACCTGCCCGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCg                                                                                        <  2:896731/138‑1 (MQ=255)
                                                       aaGGCGAGTTATAAAATGAGTTTTTTTGATGAGTTAAAACCCTCTCTGGAAGAGGCTGTCGAGATAAAACAAGGTTTAAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGaa                                                                                     <  1:904405/140‑1 (MQ=255)
                                                         ggCGAGGTATAAAATGATTTTTTTTGATGAGTTAAAACCCTCTCTGGAAGAGGCTGTCGAGATAAAACAAGTTTTGAAAAACCCTGCAGGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAAc                                                                                   <  1:838864/140‑1 (MQ=255)
                                                         ggCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAAc                                                                                   <  1:101615/140‑1 (MQ=255)
                                                          gCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACg                                                                                  <  2:266462/140‑1 (MQ=255)
                                                           gggggTAAAAAATGATTTTTTTTGATGAGTTAAAACCCTCTCTGGAAGAGGTTGTCGAGATTAAACAAGGTTTGAAAAACCCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACgg                                                                                 <  1:942355/137‑1 (MQ=255)
                                                           ggAGGTAAAAAATGATTTTTTTTGATGAGTTGAAACCCTCTCTGGAAGGGGCTGTCGAAATAAAACAAGGTTAAAAAAACCCTGCAGGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACgg                                                                                 <  2:1774171/139‑1 (MQ=255)
                                                             aGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCGACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGgtg                                                                               >  2:1144961/1‑140 (MQ=255)
                                                             aGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGgtg                                                                               >  2:1578011/1‑140 (MQ=255)
                                                                  taAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCAc                                                                          >  2:752579/1‑140 (MQ=255)
                                                                  taAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCAc                                                                          >  2:447953/1‑140 (MQ=255)
                                                                         aGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCAt                                                                   >  2:629992/1‑140 (MQ=255)
                                                                          gTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATc                                                                  >  1:917453/1‑140 (MQ=255)
                                                                              ttttGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGc                                                              >  2:1145296/1‑140 (MQ=255)
                                                                                 tgatgaGTTGAAACCCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAACCCTGCCGGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCAtt                                                           <  1:1074102/140‑1 (MQ=255)
                                                                                 tgatgaGTTAAAACCCTCTCGGGAAGAGGCTGTCGAGATAAAACAAGTTTTAAAAAACCCTGCACGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCAtt                                                           <  1:1075090/140‑1 (MQ=255)
                                                                                   atgaGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAAAAAGGTTTAAAAAACCCTGCACGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCg                                                         <  1:1522887/140‑1 (MQ=255)
                                                                                    tgaGTTAAAACCCTCTCGGGAAGGGGCTGTCGAAATAAAACAGGGTTGAAAAAACCCTGCAGGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCgg                                                        <  1:928060/140‑1 (MQ=255)
                                                                                    tgaGTTAAAACCCTCTCGGGAAGAGGCTGTGGAGATAAAACAGGTTTTAAAAAACCCTGCAGGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCgg                                                        <  1:357051/140‑1 (MQ=255)
                                                                                    tgaGTTAAAACCCTCTCGGGAAGAGGCTGTCGAGATAAAAAAAGGTTTAAAAAACCCTGCCGGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCgg                                                        <  2:612622/140‑1 (MQ=255)
                                                                                    ggaGTTGAAAGCGTCGGTAGAGGAGGCCGGCGAGATAAAAAAGGTTTTGAAAAGCCGTGCCGGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCgg                                                        <  2:1410902/139‑1 (MQ=255)
                                                                                        ttGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCt                                                    >  1:759513/1‑140 (MQ=255)
                                                                                                  ctctGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGa                                          >  2:857460/1‑140 (MQ=255)
                                                                                                  ctctGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGa                                          >  1:694273/1‑140 (MQ=255)
                                                                                                  ctctGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGa                                          >  2:1754430/1‑140 (MQ=255)
                                                                                                  ctctGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGa                                          >  1:1173558/1‑140 (MQ=255)
                                                                                                   tcggaaaGGGGATGTGGAAATAAAACAAGGTTAAAAAAACCATGCGGGGGGACCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGAt                                         <  2:195934/135‑1 (MQ=25)
                                                                                                    ggggAAGGGGCTGTGGAAATAAAACAGGGTTAAAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATg                                        <  2:88089/138‑1 (MQ=38)
                                                                                                       gAAGAGGCTGTCGAGATAAAACAGGGTTTGAAAAACCCTGCAGGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGa                                     <  1:1027009/140‑1 (MQ=255)
                                                                                                          gagGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGc                                  >  2:1708027/1‑140 (MQ=255)
                                                                                                            ggCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCtt                                >  1:195934/1‑140 (MQ=255)
                                                                                                             gCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTa                               >  1:1442348/1‑140 (MQ=255)
                                                                                                               tGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTAcc                             >  2:1653157/1‑140 (MQ=255)
                                                                                                               tGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTAcc                             >  1:50317/1‑140 (MQ=255)
                                                                                                                 tCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACcgc                           >  1:1131082/1‑140 (MQ=255)
                                                                                                                 tCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACcgc                           >  1:612511/1‑140 (MQ=255)
                                                                                                                   gaaataaaaCAGGGTTTGAAAAACCCTGCCCGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGt                         <  1:1626538/134‑1 (MQ=38)
                                                                                                                    gagatAAAAAAGGTTTGAAAAAACCCGCCCGGGTTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGt                         <  2:789162/136‑1 (MQ=38)
                                                                                                                      aTTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCtt                      >  1:369865/1‑140 (MQ=255)
                                                                                                                       ttAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTa                     >  2:990232/1‑140 (MQ=255)
                                                                                                                        tAAACAGGTTTTGAAAAACCCTGCCCGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTAt                    <  1:136521/140‑1 (MQ=38)
                                                                                                                        tAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTAt                    >  1:1304191/1‑140 (MQ=255)
                                                                                                                        aaaaaaaGGTTTAAAAAACCCTGCGCGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTAt                    <  1:1297102/139‑1 (MQ=255)
                                                                                                                          aaaaaGGTTTAAAAAAACCCTGCGGGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATc                   <  2:1173558/137‑1 (MQ=37)
                                                                                                                          aaaaaGGTTAAAAAAATCCTACCGGGTCGTCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCa                  <  1:548761/137‑1 (MQ=25)
                                                                                                                           aCAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAg                 >  1:1663400/1‑140 (MQ=255)
                                                                                                                              gggTTTAAAAAACCCTGCAGGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGcc              <  2:229783/139‑1 (MQ=255)
                                                                                                                              aGGTTTGAAAAACCCTGCACGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGcc              <  1:629992/140‑1 (MQ=255)
                                                                                                                              aGGTTTGAAAAACCCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGcc              <  1:457835/140‑1 (MQ=255)
                                                                                                                              aGGTTTAAAAAACCCTGCGGGGGAGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGcc              <  2:917453/140‑1 (MQ=37)
                                                                                                                               ggTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCt             >  2:1516069/1‑140 (MQ=255)
                                                                                                                               ggTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCt             >  2:528297/1‑140 (MQ=255)
                                                                                                                                ttttGAAAAACCCTGCAGGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTa            <  1:725813/139‑1 (MQ=255)
                                                                                                                                ttttGAAAAACCCTGCACGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTAAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTa            <  1:607956/139‑1 (MQ=255)
                                                                                                                                tttaaaaaaacCCGGCAGGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTa            <  2:520553/135‑1 (MQ=37)
                                                                                                                                 tttGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTAc           >  2:638053/1‑140 (MQ=255)
                                                                                                                                 tttGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTAc           >  2:621120/1‑140 (MQ=255)
                                                                                                                                          aCCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACa  <  1:528297/140‑1 (MQ=255)
                                                                                                                                          aCCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACa  <  1:1708027/140‑1 (MQ=255)
                                                                                                                                          aCCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACa  >  1:1229973/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                              
GGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACA  >  NC_000913/4628511‑4628788

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