Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,6640CG39.3% 30.4 / 26.3 62pseudogene (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 (22/15);  new base G (6/18);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 9.96e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGT  >  NC_000913/3862428‑3862897
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  
aaTGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTATAGTTgggg                                                                                                                                                                                                                                                  >  1:2184545/1‑239 (MQ=255)
  tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAAGGTGCgggg                                                                                                                                                                                                                                                >  1:1085342/1‑238 (MQ=255)
  tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCgggg                                                                                                                                                                                                                                                >  1:3493366/1‑238 (MQ=255)
       aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTGGGGGGaaca                                                                                                                                                                                                                                             >  1:1975488/1‑236 (MQ=255)
       aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTGAAGGGTGCggggagta                                                                                                                                                                                                                                            >  1:2886037/1‑233 (MQ=255)
             gTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGATAAGGGTTCGGGGATcaacaaag                                                                                                                                                                                                                                   >  1:3426405/1‑236 (MQ=255)
                  catcatATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGTTGGGGGGAACAACAAACCCGcc                                                                                                                                                                                                                              >  1:1924177/1‑238 (MQ=255)
                             cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTTCGGGGAACAGCAATCCCCCCACATTTTCCCGa                                                                                                                                                                                                                  >  1:800336/1‑239 (MQ=255)
                                        aaCATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTGCGGGGAACAACAAAACCCCCCCTTTTTCCCCGCTTTTCCgtgg                                                                                                                                                                                                        >  1:3400569/1‑236 (MQ=255)
                                                       ttttGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGGGGGGACAAACAAAGCCCCCCCCTTTTCCCCATTTTTCCCGTGGGGCAGGGAGTACa                                                                                                                                                                                        >  1:1114516/1‑238 (MQ=255)
                                                       tgtggggacaaaaCCTAATAAACCAGCCCTCCGCGCGTGCTCCCCCCCACAAAAAAGAAATACACTACGCGCGAAAGGGAGCCCCCCAACACGGCGTGCACCGAAAACGGCAGCGGTAAATGGAACAGGAGGTTAAATTCAAGAGTTTCGGTAATCCCCACTGGGGGGGTGTTTAAGGTGGGGGGAATAACCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACAt                                                                                                                                                                                       <  2:2468702/230‑1 (MQ=255)
                                                          tGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGAGTTTGGGAAGAAGCACACCCCCCACTTTTAACCCCTTTTTCGGTGCGGCAGAGAAGGACAtag                                                                                                                                                                                     >  2:1769382/1‑236 (MQ=255)
                                                             ggTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGGTGCGGGGAACCGCAAACCCCCCACcttttc                                                                                                                                                                                                                      >  1:3425762/1‑201 (MQ=255)
                                                                        gTCAATCACCCCCCCGCCTATGTCGCCCCCCCCACAAAAGATATACATCCGGCTCAACATTGGGCCCGCCAGCACGGCGTGTCCCAAAAACAGAAGCGGAGAAATGAACGGGAAGGTAAATTCCGGGGGTTCGTTAACCCCCCCCAGCATGGCGGTAAAGGTGGCGGAAAACACCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCa                                                                                                                                                                     <  2:500802/239‑1 (MQ=255)
                                                                                                      cacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGTTTCGGGAATCAGACAACCCGCCCCCTTTATCCGCTTTTCCCGTGCGGGCGGGAAGGAAAACGCTTAAAAAAAGCCCCCCCCACCAAAGCATTTTTAGTTACCCTGCa                                                                                                                                    >  2:1174966/1‑239 (MQ=255)
                                                                                                        caccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGGGGGGATCAGCAAAACCCCCCCCTTTACCCGATTTTTCCGAGCGGCAGGGAAGTACAATGGCCAAAAAAATCCCCACGCCCCACAAAATTTTTAGTTTACCGGCAgg                                                                                                                                  >  2:1227261/1‑239 (MQ=255)
                                                                                                         accaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGc                                                                                                                                                                                                                                                    >  2:3079785/1‑132 (MQ=255)
                                                                                                         accaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGc                                                                                                                                                                                                                                                    <  1:3079785/132‑1 (MQ=255)
                                                                                                                   tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCCAACAGCAACGGCGAAATGAACAGGAAGGCAAATTCCAGCGGTTCTGTAATTCCCACCAGCATGGCGGTTAATGTTTCGGGGATCTGACAACCCGCCCCCTTTTACCCATTTTCCGGGGCGGGCGGGAAGGAACTCACTCAAAAAAAAACCCCCGCGCCCATAACTTTTTAGTTATCCGTGCTGGGCaaaagacg                                                                                                                        >  2:1054440/1‑235 (MQ=255)
                                                                                                                     cATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGAAACAACAAGACCGCCCCCTTTACCCCGTTTTTCCGTGCGGCCGGGAAGAACAACGCGAAAGAAAAACCCACGGCACCAAAAAATTTTTAGTTTACCGGCACGGCAAAAcaaccc                                                                                                                      >  2:1794272/1‑233 (MQ=255)
                                                                                                                             tCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGTTTCGGGGAACAGCAAAACCCCCCCCTTTAACCGCTTTTTCGGGGCGGCAGTGGAgga                                                                                                                                                                                           >  1:2519771/1‑164 (MQ=255)
                                                                                                                                aCATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGTTGGGGGGATTAGCCAGCCCGGCACtttttc                                                                                                                                                                                                                      >  1:2443745/1‑134 (MQ=255)
                                                                                                                                aCATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGGAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGGTGCGGGGAACAACAAAACCCCCCCCTTTTACCCGGTTTTCCGGGCGGGCGGGAAGGAAACACCCAAAAAAAAGCCCCCCGCCCCAAAAATTTTTTGGGTTCCCGCGAGGGCAAAAACGCCTTCCGGGgaccac                                                                                                         >  2:1757724/1‑239 (MQ=255)
                                                                                                                                      aGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGTTGCGGGGATCAACCAACCCCCCCCCTTTTCCCCCTTTTTCGGTGGGGCCGGGAAAGAAAACGCTAAAAAAAAACCCCCCGCCCCCAAAATTTTTTAAGTTTCGTGGCGGGCAAAAACGCCCTCCTGGgacaaccaattt                                                                                                   >  2:2846435/1‑233 (MQ=255)
                                                                                                                                                    gtgtctACCGCAAACACCCAGGGTAAAATAAAAGGGAAGGTAAATCCCGGGGTTTCGGTAATTCCCACCGCCAGGGCGGTTGGGGTGGGGGGAATCAGAAAACCCGCCACTATTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCAc                                                                                                                                                                    <  2:368282/159‑1 (MQ=255)
                                                                                                                                                          aCCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAGGGTGCGGGGAACAACAAACCCCCCCCCTTTTACCCCTTTTTCCGGGGGGGCGGGGAGAAAAACCCCAAAAAAAATCCCCCCGCACCCAAAAATTTTTGGGTTACGCGCAGGGGAAAAACCCCCTCCCGGGACAAAGCTTTCTCCGGCTCGgcgccccaaa                                                                                >  2:1133718/1‑232 (MQ=255)
                                                                                                                                                             gaaaaCAGCAAGGTGAAAGTGAAAGGGAGGGTAACTCCCGGGGGTGCGGAAACCCCCCCAGGCTGGGGGGTTAGGGGTGGGGGAATCACCACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGaa                                                                            <  2:1089777/237‑1 (MQ=255)
                                                                                                                                                              cAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGAATTAGCCAACCCGCCACTTTTACCCGATTTTCCCGTGCGGCAGTGAAAGAAATCGCTAAAGAAAAGCCCCCCGGCCCCaaa                                                                                                                                                           >  1:368282/1‑164 (MQ=255)
                                                                                                                                                                 aCAGAAACGGTAAAATGAACAGGAAGGAAAATCCCAGGGGTTCGGTAATCCCCACCAGCAGGGGGTTTAGGGTTGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTg                                                                                                                                                                                                        <  1:2331164/117‑1 (MQ=255)
                                                                                                                                                                   aaaaaCGGTAAAATGAACGGGAAGTTAATTCCCGGCGGTTCGGTATTCCCCACCAGCTGGGCGTTTAGGGTTGGGAAAAAAACAAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCg                                                                                               <  1:437563/211‑1 (MQ=255)
                                                                                                                                                                      aaCGGGGAAATGAACAGGAAGGAAATTCCCGGCGTTTCGTTAATCCCCACCAGCATGGCGGTTAGGGTGGGGGGAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCg                                                                                          <  1:2453067/216‑1 (MQ=255)
                                                                                                                                                                          ggaaaaaGAACAGAAAGTAAAATCCCAGGGGTCCGGTAATCCCCACCACCGGGGGGGTAAGGGTGGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtg                                                                <  1:1047547/235‑1 (MQ=255)
                                                                                                                                                                          gTGAAATGAACAGGAGAGTAAATGCCAGGGTTTCGGTAACCCCCACGGGCAGGGGGTTTAAGGTGGGGGAAATCAGCAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtg                                                                <  2:279067/238‑1 (MQ=255)
                                                                                                                                                                           tGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTTCGGGGATCAGCAAACCCCCCACCTTTTCCCCATTTTTCCGTGCGGGAGGGaa                                                                                                                                                                                              >  2:2331164/1‑117 (MQ=255)
                                                                                                                                                                             aaaTGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAAGTTGCGGGGATCAGCAAGCCCCCCCCTTTTTCCCCATTTTTCGGTGCGGCAGTGAAGTACATCGCTTAAGAAATGCCCACGGCACCAAAGAtttttta                                                                                                                                                 >  1:640393/1‑155 (MQ=255)
                                                                                                                                                                             aaaTGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAGGGTGGGGGGATCAg                                                                                                                                                                                                                                          >  2:2384981/1‑73 (MQ=255)
                                                                                                                                                                             aaaTGAACAGGAAGGAAAATCCCGGGGTTTGGTTAATCCCCCCCAGCAGGGCGTTTAGGGTGGGGGAAACCGGCAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtg                                                                <  1:2194497/235‑1 (MQ=255)
                                                                                                                                                                              aaTGAACAGGAAGGAAATTCCCGGGGGTTCGGTATTCCCCCCCAGCATGGGGGTTAGGGTGGGGGAAATCAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                           <  1:153194/239‑1 (MQ=255)
                                                                                                                                                                                 tGAAAAGGAAGGAAAATTCCGGCGTTTCGTTAATCCCCCCCAGCAGGGCGGTTAGGGTGGGGGAAATAAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                         <  1:1133718/239‑1 (MQ=255)
                                                                                                                                                                                  aaaagggAAGTTAAATCCCAGCGGTTCGTAAATTCCCCCCATGAGGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGccc                                                        <  1:2376252/234‑1 (MQ=255)
                                                                                                                                                                                   aaCAGGAAGTTAAATTCCGGCGGTTGGGTAATCCCCCCCAGCAGGGCGTTTAGGGTGGGGGAAACCACCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGAt                                                                                                      <  1:1279858/192‑1 (MQ=255)
                                                                                                                                                                                          aGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGa                                                                                                                                                                                                                                                <  1:1843742/56‑1 (MQ=255)
                                                                                                                                                                                           ggaaaaTCCCAGGGTTTCGTTAATCCCCCCCAGCAGGGGGGTTAAGGTGGGGGAAATCAGCAAGCCCGCCACTTTTACCCGAttt                                                                                                                                                                                                               <  1:3561767/82‑1 (MQ=255)
                                                                                                                                                                                            gaaaaTCCCGGCGTTTCGTTATTCCCCCCCAGCTTGGCGTTTAAGGTTGCGGAAAACAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                         <  1:1769382/226‑1 (MQ=255)
                                                                                                                                                                                             tAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGTTGCGGGGGACAGGAAGGCCGCCACTTTTTCCCCATTTTTCGGTGCGGCAGGGGAGGACATCCCCAAAGAAATGCCCACGGCACCCAAGAATTTTTAGTTACCGTGCAGGGCAAAACCGCCCTTCCGGGAACAAGATTTTAGCgggc                                                                                            >  2:1279858/1‑190 (MQ=255)
                                                                                                                                                                                             aaaaCTCCGGCGGTTGAGTAATCCCCGCAAGCGTGGGGTTTGAGTTGGGGGAAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                         <  2:1975488/226‑1 (MQ=255)
                                                                                                                                                                                                          ttCGGTAATCCCCACCAGCGTGGCGGTTAAGGTGGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAAct                                                                          <  1:566194/197‑1 (MQ=255)
                                                                                                                                                                                                               aactcccccccAGCATGGCGGTTAAGGTTGCGGAAATCAGCACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAg                                                                                                                                                <  2:2023933/120‑1 (MQ=255)
                                                                                                                                                                                                                 atccccACCAGCATGGCGGTTAAGGTGGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGc                         <  1:1174966/236‑1 (MQ=255)
                                                                                                                                                                                                                  tcccccccGGCAGGGCGTTTAAGGTGGGGGGAACCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAAt                                                                                                                                                                           <  1:1717048/95‑1 (MQ=255)
                                                                                                                                                                                                                      caccaGCGGGGCGGTAAATGTGGCGGAAATCAGCAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                    <  1:1237659/239‑1 (MQ=255)
                                                                                                                                                                                                                      caccaGCAGGGGGGTTAAGGGGGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                    <  1:884707/239‑1 (MQ=255)
                                                                                                                                                                                                                       cccaGCTGGGCGGTTAAGTTTGGGGAAAACAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                           <  1:2738076/198‑1 (MQ=255)
                                                                                                                                                                                                                       cccaGCAGGGGGTTAAGGGTTGCGGAAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCt                                                                                     <  1:138546/172‑1 (MQ=255)
                                                                                                                                                                                                                        ccaGCTGGGCGGTTAGGTTGGCGGGAATAAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                  <  1:3468314/239‑1 (MQ=255)
                                                                                                                                                                                                                        ccaGAAGGGCGGTGAGGGTGGCGGAAACCAACAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                  <  1:2199944/239‑1 (MQ=255)
                                                                                                                                                                                                                         caGCATGGCGGTTAAGGGTGCGGGGATCAGCAAACCCCCCCCTTTTTCCCGCTTTTCCCGTGCGGCCGGGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTTAGGTTCCGCGCAGGGCAAAACCCCCTTCCGGGGACAACGATTTTTGCCGCTCCGCCGTCCACCCGCACTCTTGCTGAAGACGCGCCCAGTACCTATCACTGCCCCCCTCTACAGCTccccggcaca                 >  2:2199944/1‑232 (MQ=255)
                                                                                                                                                                                                                            cATGGCGGTTAAGGTTGCGGGGATCAGCAAGCCCCCCCCTTTTTCCCGATTTTCCGGGGCGGGAGGGGAGGACa                                                                                                                                                                                        >  1:3029603/1‑74 (MQ=255)
                                                                                                                                                                                                                            cATGGCGGTTAAGGTGGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCgg                                                                                                                                                                                                          <  1:448399/57‑1 (MQ=255)
                                                                                                                                                                                                                                 cGGTTTAGGGTGCGGGAATCAGCAAGCCCGCCCCTTTTTCCCCATTTTTCCGGGCGGCAGTGGAGGACATCGCTTAAGAAATGCCCCCGGCACCCAAAATTTTTTAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTTAGCGGCTCCGCGGTCCAAACGAAc                                                                           >  2:138546/1‑173 (MQ=255)
                                                                                                                                                                                                                                  ggtttAAGGTTGCGGAAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa          <  1:2142843/235‑1 (MQ=255)
                                                                                                                                                                                                                                          gagtggCGGGAAACCAGCAACACCGACACTTTTACCCGATTTTCCGGTGCGACAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa          <  2:2922638/228‑1 (MQ=255)
                                                                                                                                                                                                                                            gtggCGGAAACCAGCACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGt  <  1:3410793/235‑1 (MQ=255)
                                                                                                                                                                                                                                             ttGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCa                                                                                                                                                             >  1:2078901/1‑86 (MQ=255)
                                                                                                                                                                                                                                             tggCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGt  <  1:2089341/235‑1 (MQ=255)
                                                                                                                                                                                                                                                cGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATgccgcc                                     >  2:2274838/1‑199 (MQ=255)
                                                                                                                                                                                                                                                cGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATgccgcc                                     <  1:2274838/199‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  
AATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGT  >  NC_000913/3862428‑3862897

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