Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,6670GA31.8% 57.5 / 10.5 44pseudogene (937/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 (21/9);  new base A (0/14);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 8.12e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGA  >  NC_000913/3862430‑3862903
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGa                                                                                                                                                                                                                                                     <  2:184913/239‑1 (MQ=255)
                    atatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGTGTCGGGGATAAGCAAGCCCGccc                                                                                                                                                                                                                                   >  1:1447203/1‑234 (MQ=255)
                     tatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTGGGGGGAACCGAAAAGCCGCCACcttt                                                                                                                                                                                                                            >  1:904526/1‑239 (MQ=255)
                           cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCTGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGTTGGGGGGAACAACAAACCCGCCCCCTTTTccc                                                                                                                                                                                                                        >  1:561426/1‑237 (MQ=255)
                           cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGTGGGGATCAACAAACCCGCCCCCTTTTAccca                                                                                                                                                                                                                       >  1:114555/1‑238 (MQ=255)
                              ggTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGGGGGGGGGAAAAGCAAACCCGCCCCTTTTACCCgcgtt                                                                                                                                                                                                                   >  1:579312/1‑235 (MQ=255)
                                 tGCTGAACACCGGGAGCCAGTTTGGCGAAAAACCCTATACAATCCGCCCCCCGCCCTTGTCCCCCACACCACAAAAGAAATACATTACGGCCGACAGTGAGGCCGCCAGCACGGCGTGTACCACAAACAGCAAGGGAAAAAAAAACAGGAAGGTAAATCCCGGCGGTTGGGTACTCCCCACCCGCAGGGGGTTTAGGGGGGGGGGAATAAGCAAGCCCGCCACTTTTACCCGATTTTc                                                                                                                                                                                                                 <  2:27641/238‑1 (MQ=255)
                                                     gttgggggTAAAACCTACGCAATCCGCCCCCCCGCCAATGTGCCCCACCACACCAAAGAAATACATTCGGGGCGACATGGGGGCCGCCAGCAGGGCGTGCACCGAAAACAGCACGGGAAAAAGGAAAGGAAAGTTAAATTCCGGGGGTTGGGTAATCCCCACCAGCAGGGGGTTTAGGGTGGGGGAAATAAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACa                                                                                                                                                                                            <  2:1073521/233‑1 (MQ=255)
                                                      tttGGGGAAAAACCTAGCCATCCACCCCCCCCCCATGTTCGCCCCCCCCCCAAAAGAAAAACATTACGGCGACAATGGGGCCCGCAGGCCCGGGGTGTACCGAAACCAGCAACGTGAAAATGAACAGGAAGGAAAATCCCAGGGGTTCGGTAATCCCCCCCAGCGTGGGGGTTAGGGTGGGGGAAATAAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACAt                                                                                                                                                                                           <  2:1467764/238‑1 (MQ=255)
                                                      tttGCGGATAAACCCAAGCCAATCAGCCCCCTGCCCTTGTTGCCCCCCACACAAAAGAGATACATGCGGGACAGCATGGGGCCCCCCAGCAGGGGGGGTACCACAAACAGCAACGGTAAAATGAGCGGGAAGGTAATTTCCAGGGTTTCGTTAATCCCCACCAAGAGGGCGGTAAGGGTGGGGGAAATCAGCAGCCCCGCCACTTTTACCCGATTTTCCGGTGCGGCa                                                                                                                                                                                                       <  2:1398875/228‑1 (MQ=255)
                                                                            tCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTGCGGGGAATAGCAAACCCCCCACTTTTTACCCGTTTTTCCGTGCGGGCGTGGAGGACATCGCTAAAAAAAACCCCACCGc                                                                                                                                                                     >  1:1490662/1‑238 (MQ=255)
                                                                                                     cacacaccaAAGACTTAACTTAGGGTCGACATGGGGGCCGCCAGCGGGGCGTGACCCAAAAACAGAAAGGGAAAAAGAAACAGAAAGGAAAATCCCGGCGGTTCGTTAACCCCCCCCAGCGTGGGGGTAAGGGGTGCGGGAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTg                                                                                                                                             <  2:904526/237‑1 (MQ=255)
                                                                                                      ccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTGTGGGGAACAACAAAACCCCCCCCTTTTCCCCATTTTTCCGGGCGGCAGGGAAAATAATAGCCTAAGAAAAGACCCCCGGACCCACAGATTTTTTAGTTACCggt                                                                                                                                             >  2:1162498/1‑237 (MQ=255)
                                                                                                         caccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTATGGTTGGGGGagcaa                                                                                                                                                                                                                                                 >  2:288618/1‑138 (MQ=255)
                                                                                                          ccaaaaGATAAACATTCCGGTCGACTTGGGGCCCGCCAGCGGGGCGTGACCCAAAAACAGCAAGGGTAAAAAAAAAGGGAAGTAAAATCCCGGCGGTTGGGTAATCCCCCCCAGCAGGGCGGTTAGGGTGGGGGAAATCAAAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAc                                                                                                                                                                                              <  2:1276587/184‑1 (MQ=255)
                                                                                                                   aTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGAAAAAACAAAACCCCCCCCTTTTACCCATTTTTCGGGGGGGCAGTGAAATACATCGCTCAAGAAAAACCCACGGCCCCAAAAATTTTTTAGTTAACCTGCAGGGGAAAACCg                                                                                                                              >  2:290015/1‑238 (MQ=255)
                                                                                                                    tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGTTTCGGGGATCAGCAAACCCGCCACTTTTACCCCATTTTTCCGTGCGGCAGGGGAGGAAATCGCTtaaa                                                                                                                                                                                   >  1:1276587/1‑186 (MQ=255)
                                                                                                                         tACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGTTTCGGGGATAAGCAAACCCCCCCCCTTTTCCCGATTTTCCGGGGCGGCGGGGAAGGACa                                                                                                                                                                                            >  1:635479/1‑173 (MQ=255)
                                                                                                                              tCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTGCGGGGAACAGCAAGCCCGCCCCCTTTAACCCAATTTTCCGGGCGGCCGGGAAgga                                                                                                                                                                                               >  1:347580/1‑164 (MQ=255)
                                                                                                                                 aCATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGAATAGCAAGCCCGCCACCTTTTTCCGGATTTTcg                                                                                                                                                                                                                >  1:213663/1‑146 (MQ=255)
                                                                                                                                     tGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGGTGTGGGGAACAAAAAACCCCCCACCTTTTCCCCCATTTCCGGGGGGGCCGTGAAAAACAAACCCAAAAAAAAACCCCCCGGCCCCAAAAATTTTTGGTTTACATGCACGGCCAAAACCCCCTTCCGGGaacaa                                                                                                                >  2:1561387/1‑235 (MQ=255)
                                                                                                                                        ggCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGGTTGGGGGAAAAGGAAAACCGCCACCTTTACCCCATTTTTCGGGGCGGGAGGGGAGGAAATTGCCAaaaaaaaacccc                                                                                                                                                                          >  1:147953/1‑170 (MQ=255)
                                                                                                                                                       cgtcgaCCGAAAACAGCAGAGGAAAAATAAACAAGAGGGAAAATCCAGGGGTTTGGGTAATCCCCCCCAGCATGGGGTTTAGGGTGGCGGGAATCAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCt                                                                                           <  2:1490662/233‑1 (MQ=255)
                                                                                                                                                                   aCAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAAGTTGCGGGGATCAGCAAACCCCCCCCCTTTACCCCATTTTCCGGGGCGGCAGTGGAGGAAATCGCCAAAGAAATGCCCACGGCACCAAAAATTTTTTTGTTACCGTg                                                                                                                                             >  1:1089524/1‑178 (MQ=255)
                                                                                                                                                                            gTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTTCGGGAATCAGCAAGACCGCCCCCTTTTACCCATTTTTCGGGGCGGCCGGGAAGTACATCGCTAAAAAAAAACCCCCCGCCACAAAAATTTTTTAGTTTCCGCGCCGGGGAAAAACCCCCTTCCGGGAAAAACATTTTAACGGCCTCGCGCTCCAACAGAAcaccttccgaa                                                                       >  2:1066940/1‑227 (MQ=255)
                                                                                                                                                                                  tGAACAGGAAGATAAATCCCGGCGGTTCGGTAATCCCCACCAGCATGGGGGTTAGGGTTGGGGAAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                               <  1:224484/239‑1 (MQ=255)
                                                                                                                                                                                  agAACGGGAAGGTAAATCCCGGGGGTCGGGTAATCCCCACCAGTATGGCGGTTAAGGTGGGGGAAAGCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                               <  1:774962/238‑1 (MQ=255)
                                                                                                                                                                                    aaaggaaagaaaaaTCCCGGGGGTTCGGTATTCCCCACCAGCAGGGGGGTTAAGGTTGGGGAAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                               <  1:14256/226‑1 (MQ=255)
                                                                                                                                                                                         gaattaaaaTTGCGGCGGTGTGGTAATCCCCACCAGCGTGGGGTTAAGGGTTGGGGAAACCAGAAGCCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTa                                                          <  1:159368/231‑1 (MQ=255)
                                                                                                                                                                                          aaGTTAAATTCCAGCGGTTCGTTAATCCCCCCCAGCTTGGCGGTTAAGGTGGCGGAAATCACAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                                            <  1:290015/234‑1 (MQ=255)
                                                                                                                                                                                                 ttcccAGCGTTTCGGTTATCCCCCCCAGCAGGGCGTTTAGGGTTGGGGGAATCAGCAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATg                                                <  1:1162498/236‑1 (MQ=255)
                                                                                                                                                                                                    tccggcttttaggTAATGCCCACCAGCGTGGTGTGTATGGTTGCGGAAATCAACAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATg                                                <  1:988040/226‑1 (MQ=255)
                                                                                                                                                                                                       gCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGGGGGGATCAGCAAGCCCCCCCCCTTTACCCGATTTTTCGGGGCGGGAGGGAAGTAAATCGCCAAAAAAAAGCCCACGGCCCCAAAGAATTTTTAGTTACCGGGCAGGGGAAAACCCCCTTCCGGGAACAACGAATTCAGCGGCCCGGCGCTCCCAACGAAATCTTTTAGAAGCTGCGCCCCGTAAAACCGAATGCCGCCt                                          >  2:1040929/1‑239 (MQ=255)
                                                                                                                                                                                                               ggTAATCCCCCCCAGAAGGGCGGTTAAGGTTGCGGAAATAAGCAAGCCCGCCACTTTTAcc                                                                                                                                                                                                                         <  1:405493/61‑1 (MQ=255)
                                                                                                                                                                                                                    ccccccccACCATGGCGGTTAGGGTTGGGGAAAACAGCAGCCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGcc                              <  1:1411799/237‑1 (MQ=255)
                                                                                                                                                                                                                        cccgaggaTGGGGGTTAGGGTGGGGGGAATCAGAAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                          <  1:333249/232‑1 (MQ=255)
                                                                                                                                                                                                                           ccagagggggggTTGAGGTTGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                        <  1:1304340/229‑1 (MQ=255)
                                                                                                                                                                                                                                   cGGTTAGGTTTGGGGAAACCAGAAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTa                                                                                                                                                  <  1:1612183/110‑1 (MQ=255)
                                                                                                                                                                                                                                      gTTAGGGTGGGGGGATCCAGCATGACCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGaa                                                                                  <  2:1213647/170‑1 (MQ=255)
                                                                                                                                                                                                                                        aaaGGTGGCGGGAATCAGAAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTc           <  1:467194/238‑1 (MQ=255)
                                                                                                                                                                                                                                         gaggggggggggAACCAGCACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTc           <  1:1040929/230‑1 (MQ=255)
                                                                                                                                                                                                                                                gCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTa                                                                                                                                                                                               >  2:1062317/1‑56 (MQ=255)
                                                                                                                                                                                                                                                gCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTa                                                                                                                                                                                               <  1:1062317/56‑1 (MQ=255)
                                                                                                                                                                                                                                                 cGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGa                                                                        >  1:1213647/1‑170 (MQ=255)
                                                                                                                                                                                                                                                 cGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATg   >  1:1325144/1‑239 (MQ=255)
                                                                                                                                                                                                                                                  gggAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACa                                                                                                                                                                                            >  1:305617/1‑56 (MQ=255)
                                                                                                                                                                                                                                                  gggAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttt                                                                                                                                                          >  1:1064023/1‑89 (MQ=255)
                                                                                                                                                                                                                                                  gggAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGa  >  1:1534457/1‑239 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
TGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGA  >  NC_000913/3862430‑3862903

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: