Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,628,6490AC23.9% 173.9 / 11.8 71intergenic (+102/+206)nadR/ettADNA‑binding transcriptional repressor/NMN adenylyltransferase NadR/energy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base A (44/10);  new base C (0/17);  total (44/27)
Fisher's exact test for biased strand distribution p-value = 8.14e-10
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:546055/140‑1 (MQ=255)
   aCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTg                                                                                                                                          <  1:1184477/140‑1 (MQ=255)
     tGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCa                                                                                                                                       >  2:1107044/1‑140 (MQ=255)
      ggTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCAc                                                                                                                                      <  1:945073/140‑1 (MQ=255)
      ggTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCAc                                                                                                                                      <  1:919221/140‑1 (MQ=255)
       gTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGTTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTg                                                                                                                                          >  2:515564/1‑136 (MQ=255)
       gTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGTTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTg                                                                                                                                          <  1:515564/136‑1 (MQ=255)
        tGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACgg                                                                                                                                    >  1:1616952/1‑140 (MQ=255)
           gggAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTg                                                                                                                                 >  1:369922/1‑140 (MQ=255)
              agaTGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAACGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAACCCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAACCCTGCACGGGTGAcc                                                                                                                              <  2:987/140‑1 (MQ=255)
                atgatgGGGGAGCAGAGATACCCGTGATAAAAGTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTGGAAAACCTCTCGGGAAGGGGCTGTCGAGTTAAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCg                                                                                                                            <  2:793450/140‑1 (MQ=255)
                     gggggAGAAGAGATACCCGTGATGAACCTCCTCAAAGGCGAGGTAAAAAATGATTTTTTTTGATGATTTGAAACCCTCTCTGGAAGAGGCTGTCGAGATTAAACAGGTTTTGAAAAACCCTGCACGGGTGACCCGCCACg                                                                                                                       <  2:599943/140‑1 (MQ=255)
                        ggAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAACCCTCTCTGGAAGAGGCTGTCGAGATTAAACAGGGTTTGAAAAAACCTGCAGGGGTGCCCCGCCACGaaa                                                                                                                    <  2:468788/140‑1 (MQ=255)
                          aGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAAtt                                                                                                                  >  2:559504/1‑140 (MQ=255)
                          aGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTAAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAAtt                                                                                                                  >  1:1059622/1‑140 (MQ=255)
                            cAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGa                                                                                                                >  2:273346/1‑140 (MQ=255)
                            cAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGa                                                                                                                >  1:1120671/1‑140 (MQ=255)
                             agagATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAACCCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAACCCTGCACGGGTGACCCGCCACGAAATTGAg                                                                                                               <  1:1544831/140‑1 (MQ=255)
                             agagATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAg                                                                                                               <  1:1181789/140‑1 (MQ=255)
                              aagaTAACCGTGATGAACCTGCTAAAAGGCGAGGTATAAAATGATTTTTTTTGATGATTTAAAACCCTCTCTGGAAGAGGCTGTCGAGATAAAACAAGGTTTAAAAAACCCTGCAGGGGTGACCCGCCACGAAATTGAgg                                                                                                              <  2:602186/139‑1 (MQ=255)
                                 aTAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATg                                                                                                           >  2:1442474/1‑140 (MQ=255)
                                   aaCCGTGATAAAACTGCTCAAGGGCGAGGTAAAAAATGAGTTTTTTTGATGAGTTGAAACCCTCTCTGGAAGAGGCTGTCGAGATTAAACAGGTTTTAAAAAACCCTGCACGGGTGCCCCGCCACGAAATTGAGGATGCt                                                                                                         <  2:369922/140‑1 (MQ=255)
                                       gTGATGAAACTGCTCAAAGGCGAGGTATAAAATGATTTTTTTTGATGAGTTGAAACCCTCTCTGGAAGAGGCTGTCGAGATAAAACAAGGTTTGAAAAACCCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAgg                                                                                                     <  1:533097/140‑1 (MQ=255)
                                         gatgaAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCt                                                                                                   >  2:480302/1‑140 (MQ=255)
                                         gatgaAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCt                                                                                                   >  1:822972/1‑140 (MQ=255)
                                         gatgaAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCt                                                                                                   >  2:723004/1‑140 (MQ=255)
                                         gatgaAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGTGACCCGCCACGAAATTGAGGATGCTAAGGCtg                                                                                                  >  2:789617/1‑140 (MQ=255)
                                          atgaAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCtg                                                                                                  >  2:381098/1‑140 (MQ=255)
                                           tgaAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCtgt                                                                                                 >  1:305549/1‑140 (MQ=255)
                                           tgaAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCtgt                                                                                                 >  1:926999/1‑140 (MQ=255)
                                              accTGCTAAAAGGCGAGGTATAAAATGATTTTTTTTGATGAGTTAAAACCCTCTCTGGAAGAGGCTGGCGAGATTAAACAAGGTTTAAAAAACCCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCtgttgt                                                                                              <  2:683930/138‑1 (MQ=255)
                                                  gCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGAc                                                                                          >  2:310141/1‑140 (MQ=255)
                                                   cTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGAcc                                                                                         >  1:1423300/1‑140 (MQ=255)
                                                   cTCAAAGGCGAGGCATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGAcc                                                                                         >  1:1114249/1‑140 (MQ=255)
                                                        aGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGaaa                                                                                    >  1:1108980/1‑140 (MQ=255)
                                                                  taAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCAc                                                                          >  1:712248/1‑140 (MQ=255)
                                                                   aaaaTGATTTTTTTTGATGAGTTAAAACCCTCTCGGGAAGAGGCTGTCGAAATAAAAAAAGTTTTAAAAAACCCTGCAGGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCAcg                                                                         <  2:787916/140‑1 (MQ=255)
                                                                         aGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCAt                                                                   >  1:757660/1‑140 (MQ=255)
                                                                         aGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCAt                                                                   >  1:219011/1‑140 (MQ=255)
                                                                          gTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATc                                                                  >  1:1602857/1‑140 (MQ=255)
                                                                          gTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATc                                                                  >  2:159435/1‑140 (MQ=255)
                                                                          gTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATc                                                                  >  2:398285/1‑140 (MQ=255)
                                                                          gTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATc                                                                  >  1:787916/1‑140 (MQ=255)
                                                                          gTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATc                                                                  >  2:192078/1‑140 (MQ=255)
                                                                          gTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCCGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATc                                                                  >  1:396289/1‑140 (MQ=255)
                                                                               tttGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCa                                                             >  1:1193047/1‑140 (MQ=255)
                                                                                ttGATGAGTTGAAACCCTCTCGGAAAGAGGCTGTCAAAATAAAACAAGTTTTAAAAAACCCGGCCGGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCAt                                                            <  2:305549/140‑1 (MQ=255)
                                                                                 tgatgaGTTAAAACCCTCTCGGGAAGGGGCTGTCGAAATAAAACAAGGTTTAAAAAACCCTGCACGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCAtt                                                           <  2:1423300/140‑1 (MQ=255)
                                                                                    tgaGTTGAAACCCTCTTTGGAGGAGGCTTTCGAGATAAAACAGGTTTGAAAAACCCCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCgg                                                        <  2:747994/140‑1 (MQ=255)
                                                                                    tgaGTAAAAACCCTCTGTGGAAGAGGCTGTCGAAATAAAACAGGTTTTAAAAAACCCTGCGGGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCgg                                                        <  2:1386720/140‑1 (MQ=255)
                                                                                    ggaGTTAAAACCCTCTCGGAAAGGGGCTGTCGAGATAAAACAGGTTTTAAAAACCCCTGCGGGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCgg                                                        <  2:822972/139‑1 (MQ=255)
                                                                                     gaGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGt                                                       >  1:1255496/1‑140 (MQ=255)
                                                                                         aaaaacCCTCTCTGGAAGAGGCTGTCGAGATAAAACAGGTTTTAAAAAACCCTGCGGGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTc                                                   <  2:47008/138‑1 (MQ=255)
                                                                                            cagccTCTCGGGAAGAGGCTGTCGAGATAAAACAAGTTTTAAAAAACCCTGCAGGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAAt                                                <  2:1338908/137‑1 (MQ=255)
                                                                                                      ggacagagCTGTCAAGGTTAAAACATGGTTTAAAAAAACGTGCGGGGGCGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGc                                       <  2:1114249/133‑1 (MQ=25)
                                                                                                      ggAAGAGGCTGTCGAGATTAAACAAGGTTTAAAAAACCCTGCGGGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCg                                      <  1:480302/140‑1 (MQ=255)
                                                                                                        aaGGGGCTGTCGAGATAAAAAAAGGTTTAAAAAACCCTGCAGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACg                                   <  1:789617/140‑1 (MQ=38)
                                                                                                          gagGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGActg                                  >  1:1338908/1‑138 (MQ=255)
                                                                                                          gagGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGc                                  >  2:1517/1‑140 (MQ=255)
                                                                                                           gggCTGGCGATATGAAACAGGGTGAGAAAACCCCTGCCGGGGGGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCt                                 <  2:1616952/139‑1 (MQ=37)
                                                                                                            ggCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCtt                                >  2:1157399/1‑140 (MQ=255)
                                                                                                             gCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTa                               >  2:231816/1‑140 (MQ=255)
                                                                                                              gtgtgGAAATAAAACAGGGTTTAAAAAACCCTGCGGGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTAc                              <  2:1602857/139‑1 (MQ=37)
                                                                                                              cTGTCGAGATAAAACAAGGTTTAAAAAACCCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTAc                              <  1:1550584/140‑1 (MQ=255)
                                                                                                              cTGCCGAAATAAAACAGGTTTGGAAAAACCCTGCACGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTAc                              <  1:1583523/140‑1 (MQ=38)
                                                                                                               tGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTAcc                             >  2:1401800/1‑140 (MQ=255)
                                                                                                               ctttcgaaTTTAAAAAGGTTTGAAAAAACCTGCCCCGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTAc                              <  2:219011/133‑1 (MQ=37)
                                                                                                                 tCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACcgc                           <  1:540866/140‑1 (MQ=255)
                                                                                                                 tCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACcgc                           >  1:1269402/1‑140 (MQ=255)
                                                                                                                 tCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACcgc                           >  2:1575028/1‑140 (MQ=255)
                                                                                                                 tCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACcgc                           >  2:133866/1‑140 (MQ=255)
                                                                                                                  cGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACcgcg                          >  1:929764/1‑140 (MQ=255)
                                                                                                                   gagaTAAAACAAGTTTTGAAAAACCCTGCGCGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGt                         <  1:1693617/140‑1 (MQ=255)
                                                                                                                    gaaataaaaCAAGGTTTAAAAAACCCTGCACGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGt                         <  1:753280/134‑1 (MQ=38)
                                                                                                                      aTTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCtt                      >  2:1281020/1‑140 (MQ=255)
                                                                                                                         aaaaCAAGTTTTGAAAAACCCTGCACGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTAt                    <  1:434040/139‑1 (MQ=255)
                                                                                                                                ggTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCt             >  1:1218039/1‑140 (MQ=255)
                                                                                                                                ggTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCt             >  2:1571959/1‑140 (MQ=255)
                                                                                                                                 ttttAAAAAACCCTGCCGGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTa            <  1:723004/139‑1 (MQ=37)
                                                                                                                                 gTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTa            <  1:1401800/140‑1 (MQ=255)
                                                                                                                                 gTTTAAAAAACCCTGCACGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTa            <  1:1442474/140‑1 (MQ=255)
                                                                                                                                     gAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGt        >  1:126097/1‑140 (MQ=255)
                                                                                                                                            aCCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTa            >  1:738542/1‑130 (MQ=255)
                                                                                                                                            aCCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTa            <  2:738542/130‑1 (MQ=255)
                                                                                                                                            aCCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACa  <  1:1571959/140‑1 (MQ=255)
                                                                                                                                            aCCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACa  <  2:1193047/140‑1 (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: