Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,089,5710AG37.1% 117.3 / 22.2 70I83V (ATT→GTT) hisFimidazole glycerol phosphate synthase subunit HisF
Reads supporting (aligned to +/- strand):  ref base A (9/35);  new base G (26/0);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.26e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGCTGGCAAAACGCATAATCCCATGTCTCGACGTTCGTGATGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGACGCCGAAACCGGTAAATATCATGTGAATCAATATACCGGCGATGAAAGCCGCACCCGCGTCACTCAATGGGAAACGCTCGACTGGGTACAG  >  NZ_CP009273/2089326‑2089816
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                          
tGCTGGCAAAACGCATAATCCCATGTCTCGACGTTCGTGATGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGtt                                                                                                                                                                                                                                                          >  1:707031/1‑248 (MQ=255)
   tGGCAAAACGCATAATCCCATGTCTCGACGTTCGTGATGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTa                                                                                                                                                                                                                                                         >  2:2432/1‑246 (MQ=255)
     gCAAAACGCATAATCCCATGTCTCGACGTTCGTGATGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGtctc                                                                                                                                                                                                                                                   <  1:492692/250‑1 (MQ=255)
             cATAATCCCATGTCTCGACGTTCGTGATGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGAAGAt                                                                                                                                                                                                                                           >  1:405133/1‑250 (MQ=255)
                     cATGTCTCGACGTTCGTGATGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAGGGCTCTGGGAGATGCcccgc                                                                                                                                                                                                                                   >  2:955929/1‑248 (MQ=255)
                       tGTCTCGACGTTCGTGATGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGaaaa                                                                                                                                                                                                                                <  2:342823/250‑1 (MQ=255)
                          ctcGACGTTCGTGATGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGGGGAGGGGGGGAGTAAGGGTCGGGGAGAGGACGGCGAAAAta                                                                                                                                                                                                                              >  2:725782/1‑249 (MQ=255)
                             gACGTTCGTGATGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTc                                                                                                                                                                                                                             <  1:495407/247‑1 (MQ=255)
                             gACGTTCGTGATGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCtt                                                                                                                                                                                                                           <  2:628222/249‑1 (MQ=255)
                             gACGTTCGTGATGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCAGGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCACGAAAATTCt                                                                                                                                                                                                                            <  2:34696/248‑1 (MQ=255)
                               cGTTCGTGATGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTATGCATCTTGGAGATGCCCCCGAAATTCctttc                                                                                                                                                                                                                        >  2:929412/1‑248 (MQ=255)
                                    gTGATGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTg                                                                                                                                                                                                                    <  2:970649/249‑1 (MQ=255)
                                     tGATGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGGTTAGTCGCCGGGAGATGCCGCGACAATTTTTTCCTTTTgc                                                                                                                                                                                                                  >  2:366902/1‑250 (MQ=255)
                                         ggTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTAAGTCTCTGGAAGATGCCCCGCAAATTCTTTCCTTTGGCGCgg                                                                                                                                                                                                              >  1:19046/1‑250 (MQ=255)
                                           tCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGAt                                                                                                                                                                                                            <  2:86477/250‑1 (MQ=255)
                                             aGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATaa                                                                                                                                                                                                          <  1:851055/250‑1 (MQ=255)
                                               gtggtgAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATaaaa                                                                                                                                                                                                        <  2:207312/250‑1 (MQ=255)
                                                              cAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACtctc                                                                                                                                                                                         <  1:94414/250‑1 (MQ=255)
                                                               aGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTcc                                                                                                                                                                                        <  1:531024/250‑1 (MQ=255)
                                                                   tCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTgc                                                                                                                                                                                     <  2:486105/249‑1 (MQ=255)
                                                                      cAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTg                                                                                                                                                                                <  1:431883/250‑1 (MQ=255)
                                                                      cAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTg                                                                                                                                                                                <  1:304037/250‑1 (MQ=255)
                                                                      cAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTTCTTTGGCGCGGATAAAATTTTCATCAACTCTCCTGCGCTg                                                                                                                                                                                >  1:278930/1‑250 (MQ=255)
                                                                         ccATGAAATCATTGGCGATATCGTGCCGCTGGAAAAACGCTACGCTGAGGAAGGCGCTGACGAAATGGTGTTCTACGATATCACCGTTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCGTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCa                                                                                                                                                                             <  2:374001/250‑1 (MQ=255)
                                                                           aTGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGAAGATGCCGCGGAAAATCTTTTCTTTTGCGCGGGGAAAATTTCCCTCAAACCTCCCGCGCTGGCaaa                                                                                                                                                                           >  1:709315/1‑248 (MQ=255)
                                                                             gAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGATTAAGTCTCTGGGAGATGGCGCGAAAATTCTTTCCCTTTGCGCCGGTAAAAATTCCCTCAACCCTTCCGCGCTGGCCGAcc                                                                                                                                                                        >  1:188705/1‑250 (MQ=255)
                                                                                aTCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTTCTTTTGCGCGGATAAAATTTCCATCAACTCTCCTTCGCTGGCAGACCCa                                                                                                                                                                      >  1:350040/1‑249 (MQ=255)
                                                                                    ttGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACatt                                                                                                                                                                 <  1:508653/250‑1 (MQ=255)
                                                                                    ttGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTAAGTCTCCGGAAGAAGGCCCGAAAAATTTTTTCCTTGGCGCGGGTAAAAATTTCCTCAACTCTCCCGCGCTGGCGCACCCCACat                                                                                                                                                                  >  2:270896/1‑249 (MQ=255)
                                                                                     tggagttaccGTGCCGCGGACAAAACGCGACGGTGAAGAAGGGGCTGACGAACTGGTATTCTTCGATATCACCGTTGCCAGCGGTGGCCGTGTGGTAGATAAAAGCTGGGTGTCTCGCGTGGCGGAAGTGATCGACATTCCGTTGTGTGTGGCGGGTGGGATTAAGGCTCTGGAAGATGCCGCGAAAATTCTTTCCGATGGCGCGGATAAAATTTCCATCAACGCTCCTGCGCTGGCAGACGCAACatta                                                                                                                                                                <  2:633078/241‑1 (MQ=255)
                                                                                       gCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACattaa                                                                                                                                                               <  2:311761/249‑1 (MQ=255)
                                                                                        cGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACattaatt                                                                                                                                                             <  2:916693/250‑1 (MQ=255)
                                                                                             tCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTc                                                                                                                                                         <  1:536598/249‑1 (MQ=255)
                                                                                              cGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGATTTAATCTCTGGGAGATGCCGCGGAAATTCTTTTCTTTTGCGCGGATAAAATTTTCATTAACTCTCCTTCGCCGGCAGACCCAACATTAATTACTCGc                                                                                                                                                       >  1:667620/1‑250 (MQ=255)
                                                                                                     cTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTTGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCg                                                                                                                                                >  1:302002/1‑250 (MQ=255)
                                                                                                            aaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCCCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTTAATCTCTGGGAAATGCCGCGGAAATTTTTTTCTTTGGGGCGGGTAAAATTTTCATTAACTCTCCCTCGCTGGGCGACCCCACAATAAATTCCCCCCTTGCCGAACGcct                                                                                                                                         >  2:800244/1‑248 (MQ=255)
                                                                                                                        aagaagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAgtg                                                                                                                             <  2:59597/250‑1 (MQ=255)
                                                                                                                           aagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGcc                                                                                                                                                 <  1:376158/227‑1 (MQ=255)
                                                                                                                              gcgaTGACGACCGGGTGTGCTACGATATCAACGCTTCCAGCGATGGTCGTGTGGTAGATAATAGCTGGGTATCTCGCGTGGCCGAAGTGATCGACATTCCGTTGTGTGTGGCGGGTGGGAGTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGAATAAAGTGTCCATCAACTCTCCTGCTCTGGCAGACCGAGCATGAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATtgt                                                                                                                       <  2:950034/250‑1 (MQ=255)
                                                                                                                               cgcTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATtgtg                                                                                                                      <  2:937859/250‑1 (MQ=255)
                                                                                                                                 cTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTTAGTCTCTGGGAGATGCCGCCGAAATTCTTTCCTTTTGCGGGGGTAAAATTTTCATCAACCCCCCTGCGCTGGCAGACCCAACATTAATTTCTCCCCCGGCCCACCCCTTTTGCGTGGCGGGTAtttgggg                                                                                                                    >  2:780937/1‑245 (MQ=255)
                                                                                                                                  tGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTTGAAGATGCCGCGGAAATTCTTTCCTTTTGGGCGGGTAAAAATTTCCTCAACTCTCCCTCGCTGGCAGAACCCACATTTAATAATCGCCTGGCCGCTCGCCTTGGCGGGCCGGGGATTTTGGTc                                                                                                                   >  2:167748/1‑250 (MQ=255)
                                                                                                                                   gACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTTCTTTTGCGCGGATAAAATTTTCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCCCCCGGCCGATCCCTTTGGCGTGCAGTGTATTTTGGTCg                                                                                                                  >  1:829884/1‑250 (MQ=255)
                                                                                                                                     cGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGAAGATGCCGCGAAAATTTTTTCCTTTGGCGCGGGTAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGGGTATTTTGGTCGGt                                                                                                                >  1:125938/1‑250 (MQ=255)
                                                                                                                                      gAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTTGCGCGGATAAAAATTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTTGCGTGCAGTGTATTTTGGTCggga                                                                                                               >  1:246519/1‑248 (MQ=255)
                                                                                                                                        aCTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTAtt                                                                                                             <  1:234539/250‑1 (MQ=255)
                                                                                                                                          tgaagtTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCCGAACTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGGCTCACGGAGATGCCGCGCAAAATTTTTCCTTTTACGCGGGGAAAATTTCCCTCCACTCCCCTGCGCTGGGAGACCCAACAATAAATACAACACTCGCCGCGCCCTCTGGGGTGCAGCGCATTGTGGTCTGTATTgt                                                                                                           >  2:398121/5‑249 (MQ=255)
                                                                                                                                                   aCGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTTAGTCTCTGGAAGAAGCCCCGAAAATTTTTTTCTTTGGCGCGGGTAAAATTTTCATCAACTCTCCCTCGCTGGCAGACCCAACATTTATTACTCTCCTGGCCGGTCGCCTTGGCGTGCAGTGTATTTTGGGCGGTATTTGTACCTGGTa                                                                                                  >  2:309656/1‑250 (MQ=255)
                                                                                                                                                    cGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTa                                                                                                  <  2:374426/249‑1 (MQ=255)
                                                                                                                                                       tatCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTTAGTCTCTGGGAGATGCCCCGCAAATTTTTTCCTTTTGCGCGGGTAAAATTTTCATCAAATCTCCTGCGCCTGGCGACCCCAAATTAATTACTCGCCTCGCCCGTCGCCTTTGCGTGCCGGGTATTGTGGTCTTTGTTTTTCACTGGTacgac                                                                                              >  2:577618/1‑250 (MQ=255)
                                                                                                                                                          cACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGAAGATGCCGCGGAAATTTTTTTCTTTTGCGCGGGGAAAATTTTCATCAACTCTCCTGCGCTGGCCGACCCCACATTAATTAATCGCCCGGCCGGTCGCTTTTGGGGGCAGTGGATTGTGGGCGGTATTGATAACTGGGAACACGcc                                                                                           >  2:859944/1‑250 (MQ=255)
                                                                                                                                                            ccGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGACGCCGa                                                                                        <  2:931729/250‑1 (MQ=255)
                                                                                                                                                                  ccAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGACGCCGAAAcc                                                                                    <  1:747713/248‑1 (MQ=255)
                                                                                                                                                                         tGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGACGCCGAAACCGGTAAatat                                                                          <  2:248518/250‑1 (MQ=255)
                                                                                                                                                                                 tgGTAGATAAAAGCTGGGGATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGACGCCGAAACCGGTAAATATCATGTGaa                                                                  <  2:167511/250‑1 (MQ=255)
                                                                                                                                                                                      gATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGACGCCGAAACCGGTAAATATCATGTGAATCaa                                                              <  1:900547/249‑1 (MQ=255)
                                                                                                                                                                                       aTAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGACGCCGAAACCGGTAAATATCATGTGAATCAata                                                            <  2:261440/250‑1 (MQ=255)
                                                                                                                                                                                        tAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTTAATCTCTGGGAGAAGGCGGGAAAATTTTTTTCTTTTGGGGGGGGAAAATTTTCATCAAACGTCCTGCGCTGGGGGCACCCACCATAAATTCTCGCCTGGGCGAGCGCTTTTGGGTGCAGGGTATTTTTGTTGGGAATGGTAACTGGTACGACGCCGGAAACGGGGAATATCAAGTGAATTaaaa                                                            >  2:973823/1‑247 (MQ=255)
                                                                                                                                                                                                      ctcGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTAAGGCTCTGGGAGATGGCGCGGAAAATTTTTTCTTTGGGGCGGAGAAAAATTTCATCAAATCTCCTGCGCTGGCAGAACCAACATTAAATAATCGCCTGGCCGCTAGCCTTGGCGTGCCGGGGATTTTGGGGGGGATTTATACCTGGTACGGCGCCGGAACGGGGGAATAAACATGGAATCAATATACAGGCgcagga                                               >  2:733177/1‑243 (MQ=255)
                                                                                                                                                                                                            tGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAAGCTCTGGGAGATGCCGCGGAAAATCTTTTCTTTGGGGCGGGTAAAAATTTCATCAACCCTCCTGCGCTGGCAGAACCAAAATTAATTACTCGCCTGGCCGACCCGTTTGGGGTGGAGGGGAATGTGGGCGGGATTGAAACCTGGTAAGGCGCCGAAAACGGGAAAAATCATGTGAAACAATATACCGGCGGAGGAAGCCGCAc                                       >  2:914084/1‑250 (MQ=255)
                                                                                                                                                                                                              gCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGACGCCGAAACCGGTAAATATCATGTGAATCAATATACCGGCGATGAAAGCCGCAccc                                     <  1:615919/250‑1 (MQ=255)
                                                                                                                                                                                                               cGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTGGAAGATGCCGCGAAAATTTTTTTCTTTTGCGCGGATAAAATTTTCATCAACTCTCCTGCGCTGGCAGACCCCACATTAATTACTCGCCTTGCCGATCGCTTTGGCGTGCAGTGTATTGTGGGCGGTATTTATACCTGGTACGACGCCGAAACCGGGAAATATCATGTGAATCAATATAACGGCGATGAAAGCCGCACCcg                                    >  1:694923/1‑250 (MQ=255)
                                                                                                                                                                                                                   aGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGACGCCGAAACCGGTAAATATCATGTGAATCAATATACCGGCGATGAAAGCCGCACCCGCGTc                                <  2:465821/250‑1 (MQ=255)
                                                                                                                                                                                                                     tGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTGGAAGATGCCGCGGAAATTTTTTCCTTTGGGGCGGGTAAAATTTTCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTTGCGTGCAGGGTATTTTGGTCGGGATTTATACCTGGTACGGCCCCGAAACCGGTAAAAATCATTTGAATCAATATAACGGCGGTGAAAGCCCCCCCCCCGTccc                              >  2:499078/1‑248 (MQ=255)
                                                                                                                                                                                                                       aTCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAAGCTCTGGGAAATGCCGCGAAAATTCTTTTCTTTGGCGCGGGTAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTTATCGCCTGGCCGCTAGCTTTTGCGTGCAGTGTATTTTGGTCGGTATTTATACCTGGTACGACGCCGCAAACCGGAAATATCATGTGAATAAATATACCGGCGGTGAAAGACGCACCCCCGCCACt                             >  2:945799/1‑249 (MQ=255)
                                                                                                                                                                                                                       aTCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGGAGATGCCGCGAAAATTCTTTTCTTTGGCGCGGGTAAAATTTTCATCAACTCTCCTGCGCTGGCAGACCCCACATTAATTAATCGCCTGGCCGATCTCTTTGGCGTGCAGAGTATTTTGGGCGGTATTGATACCTGGTACGACGCCGAAAACGGAGAAAATCATGTGAATAAATATACAGGGGCTGAAAGCCCCCCCCGCGACCCTc                            >  2:84183/1‑250 (MQ=255)
                                                                                                                                                                                                                          gACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGACGCCGAAACCGGTAAATATCATGTGAATCAATATACCGGCGATGAAAGCCGCACCCGCGTCACTCaa                          <  1:769389/249‑1 (MQ=255)
                                                                                                                                                                                                                                      tgtgtgGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTTCTTTGGCGCGGATAAAATTTTCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGACGCCGAAACCGGTAAATATCATGTGAATCAATATACCGGCGATGAAAGCCGCACCCGCGTCACTCAATGGGAAACGCTc              >  1:181649/1‑249 (MQ=255)
                                                                                                                                                                                                                                                 tGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGACGCCGAAACCGGTAAATATCATGTGAATCAATATACCGGCGATGAAAGCCGCACCCGCGTCACTCAATGGGAAACGCTCGACTGGGTACAg  >  1:121730/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  gggATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGACGCCGAAACCGGTAAATATCATGTGAATCAATATACCGGCGATGAAAGCCGCACCCGCGTCACTCAATGGGAAACGCTCGACTGGGTACAg  >  1:482234/1‑249 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                          
TGCTGGCAAAACGCATAATCCCATGTCTCGACGTTCGTGATGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGACGCCGAAACCGGTAAATATCATGTGAATCAATATACCGGCGATGAAAGCCGCACCCGCGTCACTCAATGGGAAACGCTCGACTGGGTACAG  >  NZ_CP009273/2089326‑2089816

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