Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,6640CG56.1% ‑3.8 / 32.9 42pseudogene (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 (8/10);  new base G (7/16);  total (16/26)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.11e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATG  >  NC_000913/3862426‑3862902
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
ccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGc                                                                                                                                                                                                                                                      <  1:1490811/239‑1 (MQ=255)
   aTGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCggg                                                                                                                                                                                                                                                  <  1:372146/239‑1 (MQ=255)
         aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGTTGGGGGGATCa                                                                                                                                                                                                                                            >  1:291438/1‑238 (MQ=255)
          tttGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATCACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAAGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTGGGGGGATCaat                                                                                                                                                                                                                                          >  1:934835/1‑237 (MQ=255)
                        atatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGTGGCGGGAAACAGCAAACCCGCCCCttt                                                                                                                                                                                                                           >  1:223957/1‑239 (MQ=255)
                                   gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGGTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTGCGGGGGTCAGCAAGCCCGCCCCTTTTTTCCGActtt                                                                                                                                                                                                                >  2:2228996/1‑239 (MQ=255)
                                            cATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAAGGTTGGGGGAAAAGGAAAGCCCGCCCCTTTTCCCCGATTTTCCgggtgg                                                                                                                                                                                                        >  1:1905859/1‑235 (MQ=255)
                                            cATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTTCGGGGAACAGCAAACCCGCCACCTTTTTCCCCTTTTTCGCTGCg                                                                                                                                                                                                        >  1:1328351/1‑238 (MQ=255)
                                             actggggtactgTTTTGCGGAAAAAACCTAGAAAACCCCCCCCCCGCCCTTGTCCCCCCCCACCCAAAAAAAACCATTTCGGGTGAACATTGGGCCCGCCGGCACGGGGGTCCCCGAAAACAACAACGGAAAAAAAACAAGGAAGGTAAATCCCGGGGGTTCGGTAACCCCCCCCAGCAGGGCGTTTAAGGTGGGGGAAAACAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCgg                                                                                                                                                                                                       <  2:316389/228‑1 (MQ=255)
                                              gcgacaTCCATTTTTGCGAAAAACCCTAGACATCCACACCCCCGCCTATGTCCCCCCCCCCACCAAAGAGATAATTTACGGGCGACAGTGGCGCCCCCAGCCAGTCTTGGACCGAAAACGGAAACGTAAAAAAGAGGAGGAAGAAAAGTCCCGGCGTGAGGTAAATCCCCACAAGGTGGGGGGTTATGGGTGGGGGAGAAAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCa                                                                                                                                                                                                     <  2:403557/234‑1 (MQ=255)
                                                  atcccAGTTTTGCGGAAAACCATAGTCACTCACACCCCCCCCTATTTCCCCCCCCACACAAAAGGAAACATTGGCGGACGACAGTGAGGCCGCCAGCGCGGCGTGTCCCGCAACCAGCAACGGAAAAATGAAAAGAAAGGAAAATCCCAGCGGTTCGGTACTCCCCACCAGAATGGGGTTTAGGGTGGGGGAAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTg                                                                                                                                                                                                  <  2:1108314/235‑1 (MQ=255)
                                                   gatcAGTTTTGCGGAAAAAACTAAGTCACCCAGCCCCCCCCCAATGTCCCCCCCCCCCAAAAAGAAATCAATTACGGTCACAATGGAGGCCGCCAGCCGGGGGTGTACCGAAAACAGAAAGGGTAAAATGAACAGAAGGATAAATCCCGGCGGTTCGGTAATCCCCACCGGCGGGGGGTTTAGGGTGGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCa                                                                                                                                                                                                     <  2:1476639/232‑1 (MQ=255)
                                                               ggTAAAACCTAGTAAATCAGCCCTCCCCCCATTTGGCCCACCACCCCAAAGGGTACCATTGCGGTCGCCATTGAGGCCGCCAGCACGGCGTTCACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAGGTTTGCGGGAATCAGCAAGCCTGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTa                                                                                                                                                                                   <  2:1654495/238‑1 (MQ=255)
                                                                             tCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGGTGCGGGGATCAGCAAAACCCCCCCCTTTTCCCCACTTTTCCGGTCGGCATGGAGGTACATCTCTCAAAAAAATTccc                                                                                                                                                                        >  2:112242/1‑238 (MQ=255)
                                                                                         ccGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAAGGGTGGGGGAAAAGCAAACCCCCCACCTTTTCCCCAtttttcc                                                                                                                                                                                                             >  1:1512463/1‑187 (MQ=255)
                                                                                           ccccttgtgccccACAACACAAAAGAGATCAATGGCGGTCGACATGGAGGCCCCCAGCACGGGGGTTACCGAAAACAGCAACGGAAAAAGAAACAGGAAGGAAAATTCGAGCGTTTCGGTAACCCCCCCCAACAGGGGGTTTAGGGTTGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttt                                                                                                                                                        <  2:456380/229‑1 (MQ=255)
                                                                                                       caccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCggg                                                                                                                                                                                                                                                  <  1:427734/141‑1 (MQ=255)
                                                                                                          cacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGa                                                                                                                                                                                                                                                <  1:2117144/139‑1 (MQ=255)
                                                                                                          cacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAAAGAACAGGAAGGTAAATTCCAGCGGTTCGGGAATTCCCACCAGCATGGGGGTATAAGTAATGGGCAGCCGCACACCCGCCCCTTTTTCTCGAAATTCCGTGGCGGCGGGGAAGAAGAACGCTAAAAAAAATACCACGGCGACACAAATATTTGAATTAcacagac                                                                                                                                         >  2:379407/1‑233 (MQ=255)
                                                                                                                 aaGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGGTTCGGGGATCAGccaa                                                                                                                                                                                                                                       >  2:427734/1‑138 (MQ=255)
                                                                                                                    agaTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGGTTGGGGGAATAAGCAACCCCCCCCCTTTTACCCGATTTTCCGGGGGGGCGGGGAGGAAAATGCCTAAAAAAAATCCCCACACACCAAAAAATTTTTAGTTACCGAGCAGGGCcaaa                                                                                                                                >  2:2123212/1‑238 (MQ=255)
                                                                                                                    agaTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAGGGTGCGGGGAACAGCaaaa                                                                                                                                                                                                                                     >  2:2117144/1‑137 (MQ=255)
                                                                                                                                     gACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGGGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGTTGGGGGGATCCCCCACACCGCCCCCTTTTCCCCATTTTCCCGGGGGGGCGGGGAGAACAACCCTAAAAAAAATTCCCCGCGACCAAAAATGTTTTAGTTGTCGTGCACGGGGAAAACGCCCTTCCGCGacca                                                                                                                >  2:404150/1‑239 (MQ=255)
                                                                                                                                            aGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGGTTCGGGGAAAAGCCAGCCCCCCCCCTTTTTCCCCTTTTTCGGGGGGGGAGGGAAGAAACACGCCAAAAAAAAACCCCCCGGCCCCAACAAATTTTAGTTTCCCCGCAGGGGAAAAACCCCCTCCCGGGaaaaaca                                                                                                             >  2:2054666/1‑234 (MQ=255)
                                                                                                                                                                 ccGAAACCAGCACGGGTAAAATAAACAGAAAGTAAAATTCCGGCGTTTCGTTAATCCCCACCAGCTTGGCGTTTAAGTTTGCGGAAACCACCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTg                                                                                   <  1:365922/239‑1 (MQ=255)
                                                                                                                                                                               gttaaaagcaaCAGGAGGGAGAATTCCGGCGGTTGAGATGTTCCCACCGGCGTGGGTATTAGGGTGGGGGGAATCAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGGTCGGCGCTCAAACTGAACTCTTGCAg                                                                        <  2:2202672/227‑1 (MQ=255)
                                                                                                                                                                                  gAAATGAAGAGGAAGGTAAATCCCTGGGTTGCGTGAACTCCCACCAGGGTGGCGGTTAGGGTGGGGGAAACCAGAAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                   <  1:442161/238‑1 (MQ=255)
                                                                                                                                                                                  aaattGAACAGGAAGGAAACTCCCGGCGGTTCGTTAATCCCCACCAGCGTGGCGGTTAGGGTGGCGGAAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                   <  1:305533/234‑1 (MQ=255)
                                                                                                                                                                                   aaaTAAAAAGGAAGGAAAATCCCGGCGTTTCGGTACTCCCCCCCAGAAGGGCGTTTAGGTTGGGGGAAACCACCACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTg                                                                                                                                                     <  1:154012/156‑1 (MQ=255)
                                                                                                                                                                                        acagggaaggaaACTTCCGGGGGGTCGGAATCTCCCCCAGGATGGGAGGTAAGGGTGGGGTAAATCAGCAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                              <  2:1354589/233‑1 (MQ=255)
                                                                                                                                                                                         aaaggaggGTAAATGCCGGGGGTTGGATAATCCCCCCCAGCATGGGGTTAAAGGTTGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                                           <  1:627355/237‑1 (MQ=255)
                                                                                                                                                                                             gAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGGGGGAATCAGCAAGCCCGCCACTTTTTCCCGATTTTCCCGGGCGGCAGTGAAGGACATCGCTAAAGAAATGCCCCCGGCACCCAAGATTTTTTAGTTACCGTg                                                                                                                                           >  2:154012/1‑156 (MQ=255)
                                                                                                                                                                                             gAAGGAAAATCCCGGCGGTTGGGTAATCCCCACCAGCAGGGGGGTTAGGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACAt                                                                                                                                                                                         <  1:1489754/111‑1 (MQ=255)
                                                                                                                                                                                              aaGGTAATTCCCAGCGGTTCGGAATTCCCCACCAGCAGGGCGGTTAAGGTGGGGGGAATAAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                   <  1:1547480/226‑1 (MQ=255)
                                                                                                                                                                                                  aaaaTTCCGGTGGTTCGGTAACCCCCCCCAGCGGGGGGGTAAGGGTGGGGGAAACCAGAAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGt                                                          <  1:1615177/230‑1 (MQ=255)
                                                                                                                                                                                                     ttcccggCGGTTGGGTAATTCCCCCCAGTGTGGGGGTTAGGGTGGGGGAAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATg                                               <  1:126137/233‑1 (MQ=255)
                                                                                                                                                                                                     attcgggCGGTTCGGTATTTCCCACAAGTGGGGCGGTTAAGGTGGGGGAAATCACCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATg                                               <  1:139042/233‑1 (MQ=255)
                                                                                                                                                                                                      ccccgggggTTGGGTAATCCCCCCCGGAGGGGGGGTTAGGGTGGGGGGAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATgc                                              <  1:2054666/232‑1 (MQ=255)
                                                                                                                                                                                                       tCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCAGCAAACCCGCCCCTTTTACCCGATTTTCCGGTGCGGCAGGGAAGGACAACCCTAAAGaa                                                                                                                                                                              >  2:1489754/1‑111 (MQ=255)
                                                                                                                                                                                                                tcagTAATTCCCACCAGCGGGGCGTTTAAGGGTGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCaa                                     <  1:2039007/235‑1 (MQ=255)
                                                                                                                                                                                                                            ccctgccggggCGTTTGAGTTGGCGGAAATCAGCGACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                         <  2:658647/230‑1 (MQ=255)
                                                                                                                                                                                                                            ccccaGCTGGGGGGTTAAGGTTGGGGAAAACAGCACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                         <  1:2123212/237‑1 (MQ=255)
                                                                                                                                                                                                                            ccacaacgtGGCGGTTGGGGTTGGGGAACACAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                         <  1:1406220/231‑1 (MQ=255)
                                                                                                                                                                                                                                                     cGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGcc                                                                                                                                                                         >  1:808480/1‑73 (MQ=255)
                                                                                                                                                                                                                                                     cGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATg  >  1:1146295/1‑239 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
CCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATG  >  NC_000913/3862426‑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: