Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,372,4470CT32.7% 71.1 / 19.8 56E98K (GAA→AAA) menD2‑succinyl‑5‑enolpyruvyl‑6‑hydroxy‑3‑ cyclohexene‑1‑carboxylic‑acid synthase
Reads supporting (aligned to +/- strand):  ref base C (6/31);  new base T (18/0);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 9.34e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCATTTCGCCATACAGCGGTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATGTGTCTGACGCCGTGACGCGTTA  >  NZ_CP009273/2372204‑2372686
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
ccATTTCGCCATACAGCGGTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCg                                                                                                                                                                                                                                                   <  2:731072/248‑1 (MQ=255)
  aTTTCGCCATACAGCGGTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTa                                                                                                                                                                                                                                               <  1:471169/250‑1 (MQ=255)
          aTACAGCGGTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGGTTAGAGAGACAGGTTTTTTCCCGTTAACACgggc                                                                                                                                                                                                                                      >  2:20265/1‑248 (MQ=255)
               gCGGTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTACGCGGAATAGGTTTTTTTCCGGTTTCCCCgcgctcac                                                                                                                                                                                                                                  >  2:691478/1‑243 (MQ=255)
                cGGTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGGTAAGGGAATTCGGTTTTTTCCGGGTTACCCGGCTTCAAt                                                                                                                                                                                                                                 >  1:225276/1‑249 (MQ=255)
                 ggTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGCAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGGGGAGGGGATAAGGTTTTTTCCGGGTTAAcccagtgccaa                                                                                                                                                                                                                                   >  2:287776/1‑240 (MQ=255)
                  gTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTTAGAGAATCAGGTTTTCTCCGGTTTAACCGGCTTCAATCaa                                                                                                                                                                                                                              >  1:148508/1‑249 (MQ=255)
                          aaCGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGAt                                                                                                                                                                                                                     <  1:202478/250‑1 (MQ=255)
                          aaCGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGAt                                                                                                                                                                                                                     <  2:713274/250‑1 (MQ=255)
                              ggCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATagag                                                                                                                                                                                                                 <  2:741700/250‑1 (MQ=255)
                                     gATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTTAGAGAATCAGTTTTTTTTCGGGTAACCCGGCTTCAATCAGTGCCGGGTAGAGATTTGcc                                                                                                                                                                                                        >  1:175101/1‑250 (MQ=255)
                                         tGGACTCCCGCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCAccg                                                                                                                                                                                                    <  2:559002/250‑1 (MQ=255)
                                         tGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCAcc                                                                                                                                                                                                     <  1:405895/249‑1 (MQ=255)
                                           gACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCAccgcc                                                                                                                                                                                                  <  1:660810/250‑1 (MQ=255)
                                             cTCCCCCGGCATGAGGCGTACGGAGAGCGTGGTCGAGGGAAGAAACCAGCCACCGGGCGGGGATATCCGGGGTCGGGCGCGGCAATGAAATACTGTGCGGGGGGGGAGGGGCGAACATCCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGTTTAAGAGAATCAGTTTTTCTCCGGTTACCCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGt                                                                                                                                                                                                <  2:205619/250‑1 (MQ=255)
                                             cTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAGCCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATCCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGt                                                                                                                                                                                                <  1:679730/250‑1 (MQ=255)
                                                 ccccGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTTAGAGAATCAGTTTTTCTTCGGGTTACCCGGCTTCAATTAGTGCCGGACAGAGAATTTCCACCCCCGTGCCg                                                                                                                                                                                            >  1:24166/1‑250 (MQ=255)
                                                  cccGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAAACAGTTTTTTTCCGGGTTAACCGGCTTCAATCAGTGCCGGGTAGAGATTTTCCACCCCCGCGCCgg                                                                                                                                                                                           >  1:281961/1‑250 (MQ=255)
                                                         gAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGGTAAGAGGATTAGTTTTTTTCCCGTTAAACCGGCCTCAAACAGTGGCGGCGAGAGATATTTCAACCCCGCGCCGgggggacc                                                                                                                                                                                    >  2:607290/1‑246 (MQ=255)
                                                           aGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCCGTTTTTTTTCGGTTTAACCGGCTTTCATTAAGGCCGGGTAGAGAATTTGCCACCCCCGGCCCGGGGGtacc                                                                                                                                                                                    >  1:704615/1‑246 (MQ=255)
                                                              gTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTTAGAGAATCCGGTTTTTTCCGGGTAAACCGGCTTCAATCAGTGCCGGGTAGAGAATTGTCACCCCCGTGCCGGGGGGCACAATCa                                                                                                                                                                               >  1:301107/1‑250 (MQ=255)
                                                                        cGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGGTAAGAGAATCAGGTTTTTTCCGGGTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGGGCCGGGGGGCACAATCACCGCCAcccg                                                                                                                                                                     >  1:415917/1‑248 (MQ=255)
                                                                         gTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGCGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGgc                                                                                                                                                                    <  1:640785/250‑1 (MQ=255)
                                                                          tgggcggggggaGAAACCAGCGAACGTGCGGGGATATCCGGGGGGGGGCGGGGCAATGAAAGACTGTGCGGGGGGTGAGAGAGAAATATTCCCGGCGGGCGAATTGCCTGATTCGCGCGGCGGTCAATTAGCGCCGGCGGGCGATGGGCGGTTAGGAGAATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGAATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGgct                                                                                                                                                                   <  2:251297/239‑1 (MQ=255)
                                                                           ggTCGATGGTAGAAACCAGCCACCGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGgctg                                                                                                                                                                  <  1:213171/250‑1 (MQ=255)
                                                                            gTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGgctg                                                                                                                                                                  <  1:168630/249‑1 (MQ=255)
                                                                                    tAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTTAGAGAATTAGTTTTTTTTCGGGTTACCCGGCCTCAATCAGTGCCGGATAGAGAATTTCCACCGCCCTGCCGGAGGTCACAATCACCGCCACCCGCTGCTTTCTGAc                                                                                                                                                         >  1:413378/1‑250 (MQ=255)
                                                                                          ccAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGc                                                                                                                                                   <  1:625917/250‑1 (MQ=255)
                                                                                                  cGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTTTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCTCCGCCGTGCCGGGGGGCCCAAACACCGCCACCCGCCGCTTGCTGGCTTTTCCCAAccccc                                                                                                                                           >  1:120987/1‑249 (MQ=255)
                                                                                                          gATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGGTTAGAGAAACCGGTTTTTTCCCGGTTAACCCGCCTCCATCCATGCCGGGTAGGGATTTTCCCCCCCCCTGCCCGGGGGCACCATCACCCCCCCCCGCCGCCTGCTGACTTTCGCCCGcccccccccccgg                                                                                                                                   >  2:322710/1‑241 (MQ=255)
                                                                                                                           gcgGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATGGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTc                                                                                                                  <  1:545702/250‑1 (MQ=255)
                                                                                                                            cgGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCa                                                                                                                 <  1:105660/250‑1 (MQ=255)
                                                                                                                            cgGCAATGAAAGGCGGTGCGGGGGGTGAGGGGCGAACATTCCCGGCTGGCAAATTGCCTGATCCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCa                                                                                                                 <  2:697459/250‑1 (MQ=255)
                                                                                                                                    aaaTACTGTGCGTGGGGTGAGAGGCGACCATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAAt                                                                                                         <  1:458377/249‑1 (MQ=255)
                                                                                                                                    aaaTACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTAAATTAGCTCCGGCGGGCGATCGGCGGTAAAGAGAATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATg                                                                                                        <  2:529147/250‑1 (MQ=255)
                                                                                                                                         cTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGgtg                                                                                                    <  1:520763/249‑1 (MQ=255)
                                                                                                                                                    gggagaGGCGAACATTCCCGTCTGGCAATTTGCCTGATTCGCGCCGCAGGAAATTACCCCCGGCGGGCGATCGGCGGTTAAGAAAATCAGTTTTTCTCCGGTTACCCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTgaatgaat                                                                                        <  2:237983/248‑1 (MQ=255)
                                                                                                                                                               aCATCCCCGGCGGGCGAATTGCCTAATTCGCGCCGCAGTAAATTAGCGCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTACCCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATtctc                                                                             <  2:239567/250‑1 (MQ=255)
                                                                                                                                                                      cGGCTGGCGAATTCCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGTTTAACCCGGCTTCAACCAGTGCCGAATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTccgccgcc                                                                      <  2:320856/250‑1 (MQ=255)
                                                                                                                                                                                 ttGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTTGCCAGCCCCAGCGCCAGATGCCCCCACCCCCGCTCATCGAAAAGGGGGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTaa                                                          >  1:233363/1‑250 (MQ=255)
                                                                                                                                                                                  tgccgtaTTCCCCCCGCGGTCATTTAGCCCCGGCGGCCGACCGGCGTTTAAGAGAATCATTTTTTCTCCGGTTACCCCGGCTTCAATCAGTGCCGAATAGAGTTTCCCCCCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTaa                                                          <  2:44621/243‑1 (MQ=255)
                                                                                                                                                                                    ccTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTTTCCGGTTAACCCGGCTTCAATCAATGCCGGATAGAGATTTGCCACCGCCCTGCCCGAGGTCACAAACACCGCCCCCGGCTGCTTGCTGGCTTTCGCCCGCCCCCCCCCCCGGAGGCCCCACCCCCCCCCCTCCGAATGGGTGTGGTGGATTAATGCGGAATTCCCCCCCGCCCCCAACCTTTAccg                                                       >  2:637203/1‑248 (MQ=255)
                                                                                                                                                                                      tGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTTTCCGGTTAAAACGGCCTCAAACAGTGCCGGGGAAAGAATTGTCACCCCCGGGGCGGGGGGCACAATAACCGCCACCCGCTGCCTGCTGACTTTCGCCAGCCCCCGCGCCCGAAGCCCCAACCCCCCCGCACCGCAATGGGGGGGGGGAATGCAGTCGGGATTCTCCCCCTCCGCTAACGtgtaccgtg                                                     >  2:718010/1‑242 (MQ=255)
                                                                                                                                                                                              gccgcAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTTAGAGAAACAGTTTTTTTCCCGGTTAACCGGCCTCCATCAGTGCCGGGTAGAGATTTTCCCCCCCCGCGCCGGGGGGCACAAACACCGCCCCCCGGCGCTTGCTGACCTTCGCCAGCCCCCACGCCCGATGCCCCAACCCCCCCCCATCGAAATGGGGGTGGGGGAAGAAAGCGGGGATCTCCGCCCCCCCTAAAGTTTACGgggtgggaaac                                             >  2:698590/1‑243 (MQ=255)
                                                                                                                                                                                                cgcAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGa                                           <  1:20265/250‑1 (MQ=255)
                                                                                                                                                                                                     tCAATTAGCTCCGGCGGGCGATCGGCGGTTTAGAGGATCCGTTTTTTTCCCCTTAACCCGGCTTCAATCAGTGTCGGGAAGAGATTTTTCCCCCCCCTGTCGGGGGGCACCATTACCCCCACCCGGCCCCTGCCGGCTTTCCCCCGCCCCCACGCCCCATGGCCCCACCCCCCCCCCTCGCAACTGGGGGGGGGGAATAAATCCGGAATTTCCGCCCCCCCCAACCCTTAACGGTGAGGAAcgcggcgtg                                       >  2:108119/1‑245 (MQ=255)
                                                                                                                                                                                                      cAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGGTAAACCCGCTTCAATCAGTGCCGGATAGAGAATTGCCACCCCCGTGCCGGCGGGCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCCGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGGGTGGTGAATGAAAGCGGAAATCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTgg                                     >  1:696743/1‑250 (MQ=255)
                                                                                                                                                                                                       aaTTTGCCCCGGGGGGCCGTCGGGGGGTAAGAAAAACAGTTTTTCCCCGGTTAACCCGGCTTCAAACCGGGCCGGAAAAAGAGTTGCCCCCCGCCGTGCCGGAGGTCACAACCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCACTAACGTTAACGGTGTAGAACGCGAGCCTgg                                     <  2:175101/250‑1 (MQ=255)
                                                                                                                                                                                                           cgccccGGGGGGCGATCGGCGGTTAAGAGAATCATTTTTTCTCCGTTTACCCCGCCTTCAATCAGTCCCGGATAGAGATTTCCCCCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCg                                 <  2:415917/245‑1 (MQ=255)
                                                                                                                                                                                                           aGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTTTCCCGTTAACCCCGCCTCAATCAGTGGCGGGTAGAGATTTTCCACCGCCCTGCCGGGGGGCACAATCACCCCCACCCGCCGCTTTCTGACCTTTGCCCGCCCCCGCGCCAGATGGCCCCACCCCCGCTCCTCGCAAAGGGGGTGGGGGATGAATGCGGGAATCTCCCCCCCCCCCAACGGTTACGGGGTAGAACGCGAGGCTGGGGCgg                                >  2:80814/1‑249 (MQ=255)
                                                                                                                                                                                                                   cgggcgccgggcgtTAAAGAGAATGAGTTTTTGCCCGTGTACCCCGGCTCAATTAGGTGCCGAATAGAGTTTGCCCCCCGCCGTGCCGGAGGTCACAATCCCCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGAtg                        <  2:586569/241‑1 (MQ=255)
                                                                                                                                                                                                                                gTTAAGAGAATCAGTTTTTCTCCGGTTACCCCGGCTTCAATCAGTGCCGGATAGAGATTTCCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATGTGTCTGACGCCGt           <  1:737650/250‑1 (MQ=255)
                                                                                                                                                                                                                                    agagAATCGTTTTTTCCCCGGTTAACCCGGCTTAAATCAGTCCCGGATAGAGATTGGCCCCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATGTGTCTGACGCCGTGAcg       <  2:185003/250‑1 (MQ=255)
                                                                                                                                                                                                                                     gagaATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATGTGTCTGACGCCGTGAcgc      >  2:348376/1‑250 (MQ=255)
                                                                                                                                                                                                                                          aaTCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATGGGTCTGACGCCGTGACGCGt    >  2:708643/1‑249 (MQ=255)
                                                                                                                                                                                                                                            accatttttCCCCGGGTCACCCGGGCTCAAACCGGCGCCGAATGGAGTTTTGCCACCGCCGTGCCGGAGGTCCCAATCCCCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATGTGTCTGACGCCGTGACGCGTTa  <  2:18156/246‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
CCATTTCGCCATACAGCGGTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATGTGTCTGACGCCGTGACGCGTTA  >  NZ_CP009273/2372204‑2372686

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: