Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,6680AG30.5% 80.1 / 12.8 59pseudogene (936/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 A (8/33);  new base G (18/0);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 2.41e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATG  >  NC_000913/3862435‑3862902
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTGGGGGGaaca                                                                                                                                                                                                                                                  >  1:1975488/1‑236 (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)
                                                                                                            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)
                                                                                                                                                                                                                                     gtgggggaaaTAAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGt       <  1:1496621/230‑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)
                                                                                                                                                                                                                                          ggaaaTCAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATg  <  1:2799748/235‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
ATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATG  >  NC_000913/3862435‑3862902

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