Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,6640CG35.7% 25.1 / 18.7 43pseudogene (940/1617 nt)glvCpseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component
Reads supporting (aligned to +/- strand):  ref base C (13/14);  new base G (6/9);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGAT  >  NC_000913/3862427‑3862901
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
cAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCg                                                                                                                                                                                                                                                     <  2:615515/239‑1 (MQ=255)
        aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCACCGGGGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAAGTGTCGGGGATccg                                                                                                                                                                                                                                             >  2:262165/1‑237 (MQ=255)
        aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGATGCGGGGATcc                                                                                                                                                                                                                                              >  1:1768040/1‑237 (MQ=255)
        aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTGgggggaaca                                                                                                                                                                                                                                              >  1:252301/1‑233 (MQ=255)
        aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGCAATTCCCACCAGCATGGCGGTTAAGGGTGTGGGGAtacg                                                                                                                                                                                                                                             >  1:1797270/1‑236 (MQ=255)
                             gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACAAGCCTGGCGGTTTAATGTGCGTGAAACAACAAACGCACCACCTTTAccc                                                                                                                                                                                                                        >  2:1668528/1‑238 (MQ=255)
                              cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGGTGGGGGAGACAGCAAACCCGCCCCTTTTAccccc                                                                                                                                                                                                                      >  1:209113/1‑237 (MQ=255)
                              cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGGAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGTTGCGGGGAAAAGAAAACCCCGCACCTTTAc                                                                                                                                                                                                                          >  1:462016/1‑235 (MQ=255)
                                tGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTGGGGGGAACAGCAAAACCCCCACTTTTAccccaat                                                                                                                                                                                                                    >  1:1234967/1‑235 (MQ=255)
                                  gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGGGTCGGGGAGCAGGCAGCCCGCCCCCTTTAccccgatt                                                                                                                                                                                                                   >  2:696285/1‑233 (MQ=255)
                                     gCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTTCGGGAAATCACAAACACCCCACCTTTTACCCGATTTGcc                                                                                                                                                                                                               >  1:923975/1‑239 (MQ=255)
                                        gaaaggcgtgATACAGTGTCGCGACAAACCCGAGTCACCCAGCCCTCCCCCGATGGCCCCCACAACACAAAAGAGATACATTGGGGTGAAAATGGAGGCCGCCAGACGGGGGTGTCCCGAAACGAGAAGCGGAGAAAGAAACAGGAGGGTAAATCCCAGGGGTTCGGTATTCCCCACAAGAAGGGCGTTTAAGGTGGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTg                                                                                                                                                                                                           <  2:62530/230‑1 (MQ=255)
                                             tCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAATATGGCGGAATCAGAAAACCCGCCACATTTACCCCACTTTcagtctcggcg                                                                                                                                                                                                      >  1:2257592/1‑229 (MQ=255)
                                                     catttagcgggTAACCCCTATCCAACCAGCCCCCCCCCCATGTTCCCCACCACCACAAAGAAATAAATCGCGGTCGTGATAGGGGCGCCCAGCACGGTTGTCACCGCAAACAGAAACGGAGAAATAAAAGAGAAGGAAAATTCCGGGGTTTCGAAATTCCCCCCCAGCAGGGTGGTTAGGGTGGGGGGATTAAACAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTa                                                                                                                                                                                             <  2:234818/230‑1 (MQ=255)
                                                                                                cgCCCACAGCAACAAAGAGATATATGGCGGTCGACATAGGGGCCGCCTGCAGGGGTTGTCCGGAAAACAGCGAGGGGAAAAGGAGAGGGAAGTAAAACCCCGGGGTTTGGGTATTCCCCAGCAGCGTGGCGGTTAGGGTGGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACAACGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTa                                                                                                                                               <  2:2174367/238‑1 (MQ=255)
                                                                                                       cacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGATGCGGGGAAAAGGAAAACCCCCCCCTTTTCCCCCGTTTCCGGTGCGGGACGGATAGAAATCGCTCAAAAAAAAGCCACGGCGCCCAAAGTTTTTTAGTTatagcgcg                                                                                                                                        >  2:1087858/1‑232 (MQ=255)
                                                                                                       cacaccaAAGAAAAACATTCCGGGAAAAATTGAGCCCGCCAGCACGGCGTATACCGAAAACAGCAACGGTAAAATGAACAGGAAGGTAAATCCCAGCGGTTCGGTAATCCCCACCAGCAGGGGGGTTAGGGTTGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGc                                                                                                                                                                                                          <  2:806283/177‑1 (MQ=255)
                                                                                                                 agaTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGGGGGGAACAGCAAACCCGCCACCTTTAACCGCATTTTCGGGGCGGGCGTGAAg                                                                                                                                                                                               >  1:806283/1‑177 (MQ=255)
                                                                                                                                     tGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGATAAGGTTGCGGGGAACAGCAAAACCGCCCCTTTTAACCCATTTTTCg                                                                                                                                                                                                              >  1:1838856/1‑143 (MQ=255)
                                                                                                                                        ggCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAGGGTTGGGGGATTAACAAACCCCCCCCCTTTTACCGAATTTTc                                                                                                                                                                                                                >  1:1204193/1‑139 (MQ=255)
                                                                                                                                            gccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTTCGGGGAATAACAAAACCCCCCCCTTTAACCCCTTTTTCCGTGCGGCCGGGGAAGACAACGCCAAAAAAAAAGCCCCCGGCCCCAAAAATTTTTAGTTAACCGTGCGGGGACAACCGCCCTCCCGGGAAAACGAATTTCAGCg                                                                                                   >  2:796116/1‑239 (MQ=255)
                                                                                                                                            gccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTTCGGGAAAAAACAAACCCGCCCCCTTTACCCGACTTTTCGGGGGGGGCGGGGAGTAAAATGCCTAAAAAAAGGCCCCCGCCCCCAAAATTTTTTGTGTTACGTGGAGGGGAAAAACCCCCTCCGGGGAAAACCACTTTAGccg                                                                                                  >  2:2072152/1‑237 (MQ=255)
                                                                                                                                            gccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTGCGGGGAACAGCAAACCCGGCAATTTTTTCCGAGTTTTCCGCGCGGGAGGGAAGGAAAAAGGGCAAAAAAAAGCCCAGGCACCAAAAAGTTTTTGATTTACGGGGCAGGGGAAAACGCCTTCCCGGGACCACgcattttgca                                                                                                    >  2:1966635/1‑229 (MQ=255)
                                                                                                                                                                 aaaCAGAAAGGGTGAAAGGAACAGGAAGGTAATTCCCAGCGTTGCGGTAATCCCCCCCAGCAGGGCGGTTAGGGTGGCGGGAACCAGAAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGTGGCTCGGCGCTCAAACTg                                                                                  <  1:1261233/235‑1 (MQ=255)
                                                                                                                                                                                  tGAACAGGAAGGTAAATTCCGGCGGTTCGGTATTCCCCCCCAGCTGGGGGGTTAGTGTGGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                             <  1:1259036/239‑1 (MQ=255)
                                                                                                                                                                                  agAACAGGAAGGTAAATTCCAGGGGTGCGGTAATCCCCACCAGCAGGGCGGTTAAGGTTGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                             <  1:262165/238‑1 (MQ=255)
                                                                                                                                                                                        gggaattAAATCCCACGCGGTGCTGTCATCCCCACCAGCATGTGTGTTGAGTTGGCGGAAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                  <  2:252301/223‑1 (MQ=255)
                                                                                                                                                                                         gAAGGTAAATTCCGGGGGTTCCGATATCCCCCACGAGCTGGGCGTTTAGGGTTGGGGAAACAACCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                             <  2:2298484/233‑1 (MQ=255)
                                                                                                                                                                                          aagaaataTTCCCGCGGTTCGTTTATTCCCCCAAGCAGGGCGGTGAGGGTTGCGGTAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                  <  1:1401219/219‑1 (MQ=255)
                                                                                                                                                                                           aagTAAATTTCAGGGGTTTCGGTAACCCCCCCCGGCAGGGCGTTTAGGGTGGCGGAAAGACGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                             <  1:2054818/229‑1 (MQ=255)
                                                                                                                                                                                                   tccccacGGTTCGGGAAATTCCCACGACCTTGGCGGTTAGGGGTCCGGGAACCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATg                                              <  1:1504087/233‑1 (MQ=255)
                                                                                                                                                                                                               ttggttAATCCCCCCCAGCGTGGCGGTTAAGGTTGGGGAAACCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTa                                                                                                                                               <  1:2207401/128‑1 (MQ=255)
                                                                                                                                                                                                                    tAATCCCCACCAGCAGGGCGGTTAAGGTGGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCaa                                                                                                                                                               <  2:1099261/114‑1 (MQ=255)
                                                                                                                                                                                                                          ccacaaGCAGGGCGGTTAAGTTGGCGGGAACAAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGaa                                                                                <  1:1773244/184‑1 (MQ=255)
                                                                                                                                                                                                                           ccgcgatgtGGGGGTTGGGGTCGCGAGAACCAGAGACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                        <  2:1234967/231‑1 (MQ=255)
                                                                                                                                                                                                                           ccccaGCAGGGCGTTTAGGGTTGCGGGAAACAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                        <  1:1583587/237‑1 (MQ=255)
                                                                                                                                                                                                                           caccaGCAGGGCGGTTAAGGTGGCGGGAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                        <  1:796116/239‑1 (MQ=255)
                                                                                                                                                                                                                           caccaGCAGGGCGGTAATGGTGGCGGAAATCAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCCCCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                        <  1:2201794/239‑1 (MQ=255)
                                                                                                                                                                                                                              caGCATGGCGGTTAAGGTTGCGGGGAACAACAAACCCCCCCCCTTTACCCGGTTTTCCGGGGCGGCAGTGAAGTACACCGCTAAAGAAAATCCCACGGCACCAAAAAttttttg                                                                                                                                                    >  1:1099261/1‑112 (MQ=255)
                                                                                                                                                                                                                                      cGGTAAAGGTAGCGGAACACAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa              <  1:1668528/238‑1 (MQ=255)
                                                                                                                                                                                                                                      cGGATAAGGTTGGGGGAACCGGAAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa              <  1:1087858/238‑1 (MQ=255)
                                                                                                                                                                                                                                             agtgtggCGGAAATCAGCAGCCCCGCCACTTTTCCCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTc         <  1:1933047/232‑1 (MQ=255)
                                                                                                                                                                                                                                                   gCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCaa                                    >  1:1997167/1‑205 (MQ=255)
                                                                                                                                                                                                                                                   gCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGAt  >  1:1971669/1‑239 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
CAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGAT  >  NC_000913/3862427‑3862901

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