Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,855,9980TC36.2% 68.9 / 17.1 48pseudogene (966/1617 nt)BW25113_RS19145PTS sugar transporter subunit IIB
Reads supporting (aligned to +/- strand):  ref base T (22/8);  new base C (0/17);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 4.04e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAC  >  NZ_CP009273/3855751‑3856239
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                            
ggAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGCGAAATGAACAGGAAGGCAAATTCCAGCGGTTCGGTAATTcc                                                                                                                                                                                                                                                            >  1:675314/1‑250 (MQ=255)
                 agcagGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGtt                                                                                                                                                                                                                                          >  1:18016/1‑249 (MQ=255)
                     gggTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGt                                                                                                                                                                                                                                     <  2:437594/250‑1 (MQ=255)
                      ggTAAAGCGTAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGtt                                                                                                                                                                                                                                    <  2:103037/250‑1 (MQ=255)
                       gTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTg                                                                                                                                                                                                                                   <  1:458379/250‑1 (MQ=255)
                        tAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGGTGc                                                                                                                                                                                                                                  >  1:674220/1‑250 (MQ=255)
                          aaGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCgg                                                                                                                                                                                                                                <  1:365545/250‑1 (MQ=255)
                          aaGCACAACACAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCgg                                                                                                                                                                                                                                <  1:573096/250‑1 (MQ=255)
                               caACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATc                                                                                                                                                                                                                           <  1:130964/250‑1 (MQ=255)
                                 aCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGCAAGGGTGCgggcaaccg                                                                                                                                                                                                                         >  2:246964/1‑244 (MQ=255)
                                      aTGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGGGGGAATAAGCaaac                                                                                                                                                                                                                  >  1:471786/1‑248 (MQ=255)
                                         gCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGccc                                                                                                                                                                                                                <  1:693552/249‑1 (MQ=255)
                                              ttGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGGAAATTCCAGCGGTTCGGTAATTCCCAACAGCATGGCGGGTTAGATTGGGGGGAGCAGCAAAGCCCCCAcc                                                                                                                                                                                                         >  2:690830/1‑249 (MQ=255)
                                                       catcatATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTGAGGGAAGAAGCAAAcccccccccctttacccc                                                                                                                                                                                                >  2:70413/1‑238 (MQ=255)
                                                          catATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGTGGCGGGGAATAGCAAACCCCCCACCTTTTcccccttt                                                                                                                                                                                             >  1:336299/1‑245 (MQ=255)
                                                          catATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTACCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTCAGGGTGCCGGGATAGTAAAACCACGCACTTTTACCCCAttt                                                                                                                                                                                             >  2:502335/1‑250 (MQ=255)
                                                                          cTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGGTGCGGGGAACAGCAAGCCCCCCACTTTTACCCCGTTTTccagagagggacgg                                                                                                                                                                              >  1:415844/1‑237 (MQ=255)
                                                                          cTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTAGGTAATTCCCACCAGCATGGCGGGTAAAGGTGGGGGGAACCACAAGGCCCGCACATTTACCCCATTTTCCGGTGCGGccggc                                                                                                                                                                             >  2:296048/1‑245 (MQ=255)
                                                                                                gCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCAGCAAACCCCCCACTTTTACCCGATTTTTCCGTGCGGCAGCGAAGTAAATCGCTTAAGAAAAGcc                                                                                                                                                     >  1:660052/1‑250 (MQ=255)
                                                                                                 cGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAACCAGCAAGCCCGCCCCTTTTTCCCGATTTTTCGGTGCGGGCGGGAAGGACATCGCCACAAAAAAGccc                                                                                                                                                    >  2:627285/1‑250 (MQ=255)
                                                                                                       aaccgactcatccAGCCCCCCCCCCGTGTTCCCCCCCCCACTAAAGAGAAACATTACGGTGAACACTGAGACCGACAGCCGGGCGTGTCCCGAAACCCGAAAGAGAAAAATGAAAGAGAGGGTAAATTCGGGGGGTTCGGTAACCCCCCCCAGCGTGGGGGTTAAGTGAGCGGAAATCATCACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCa                                                                                                                                              <  2:465621/238‑1 (MQ=255)
                                                                                                         ccTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTTGGGGGAACAGCCAACCCCCCCCCTTTTCCCCGTTTTTCCGTGCGGCCGGGAAATAAATCCCCTAAAAAAATCCCCCCgcacc                                                                                                                                             >  1:505912/1‑246 (MQ=255)
                                                                                                               aatcccgccgtCCGCACTTGTGCCCCACCAAACAAAAGAGATAATTTGCGGTCGACTGTGGGGCCGCCAGCGGGGGTTTCACCAAAACCAGAATGGTGAAAATGAAGAGGAGGGAAATTCCCGGCGTTGAGGTAATACCCACAAGCTTGGGTGTTAGGGTGGGGGGAAAAAGAAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAt                                                                                                                                     <  2:667943/240‑1 (MQ=255)
                                                                                                                aaacacaCCCCCGCCGTTGTTCCCCACCCCAAAAAAGAAAAACATGCGGGTGAAAAGGGAGGCCGCCAGCAGGGCGTGTCCCGAAAACAGAAAGGGAAAAAAAAACAGGAAGGAAAATCCCAGCGGTTCGGAAACCCCCACCAGCAGGGGGGTTAGGGTGGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAtt                                                                                                                                    <  2:158080/244‑1 (MQ=255)
                                                                                                                aaTCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGGAACAGCAAACCCCCCCCCTTTTCCCCATTTTTCCGTGCGGCAGGGAAGTACAACCCTAAAAAAATGGCACCCGCACCACAAAtt                                                                                                                                    >  2:568095/1‑250 (MQ=255)
                                                                                                                    aGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAAGGTGCGGGGAACAGCCAGCCCCCCACTTTTTCCCCATTTTTCGGTGCGGGAGTGGAAGACATCGCTAAAAAAAAGCCCCCCGCACCAAAAATTTTTg                                                                                                                                >  1:602319/1‑250 (MQ=255)
                                                                                                                        cTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGTTGCGGGAATCAGCAAGCCCGCCCCCTTTTCCCGGTTTTTCGGGGCGGCCGGGGAGGACAACGACTAAGAAATTCCCACCGCACCCAAAGAtttttaggt                                                                                                                             >  1:226887/1‑246 (MQ=255)
                                                                                                                           cGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCAGCAAACCCCCCACCTTTTTCCGATTTTCCCGTGCGGCAGTGAAGTAAAACGCCAAAAAAAAACCCACGGGACCAAAGATTTTTTAGGTAcc                                                                                                                         >  1:78446/1‑250 (MQ=255)
                                                                                                                            gcccttttcccccACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGAAACGGTAAAATAAACAGGAAGGTAAATTCCAGCGTTTCGGTAATCCCCACCAGCAGGGCGGTTAAGGTGGCGGGAACCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCg                                                                                                                        <  1:544864/240‑1 (MQ=255)
                                                                                                                                atgtcccccACCACAACAAAGAGATACATTACGGTCGACAGGGAGCGCGCCAGCGCGGTGTATACCGAAAACAGCAGCGGAAAAATGAACAGGAAGGTAAATCCCAGCGGTTCGGTAACCCCCACCAGCGGGGGGGTTAGGGTTGGGGAAATCAGCACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCa                                                                                                                    <  2:120611/244‑1 (MQ=255)
                                                                                                                                      cccACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAAGTTGCGGGGAACACCAAAACCCCCCCCTTTTCCCCCTTTTTCCGTGCGGCCGGGGAGAAAATCCCCTAAAGAAACCCCCCCGCCCCAACAAATTTTTGAGTTACCGtccgggcc                                                                                                                >  2:101972/1‑242 (MQ=255)
                                                                                                                                            cccaaaaaGAAAAACATAGCGGTGGAGAGTGAGCCCGCCGGCCGGGGGGGCACCGAAAAAAGCAGCGGAGAAAGGAACAGGAAGATAAATCCCAGGGGTGGTGAAATCCCCACCAGCGTGGGGGTGAGGGTGGGGGAAAGCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGc                                                                                                       <  2:278578/245‑1 (MQ=255)
                                                                                                                                                caAAGAAATACATCCGGGCCGCAATGGAGGCCCCCCCCGCGGGGTGTCCCGAAAACAGAACGGGTAAAATGAAAAGGAAGATAATTCCCGGGGGTCCGGTATCCCCCCCCTGCTGGGGGTTTAGGGTGGCGGAAATCACCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCGGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTc                                                                                                   <  2:71866/250‑1 (MQ=255)
                                                                                                                                                    gagaTACATTCCGGTCAACATGGAGGCCCCCAGCAGGGCGGGTACCGAAAACAGAAACGGTAAAATGAAAAGAAAGGTAAATCCCGGCGGTTGGGTAATCCCCCCCAGCGTGGGGGTAAAGGTGGGGGGAATCAGAAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCggg                                                                                               <  1:368040/250‑1 (MQ=255)
                                                                                                                                                      gaTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGGAAATTCCAGCCGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGGTGCGGGAAAACGCAAAGCCCCCCCTTTTTCCCCACTTTTCGGTGGGGGCGGGGAAGAAAACGGCTCAAGAAAGCCCCCGGTACCCAAAAATTTTTGGGTGCCTGTGCGGGCGAAACCCCCTTCCcggg                                                                                              >  2:240302/1‑248 (MQ=255)
                                                                                                                                                         aCAGTAGGGGCGACATGGAGGCCGCCAGCGGGGGGTGTCCCGCAAACAGAAAGGGAAAAAGAAACAGGAAGGAAAATCCGGGGGGTTCGGTAATCCCCCCCAGCGTGGGGGTAAGGGTGGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGaacaa                                                                                          <  1:120486/250‑1 (MQ=255)
                                                                                                                                                                               cgccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAGGGTGCGGGGAACAGCAAGACCGCCACTTTTACCCCAGTTTTCCGTGCGGCAGTGGAGTACATCGCCCAAGAAATGCCCCCCGCCCCAAAGAATTTTGTGTTTCCCGGCAGGGGGAAAACCCCCTCTCGGAAGAAACATTTTAGCGGCTCTGgcgcc                                                                    >  2:519812/1‑249 (MQ=255)
                                                                                                                                                                                    gccgggCGGGTACCGAAAACAACAAGGGAAAAAAGAACGAGAGGGAAAATCCCGGCGGTCCGGAATCCCCCACCAGCGGGGGGTTTAGGGGTGCGGAAACCACCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACt                                                              <  2:674220/246‑1 (MQ=255)
                                                                                                                                                                                     cACGGCGTGTACCGCAAACAGCAGCGGTGAAAAGAACAGGAAGGTAAATTCCGGCGGTTCGGTAATCCCCACCAGCTGGGGGGTTAAGGTTGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTg                                                             <  1:340713/250‑1 (MQ=255)
                                                                                                                                                                                                            agcgagaaaaaGAAGAAGAAGAAAAATCCCGGCGGTCGATAACTCCCACCCAGCATGGCGGTAAAGGTGGCGGGAAACAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGccc                                       <  2:303989/243‑1 (MQ=255)
                                                                                                                                                                                                                   aaaTGAACAGGAAGGTAAATCCCGGCGGTTCGGTAATCCCCACCAGCATGGCGTTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACAt                                <  1:105188/250‑1 (MQ=255)
                                                                                                                                                                                                                         aCAGGAAAGTAAATTCCAAGGTGTTGGTTAATCCCCCCCACCAGGGCGTTTAGGGTGGGGGAAATAAACAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGaa                           <  2:558256/250‑1 (MQ=255)
                                                                                                                                                                                                                           aGGAAGTTAAATTGCAGCGTTCCGCTAACCCCCACCAGAGGGGGGGTGAAGGTGGGGGAAATCACCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATgc                        <  2:17293/250‑1 (MQ=255)
                                                                                                                                                                                                                                gtaaaaTTCCGGGGTTCGGGAAATCCCCCCCAGCTTGGGGGTAAAGGTGGCGGGAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCt                   <  1:610912/247‑1 (MQ=255)
                                                                                                                                                                                                                                   aaaTCCCAGCGGTTCGTAACTCCCCCCCAGCAGGGGTGTTAAGGTTGGGGAAAACAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTAAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCa                <  1:690830/250‑1 (MQ=255)
                                                                                                                                                                                                                                     atcccAGCGGTTCGAAAATCCCCACCAGCATGGCGGTTAGGGTGGCGGAAATAAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAAc              <  1:143860/247‑1 (MQ=255)
                                                                                                                                                                                                                                                 ttcgtaaaTCCCCACCAGCATGGGGTTAAAGGTTGCGGAAATCACCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGa   <  1:656986/244‑1 (MQ=255)
                                                                                                                                                                                                                                                    ggTAATTCCCACCAGCATGGCGGTTAAGGGTTCGGGGATCAGCAAGCCCGCCCCTTTTACCCGATTTTTCGGGGCGGGAGTGAAGGACATCGCTAAAGAAATGCCCCCGGCACCCAAAATTTTTTAGTTACCGTGCCGGGCAAAACCGCCCTTCCGGGAACACGATTTTAGCGGCTCGGGGCGCCAAACGAACTCCTGCCGAAGCTTGGCCCCGTACATCTGGATGGCGCCCTCCACAGCCTGCCGGAc  >  2:139669/1‑249 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                            
GGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAC  >  NZ_CP009273/3855751‑3856239

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