Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,6770GA34.2% 8.2 / 17.9 47pseudogene (927/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 G (6/19);  new base A (13/0);  total (20/27)
Fisher's exact test for biased strand distribution p-value = 1.00e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.75e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTG  >  NC_000913/3862456‑3862910
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                
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)
                                                                                                                                                                                                                agtgttaggggAATAACCAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTcc                 <  1:2072152/231‑1 (MQ=255)
                                                                                                                                                                                                                agtgtggCGGAAATCAGCAGCCCCGCCACTTTTCCCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTc                  <  1:1933047/232‑1 (MQ=255)
                                                                                                                                                                                                                agggggggggaaaTCAGGAACGCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACATCTGCCGGACCAAAGATAAACTGTcc                 <  1:295286/228‑1 (MQ=255)
                                                                                                                                                                                                                  gtgggggggAACAAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAAc                                                                                             <  2:940041/155‑1 (MQ=255)
                                                                                                                                                                                                                      gggggAACCAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGa            <  1:1380092/236‑1 (MQ=255)
                                                                                                                                                                                                                      gCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCaa                                             >  1:1997167/1‑205 (MQ=255)
                                                                                                                                                                                                                      gCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGAt           >  1:1971669/1‑239 (MQ=255)
                                                                                                                                                                                                                        gggAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGa         >  1:1854644/1‑239 (MQ=255)
                                                                                                                                                                                                                          gAACCAACAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGaa        <  1:1290091/238‑1 (MQ=255)
                                                                                                                                                                                                                            aTCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTg                                                                                   >  1:940041/1‑161 (MQ=255)
                                                                                                                                                                                                                              cAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAACGATAAACTGTCCGTAGATGTAgggg  >  1:1078049/1‑238 (MQ=255)
                                                                                                                                                                                                                                ggacgCCCGCCACTTTTACCCGATTGTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCACACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAgtggtg  <  2:748894/234‑1 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                
GCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTG  >  NC_000913/3862456‑3862910

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