Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,621,9790TA28.3% 106.0 / 12.3 54T107S (ACT→TCT) pdeFcyclic‑guanylate‑specific phosphodiesterase PdeF
Reads supporting (aligned to +/- strand):  ref base T (30/8);  new base A (0/15);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 7.84e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACCCGAATGATGAAGTAGCACAGCGGGACACCGATCAGATTACCCACCAGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGC  >  NZ_CP009273/2621732‑2622220
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
aCCCGAATGATGAAGTAGCACAGCGGGACACCGATCAGATTACCCACCAGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTa                                                                                                                                                                                                                                                    >  1:228060/1‑249 (MQ=255)
  ccGAATGATGAAGTAGCACAGCTGGACACCGATCAGATTACCCGCCAGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAg                                                                                                                                                                                                                                                 <  2:385989/250‑1 (MQ=255)
    gAATGATGAAGTAGCACAGCGGGACACCGATCAGATTACCCACCAGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑c                                                                                                                                                                                                                                              >  1:474822/1‑250 (MQ=255)
       tgatgaAGTAGCACAGCGGGACACCGATCAGATTACCCACCAGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCt                                                                                                                                                                                                                                           >  2:299450/1‑250 (MQ=255)
           gaAGTAGCACAGCGGGACACCGATCAGATTACCCACCAGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTaa                                                                                                                                                                                                                                      >  1:42990/1‑250 (MQ=255)
              gTAGCACAGCGGGACACCGATCAGATTACCCACCAGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCg                                                                                                                                                                                                                                   <  2:410348/250‑1 (MQ=255)
                                 aTCAGATTACCCACCAGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGaatgataa                                                                                                                                                                                                                >  1:560724/1‑250 (MQ=255)
                                         aCCCACCAGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATgataatgataat                                                                                                                                                                                                         >  1:637739/1‑249 (MQ=255)
                                            cccaagaaaGGCCAGATAATTAATTAAAGTCCCAAGTTAAAAAGCCATAACACCAACCAGATTTTCTCTGCGCGCCAGAAATCCTACAAATGCAACAAACGGAAAAAGTATCAAAAACAGCTTTGCAGGAAACACAATCTGCCAAAATATACGCTGGGAAAATAAAC‑GGTTATCCCCATGTGAAACGTTTTTACGCCGGGTAGTAAAGC‑CCCT‑GTATCCGCCCCAACAGAGAATGATAATGATAATaaaa                                                                                                                                                                                                     <  1:507961/243‑1 (MQ=255)
                                               caGTAAGCCCTGATAATTGATTAAAGCCCCGAGTTAAAAAGGAATACCCCCGACCATATTTTCTCTGCTCGCAAGTAATCCAAAAAATACAGAAAACTGAAAAATTAAAAGAAACAGCGTTGCAGGAAACACATTCTGCAAAAATATACGC‑GGGAAAATAAACGGGGTATCACCTGGTAAAACGTTTTTACGCGGGGGATAAAAGC‑CCTT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTa                                                                                                                                                                                                  <  2:578468/250‑1 (MQ=255)
                                                   aaGGCCTGATAATTGATTACAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATCCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCCGCCAGAATATACGCTGTGAAATTAAAC‑GGGTATCACCATGAGAAACGTTGATACGCCGGGCAGTAAAAA‑CCCT‑CTAACCGCCCCTACAGTGAATGTTAATGTTAATAAAATAAGcag                                                                                                                                                                                              >  2:491078/1‑248 (MQ=255)
                                                          gATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAACCAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCAACATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAAGTAGCCGTTAATGa                                                                                                                                                                                       >  2:136032/1‑250 (MQ=255)
                                                           aaaaTTAATAAAAGTCCCGAGGTAAAAAGCCATGACACCGACCAAATTTCCTCTGCCCGCCAGAAATCCTACAAAGGAAACAAACTGAAAAAGTAACAAAAACAGCGTGGCAGAAAACAAAATCTGCAAAAATATCCGCTGGAAAATAAAAC‑GGGTACCACCATGGAAAACGTTGTTACGCGGGGGAGAAAAGC‑CCCT‑GAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATg                                                                                                                                                                                        <  1:30354/246‑1 (MQ=255)
                                                             aaTTGATTAAAGCCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCCCGCCAGTAATCCTACAAATGCAGAAAACTGAAAAAGTATCAAAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATACGCGGGAAAATAAAAC‑GGGTACCCCCATGTGAAACGTTGTTACGCGGGGGAGTAAAGC‑CCCT‑GAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAg                                                                                                                                                                                    <  1:472021/250‑1 (MQ=255)
                                                                 gATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGtt                                                                                                                                                                                >  1:370304/1‑250 (MQ=255)
                                                                     aaaGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCg                                                                                                                                                                             >  1:631105/1‑249 (MQ=255)
                                                                       aGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTa                                                                                                                                                                           >  1:261454/1‑249 (MQ=255)
                                                                               ggTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                                                                                                                                                  >  1:649559/1‑250 (MQ=255)
                                                                                gTTAAAAGGCATGACACCGCCCAATTTTTCCCTGCTCGCCAGTAATCCAAAAAATGCAGCAAACTGAAAAAGTATAAAAAACAGCGTGGCAGAAAACACAATCTGCCAGAATATACGTTGGAAAATAAAAC‑GGGTACCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGC‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGc                                                                                                                                                                 <  1:573534/250‑1 (MQ=255)
                                                                                            gACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTg                                                                                                                                                     >  1:530891/1‑250 (MQ=255)
                                                                                                ccGACAAGATTTCCTCGGCTCGCAAGTAACCCTAAAAATACAGCAAACTGAAAAAGTAACAAAAACAGGGTTGCAGAAAACACAATCTGCAAGAATATCCGCGGGAAAATAAAAC‑GGTTATCACCAGGAAAAATGTTGTTACGCGGGGGATAAAAGC‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGa                                                                                                                                                 <  2:228060/250‑1 (MQ=255)
                                                                                                       gATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAAt                                                                                                                                          >  1:60328/1‑250 (MQ=255)
                                                                                                             cgtgccTCGCAAGAAATCCAAAAAATCAAGAAAACGGAAAAAGTAACAAAAACAGCTTTACAGAAAACACAATCTGCAAGAATATACGCTGAAAAATAAAGC‑GGGTATCACCATGTAAAACGTTGTTACGCCGGGGAAAAAAGA‑CCTG‑TTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTa                                                                                                                                     <  2:114278/244‑1 (MQ=255)
                                                                                                                       cAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGAGCCCt                                                                                                                          >  1:290804/1‑250 (MQ=255)
                                                                                                                                 aCAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATgccg                                                                                                                >  2:631131/1‑250 (MQ=255)
                                                                                                                                  cAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATgccgc                                                                                                               >  1:541563/1‑250 (MQ=255)
                                                                                                                                  cAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GAGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATgccgc                                                                                                               >  1:354632/1‑250 (MQ=255)
                                                                                                                                        cagcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCt                                                                                                         >  1:175059/1‑250 (MQ=255)
                                                                                                                                             aaCGGGAAAAGTAAAAGAAACAGCGTGGAGGAAAACACAATCGGCAAGAATATCGGCTGGAAAATAAAAG‑GGGTATCACAATGTAAAATGTTGTTCGGCGGGGAAAAAAAGC‑CCCG‑GAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaaa                                                                                                    <  2:210992/250‑1 (MQ=255)
                                                                                                                                              aCTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCAAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaaaa                                                                                                   >  1:638850/1‑250 (MQ=255)
                                                                                                                                                    aaaGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatcat                                                                                             >  2:420417/1‑250 (MQ=255)
                                                                                                                                                               aaCAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTg                                                                                                                                                             >  1:571363/1‑175 (MQ=255)
                                                                                                                                                                acggCGTGGAAGAAAACACAAGCTGCAAAAATACACGCGGGAAAATAAAAC‑GGGTATCCCCATGAAAAAGTTTGTTACGCCGGGGAGAAAAGC‑CCTT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGatcat                                                                                 <  2:151181/247‑1 (MQ=255)
                                                                                                                                                                 cAGCGTTGAAGGAAACACAATCTGCCAGAATATACGTTGGGAAATTAAAC‑GGTTATCCCCATGTGAAACTTTGTTACGCCGGGGATTAAAGA‑CCCA‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGatcatc                                                                                <  2:281333/250‑1 (MQ=255)
                                                                                                                                                                  aGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGatcat                                                                                 >  1:317277/1‑248 (MQ=255)
                                                                                                                                                                          ggaaaaCAAAACCTGCAAGAATAGACGCGGGAAAATAAAAC‑GGGTATCACCTGGTAAAACGTTGTTCCGCGGGGGAG‑AAAAA‑CCCTGGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCaga                                                                       <  2:225419/247‑1 (MQ=255)
                                                                                                                                                                                aaaaaTCTGCCAAATTAT‑CAGTGGGAAAATAAAA‑CGGTCTCCCCATTGGAAAAGGTGTTAACGCGGGGAATAAAAACCCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGa                                                                 <  2:458194/247‑1 (MQ=255)
                                                                                                                                                                                    aTCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCCGATaa                                                              >  2:466263/1‑249 (MQ=255)
                                                                                                                                                                                      cTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATaa                                                              >  1:65456/1‑247 (MQ=255)
                                                                                                                                                                                         ccAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAagag                                                         >  2:219092/1‑249 (MQ=255)
                                                                                                                                                                                         ccAGAATATACGCGGAGAAATAAAAC‑GGATATCCCAAGGTGAAACTTTTTTACGCGGGGGAGAAAACC‑CCCT‑ATAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAgaga                                                        <  2:598124/250‑1 (MQ=255)
                                                                                                                                                                                                 aaggcgggaaaaTAAAAC‑GGTTATCACAATGTGAAACGTTTTTCCGCCGGGGAAAAAAGC‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTc                                                <  1:151098/241‑1 (MQ=255)
                                                                                                                                                                                                         aaaaTAAACC‑GGTTACCCCCATGGAAAACGTTGTACCGCGGGGAAGAAAAAA‑CCTT‑TAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGgaag                                        <  2:198025/249‑1 (MQ=255)
                                                                                                                                                                                                          aaaTTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGgaag                                        >  2:190092/1‑249 (MQ=255)
                                                                                                                                                                                                           aaTTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGgaaga                                       >  2:236168/1‑249 (MQ=255)
                                                                                                                                                                                                            aTTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGgaagaa                                      >  2:451528/1‑249 (MQ=255)
                                                                                                                                                                                                               aacc‑GGGTACCCCCATGAAAAACGTTTTTACGCCGGGGAGAAAAGC‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAatat                                  <  1:158497/247‑1 (MQ=255)
                                                                                                                                                                                                                  c‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGGGTAAAGA‑CTCT‑GTAATCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTTGTAAAAACCCAGCTCTGCATATTTGGGGACAAAAATCTCTAAGGCGATCCCTGGCAATGCCGCTCAGCTGAAAAACATCATCTGGCTGATCATCAGCGCCAGAGGTGTTCAAAAGTAAAAAACCTTACTGGTAGAAATATgag                               >  2:2216/1‑248 (MQ=255)
                                                                                                                                                                                                                                   aaaGGTTGTTACGCGGGGGAAAAAAGA‑CCCT‑AAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAAtgtg               <  2:492743/250‑1 (MQ=255)
                                                                                                                                                                                                                                    agcGTTGTTACGCCGGGGAGAAAAGC‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTgg              <  1:502029/248‑1 (MQ=255)
                                                                                                                                                                                                                                        tttttACGCCGGGAAGAAAACC‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCg           <  1:361717/246‑1 (MQ=255)
                                                                                                                                                                                                                                         tgtaccGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGa          <  1:210330/244‑1 (MQ=255)
                                                                                                                                                                                                                                           tgaCGCCGGGGGAAAAAGC‑CCCT‑GAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGaaaa       <  2:7467/248‑1 (MQ=255)
                                                                                                                                                                                                                                                cgggggAGAAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGc  <  2:42990/248‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
ACCCGAATGATGAAGTAGCACAGCGGGACACCGATCAGATTACCCACCAGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGA‑CCCT‑GTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGC  >  NZ_CP009273/2621732‑2622220

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