Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,745,0190CG43.8% 17.3 / 26.8 48G72G (GGC→GGGyiaN2,3‑diketo‑L‑gulonate TRAP transporter large permease protein
Reads supporting (aligned to +/- strand):  ref base C (7/20);  new base G (9/12);  total (16/32)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGC  >  NC_000913/3744783‑3745245
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                   
gCCTCCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCGCTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATgcgggcgg                                                                                                                                                                                                                                   <  1:818867/239‑1 (MQ=255)
gCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGAATgggggcgg                                                                                                                                                                                                                                   >  1:946984/1‑239 (MQ=255)
        aaGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGCGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCaa                                                                                                                                                                                                                          <  1:768915/239‑1 (MQ=255)
               ttAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCAt                                                                                                                                                                                                                   <  1:454309/239‑1 (MQ=255)
                      tGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATTCGGGCGGCTTGTCAAAACGCAttgttta                                                                                                                                                                                                             >  1:1228771/1‑236 (MQ=255)
                      tGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACATTGATTCGGGGGGCCTTGTCAAACCGcttg                                                                                                                                                                                                                  >  2:794543/1‑231 (MQ=255)
                      tGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGAATGCGGGGGGCTTGTCCAAACGCATTGTTGa                                                                                                                                                                                                             >  1:180324/1‑238 (MQ=255)
                         cTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGGGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGGGGCTTGTCAAAACGCATTGTTGACCt                                                                                                                                                                                                          <  1:794543/238‑1 (MQ=255)
                         cTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCGGGACATGTTTGATGCCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCt                                                                                                                                                                                                          <  2:391720/238‑1 (MQ=255)
                           gtgCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATCATGCATGCGGGGGGCCTGTTCAAACGCATTTTTTACCttc                                                                                                                                                                                                        >  1:871643/1‑236 (MQ=255)
                                        tgggtggtctgtTGGGTGGTATGGCTATCGGTTTGCTTTTTGCCGGGGCTCTGTGGTTGTGCGGGGGGGCGTTAATGTTCGGGCGGGACGTGTTTGATGTCCAGATCATGGCACAACCGCTGGTGAGCGGAGCCGATAGCTTCTCCCTGCGGGCGATTCCGTTCTTTGTGCTGGGGGGTGAAATCATGAATGCGGGCGGTTTGTAAAAACGCATTGTTGACCTGCCGATGAAACtggtg                                                                                                                                                                                          <  2:1338211/233‑1 (MQ=255)
                                               tctgtTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGAATTCGGGGGGGTTGTCAAAACGCCCTGTTGTCCTGCCGGCGAAACTGGGGGGGCCTa                                                                                                                                                                                   >  1:310541/1‑239 (MQ=255)
                                                 gtgggggTGGTGTTGCTTGTGGTTGTCCGTTTGCCGGGGCCTTGTTGTTGTGGGGGGCGGGGTTAATGTTCTGGCTGGACATGTTGGAGGTCCAATTCAGGGCAAAAACGTTGGTGAACGGGGCAGATAGCTTCTCCCGGCGGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGAATGCGGGGGGCTTGTAAAAACGAATTGTTGACCTGCCGATGAAACtggt                                                                                                                                                                                           <  2:1281082/225‑1 (MQ=255)
                                                   ttGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCt                                                                                                                                                                                                                                <  1:499809/190‑1 (MQ=255)
                                                      ggtggtATTGCTCTCGGTTTGCGTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCGTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAAGGCGGGCGGCTTGTCa                                                                                                                                                                                                                           <  1:1240614/192‑1 (MQ=255)
                                                           tatggatATCGGTTTGCCTTGGGCCTGGGCTCGTGTTTTGTGCGGGGCGGGGTTAGTGTTGGGGCGGGACATTTTTGATGTCCAAGTCAGGGCACAACCGTGGGGGAGCGGGGCCGATAGCTTCTCCCGGCGGGGGATCCGGTTTTTTGTGCTGGCGGGGAAAATCAGGAATGCGGGGGGTTTGTCAAAACGCATTGTGGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCggg                                                                                                                                                                        <  2:1234661/232‑1 (MQ=255)
                                                           tATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCGGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTAAAAACGCATTGTTGGCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCggg                                                                                                                                                                        <  1:557543/238‑1 (MQ=255)
                                                           tATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACAAGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCCTGAATTCGGGCGGCTTGTCAAAACGCATTGTTTACCTGCCGGTTAAACATGTgggggctaca                                                                                                                                                                                 >  1:1281082/1‑223 (MQ=255)
                                                                    cGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATTCTGAATTCGGGCGGCCTGTCCACACCCCTTGTTTACCCGCCGAGTGAACTGGTGGGGGATAAACACGGCgggccggggtcc                                                                                                                                                               >  2:1335160/1‑234 (MQ=255)
                                                                     ggTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACATGAATTCGGGGGGGTTTGCCAAAAGCACTTTTTTGCTTGCGCTGTAACTGGGGGGGGCGAAAACAGGCggggtgggggtcc                                                                                                                                                              >  2:963882/1‑233 (MQ=255)
                                                                                tGCCTGGGCACTTTTGTGGTGCGGGGGGGCGTTATTGTTCTGGCTGGACATGTTTGTTGTCGAGATCGTGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGTTTCCGTTCTTTGTGCGGGCGGGTGAAATCATGAATGCGGGCGGTTTGTCAAAACGATTTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTgg                                                                                                                                                  <  2:946984/239‑1 (MQ=255)
                                                                                 gCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCTCAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGGTTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCa                                                                                                                                                                                                                           <  1:528530/165‑1 (MQ=255)
                                                                                    tGGGCACTGTGTTTGTGGGGGGCGGCTTTAATTTTCTGGCTGGACATTTTTGATGTCCAGATCATGGCACAACCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGc                                                                                                                                              <  1:954289/239‑1 (MQ=255)
                                                                                            tgttgttggGGGGGGGGGCTTTAATGTTCGGGCGGGACATGTTTGAGGCCCAGATCGTGGCACAAACGCTGGTGAACGGGGCCGATAGCTTCTCCCTGCTGGCGTTTCCGTTTTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGTTTGTCAAAACGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTgg                                                                                                                                                  <  2:776777/228‑1 (MQ=255)
                                                                                                    ggcggggcgggTTTATTGTTTGGGCTGGCCTTGTTTGTGGTCCAGATCATGCCACAAGGGCTGGGGAAGGGGGCTGCTAGTTTGTCCCTGCGGGCGATTCCGTTCTTTGGGCGGGCGGGGGAAATCGGGAATGGGGGGGGGTGGTAAAAACGCTTTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAg                                                                                                                               <  2:606829/238‑1 (MQ=255)
                                                                                                      cGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGAATTGGGGGGGCCTGTCAAAACGCATTGTTTACCTGCCGATGAAACTGGGGGGGGATAAAACAGGCGGGCTGGGCTACGCCGGCGTGCTGGGCGCTATGAt                                                                                                                                        >  1:776777/1‑228 (MQ=255)
                                                                                                        ggggtgggTTAATGTTCGGGCGTGACATGTTTGTTTTCCAGTTCATGGCCAAACCGCTGGTGACCGGGGCCGGTAGCTTCTCCCTGCTGGCGTTTCCGTTCTTTGTGCTGGCGGGTGAAATCATGATTGCGGGGGGTTTGTCAAAACGCTTTGTTGCCCTGCCGATGAAACtggt                                                                                                                                                                                           <  2:104008/168‑1 (MQ=255)
                                                                                                                  aaTGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCAAGAATGCGGGGGGCTTGGCAAAACGCATTGTTTAACTGCCGATGAAAATGGTGGGGGATaaa                                                                                                                                                                                 >  1:104008/1‑175 (MQ=255)
                                                                                                                           gctggACATGTTTGGTGTCCTGATCATGGCACAAACGCTGGTGACCGGAGCCGATAGCTTCTCCCTGCTGGCGTTTCCGTTCTTTGTGCTGGCGGGTGAAATCATAAATGCGGGCGGCTTGTCAAACCGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTa                                                                                                                                                                <  2:665362/183‑1 (MQ=255)
                                                                                                                                     tttGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAAACATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCCGGGGTACGgcggcggg                                                                                                                                                      >  1:665362/1‑181 (MQ=255)
                                                                                                                                            tCCAGATAATGGCACAACCGTTGTTGAAGGGGGCCGATAGCTTCCCCTTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTAAAATCATGAATGCGGGGGGCTTGTCAAAACGCTTTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCt                                                                                                                                             <  2:837984/185‑1 (MQ=255)
                                                                                                                                              cAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATTATGAATGCGGGCGGCTTGGCAAAAAGCACTTTTGTGCTCCCGGCGAAAATGGGGGGGGCTaaa                                                                                                                                                                                 >  1:33845/1‑147 (MQ=255)
                                                                                                                                                       gCACAACCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTgg                                                                             <  1:963882/238‑1 (MQ=255)
                                                                                                                                                         acaAGCGCTGGTGAACGGAGCCGATTGCTTCCCCCGGCTGGCGTTTCCGTTTTTTGTGCGGGGGGGTAAAATCATGAATGGGGGGGGCTTGTCAAAACGCTTTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGc                                                                          <  1:240723/239‑1 (MQ=255)
                                                                                                                                                           gataacgggTGAACGGGGCCGATAGCTTCTCCCTGCGGGCGATTCCGTTTTTGGTGCGGGCGGGTGAAATAATAAAGGCGGGCGGCTTGTAAAACCGCATTGTGGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTg                                                                                                                                                      <  2:36200/154‑1 (MQ=255)
                                                                                                                                                                     gAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAAACATTAATGCGGGGGGCTTGTCAAAACGCATTGTTTACCTGCCGGTGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGGGCTGGGAGCTa                                                                                                                                            >  1:36200/1‑161 (MQ=255)
                                                                                                                                                                          gggCCGATAGCTTCCCCCTGCTGGCGATTCCGTTCTTTGTGCTGGGGGGTGAAACCATGAATGCGGGGGGCTTGTCAAAACGCATTTTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAgcgc                                                          <  1:993941/236‑1 (MQ=255)
                                                                                                                                                                           agcggttggCTTCTCCCTGCGGGCAATTCCGTTTTTTGTGTGGGGGGGTAAAATCAGGAATGCGGGGGGCTTGTAAAAACGCATTGTTGGCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAgc                                                            <  1:462104/227‑1 (MQ=255)
                                                                                                                                                                             ccGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATAccgccgc                                                                                              <  1:141728/199‑1 (MQ=255)
                                                                                                                                                                                   gtttCTCCCTGCGGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGTTTGCAAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCg                                                                                                                                                           <  2:741863/130‑1 (MQ=255)
                                                                                                                                                                                    cTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCAGGAAGGCGGGCGGTTTGTAAAAACGCTTTGTGGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCgggc                                                                                                                                                                       <  2:394246/119‑1 (MQ=255)
                                                                                                                                                                                     ttCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGc                                                                                                                                                                                                        <  1:984149/85‑1 (MQ=255)
                                                                                                                                                                                             gctgGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATCATGAATTCGGGCGGCTTTTCAAAACACATTGTTGACCTGCCGGTGAAAATGGGGGGGGCTAAAACAGGCGGGGTGGGCTACGTTGGCGGGCttgg                                                                                                                                                 >  1:741863/1‑129 (MQ=255)
                                                                                                                                                                                              ctgGCGTTTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGGGGCTTGTCAAAACGCATTGTTGACCTGCCGAt                                                                                                                                                                                                    <  1:666470/80‑1 (MQ=255)
                                                                                                                                                                                              ctgGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCCTGAATGCGGGCGGCTTTTCCAAACGCCTTGTTGACCTGCCGGTGAAACTGGGGGGGGATAAACCAGGCGGGCTGGGGTACgg                                                                                                                                                             >  1:394246/1‑118 (MQ=255)
                                                                                                                                                                                               tgGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAAAATGAATGGGGGGGGCTTTTCAAAACGCATTGTTGGCCTGCCGATGaaaat                                                                                                                                                                                              >  2:984149/1‑83 (MQ=255)
                                                                                                                                                                                                   gATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGTTTGTCAAAACGCTTTTTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGGGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGgcggcg                                 <  1:833124/238‑1 (MQ=255)
                                                                                                                                                                                                                             aaaTCAGGAAGGCGGGGGGTTTGTAAAACCGCATTGTTACCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTAt       <  1:1069062/238‑1 (MQ=255)
                                                                                                                                                                                                                                cggagagtGCGGGGGGTTTGTAAAAGCGGTTGTTTGCCCTGCGGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATcgcg   <  2:180324/232‑1 (MQ=255)
                                                                                                                                                                                                                                 cATGAAGGGGGGCGGTTTGTCAAACCGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCgcgc  <  1:155206/239‑1 (MQ=255)
                                                                                                                                                                                                                                             cggcTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAgg                                                                                                                                                                            >  2:683163/1‑57 (MQ=255)
                                                                                                                                                                                                                                             cggcTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAgg                                                                                                                                                                            >  2:652678/1‑57 (MQ=255)
                                                                                                                                                                                                                                             cggcTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAgg                                                                                                                                                                            <  1:683163/57‑1 (MQ=255)
                                                                                                                                                                                                                                             cggcTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAgg                                                                                                                                                                            <  1:652678/57‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                   
GCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGC  >  NC_000913/3744783‑3745245

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: