Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,652,4200CT33.9% 105.3 / 13.2 59S32F (TCT→TTT) ynfEDmsA/YnfE/YnfF family dimethyl sulfoxide reductase
Reads supporting (aligned to +/- strand):  ref base C (10/29);  new base T (20/0);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.43e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCA‑GGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTGAATTAC  >  NZ_CP009273/1652175‑1652667
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                           
aTACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTttttttttt                                                                                                                                                                                                                                                         >  2:331310/1‑249 (MQ=255)
      caTAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGgtt                                                                                                                                                                                                                                                   >  2:439013/1‑250 (MQ=255)
       aTAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTGCCGGTT‑a                                                                                                                                                                                                                                                 >  1:805101/1‑250 (MQ=255)
       aTAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGGTT‑a                                                                                                                                                                                                                                                 >  1:321818/1‑250 (MQ=255)
                  ttattaCTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGCTTTTTTTTGGCGTTT‑aacgctctgcag                                                                                                                                                                                                                                      >  2:382956/1‑239 (MQ=255)
                      taCTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATgcag                                                                                                                                                                                                                                  >  1:231651/1‑250 (MQ=255)
                        cTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGcc                                                                                                                                                                                                                                                       >  1:51979/1‑228 (MQ=255)
                         tCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCg                                                                                                                                                                                                                                                      <  2:51979/228‑1 (MQ=255)
                           aTACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATgcagcag                                                                                                                                                                                                                               <  2:883122/248‑1 (MQ=255)
                                         atatGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTTCGTTT‑TTCCTGCGCAATGCAGGAGCTGCGGTATGACA‑GGc                                                                                                                                                                                                              >  2:644773/1‑250 (MQ=255)
                                         atatGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTT‑TCCCTTCGCCATGACGCAGCCGCGGTCCAACA‑GGc                                                                                                                                                                                                              >  2:592011/1‑250 (MQ=255)
                                            tGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGccc                                                                                                                                                                                                            <  2:400581/249‑1 (MQ=255)
                                              gATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTT‑TCCCTTCGCAAATCATCAGCCGCGGTtcacccagcacccc                                                                                                                                                                                                          >  2:404436/1‑239 (MQ=255)
                                               aTTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGa                                                                                                                                                                                                        <  1:344765/250‑1 (MQ=255)
                                                   ataatTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTGCCGTTT‑TCCCTGCGCAATGCAGCCGCAGCGGGACAACA‑GGCCCGCGaaa                                                                                                                                                                                                      >  1:695606/1‑248 (MQ=255)
                                                   ataatTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTT‑TACCTGCGCCAAGCCGCAGCACCGGTACAACC‑AGCCcgcgcaca                                                                                                                                                                                                     >  2:723809/1‑245 (MQ=255)
                                                     aatTCTTATCCCAAAATGTGTAGACTCCCCCTATGGTAAGCGCGCTAAAGATTCAATATATAAACTTTAATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAGCGATGGGTGGGCATCAGCCGCAGAACACTCGTTAAAGCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAgg                                                                                                                                                                                                  <  2:528276/250‑2 (MQ=255)
                                                      atTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCCGGCGGTTTTTTTTTGCCGTTT‑ACTCTGCCCAATGCACCAAGCGCCGTGAAAAA‑GGCCCGCCCAAAAATg                                                                                                                                                                                                 >  1:797433/1‑250 (MQ=255)
                                                           tAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTg                                                                                                                                                                                            >  1:553759/1‑250 (MQ=255)
                                                                  aaaTGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGGACAACA‑AGCCCGCGAAAAAGAGGTCTGGGGGGcc                                                                                                                                                                                     >  1:414323/1‑250 (MQ=255)
                                                                         tAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCg                                                                                                                                                                              <  1:2901/250‑1 (MQ=255)
                                                                            aCTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTT‑ACCCTGCGCAATGCAGCAGCACCGGGACACCA‑AGCCCGCGCAAAAGAGGGTTGGGGGGCCCGTTcccgc                                                                                                                                                                            >  2:398054/1‑246 (MQ=255)
                                                                             cTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTT‑ACCCTGCGCAATGCAGCCGCAGCGGGACAACA‑GGCCCGCGGAAAAGTGGGCTGGGGGGCCTGTTCCGGCaa                                                                                                                                                                          >  1:770760/1‑250 (MQ=255)
                                                                                 ctctATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGGACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTAACCCTGCGCAATG‑AGCAGCAGCGGGTCAAAAGGGCCCGGGAAAAGGGGGTTTGGGGTGCCT‑TTCCggtaaccgt                                                                                                                                                                      >  2:787060/1‑242 (MQ=255)
                                                                                      ttGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGt                                                                                                                                                                   <  2:666608/248‑1 (MQ=255)
                                                                                            gcgcgcTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCtg                                                                                                                                                           <  1:787618/250‑1 (MQ=255)
                                                                                            gcgcgcTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCt                                                                                                                                                            <  2:684001/249‑1 (MQ=255)
                                                                                                   aaaTATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTT‑ACCCTACGCAATGCAGCAGCAGCGGGACAACA‑GGCCCGCGGAAAAGCGGGCTGGGGGGACTGTTCCGGCAACTGTGGGAGCCGCCTTGCACtt                                                                                                                                                    >  1:77358/1‑250 (MQ=255)
                                                                                                     atatTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCg                                                                                                                                                  <  1:827151/250‑1 (MQ=255)
                                                                                                     atatTCAATATATAAACTTTTATATAACGATAAAGAAAAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCg                                                                                                                                                  <  2:630818/250‑1 (MQ=255)
                                                                                                           aaTATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTT‑ACCCTGCGCAATTCAGCCGCCGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGGGCCCGTTCCCGCCACTGTGGTAGCCGCTGTGGCCTTCGTATAcc                                                                                                                                            >  2:883225/1‑249 (MQ=255)
                                                                                                               tataAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGt                                                                                                                                         <  1:486378/249‑1 (MQ=255)
                                                                                                                   aaCTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAg                                                                                                                                    <  2:543983/250‑1 (MQ=255)
                                                                                                                    aCTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGa                                                                                                                                   <  2:93494/250‑1 (MQ=255)
                                                                                                                     cTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGAt                                                                                                                                  <  1:33129/250‑1 (MQ=255)
                                                                                                                        ttATATAACGATTAAGAACAGGTAGTGAGTTAAGTCCAAAAATGAACGAATGGTATCCATCAGCTACCGAACACGCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCGGACGGTTTTTCTTAACCGTTT‑ACCCAGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGGCTGGGGTGCATGTTGCGTCAACTGTGGTAGCCGCGATGAACTTCGTCTACATGATAAAGATaa                                                                                                                                <  2:355604/249‑1 (MQ=255)
                                                                                                                         tatataACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTT‑TACCCGCGCAATGCAGCAGCCGCGGGACCACC‑GGCCCGCGAAAAAGGGGGCTTGGGGGGCTGTTCCGTCCACCGTGGGTGCCGCTGGGCACCTCGGCTACCTGTTTAAGATAAt                                                                                                                               >  2:221821/1‑249 (MQ=255)
                                                                                                                           tataACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGaa                                                                                                                            <  1:44503/250‑1 (MQ=255)
                                                                                                                              aaCGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTGCCGTTT‑AACCCGCGCAATGCAGCAGCAGCGGGACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGGGGCTGTTCCGGCAACTGTGGTAGcc                                                                                                                                                               >  2:477686/1‑212 (MQ=255)
                                                                                                                               aCGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCg                                                                                                                                                              <  1:477686/212‑1 (MQ=255)
                                                                                                                                        aaCAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTGCCGTTT‑TACCTGCGCAATGCAGCAGCAGGGGGACAAAA‑GGCCCGCGGAAAAGGGGTCTGGGGGGCCTGTTCCGGCAACTGTGGTAGCCGCCGTGCACATCGTATACCTGTTTAAAGTAATGAAGAGA‑GCAgggggg                                                                                                              >  2:307790/1‑250 (MQ=255)
                                                                                                                                           aGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTGCCGTTT‑ACCCTGTGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGGCTGGGGGGCCTGTTCCGTCAACTGTGGTAGGCGCTGTGCCCTTTGGCTACATGTTAAAGATAATGAAGTGA‑ACTggggggga                                                                                                            >  2:484564/1‑247 (MQ=255)
                                                                                                                                                   agtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTGCCGTTT‑ATCCTGCGCAATGCAGCAGCAGCGGGACAACA‑AGGCCCCGAAAAAGTGGTCTGGGGGGCCTGTTCCCTCCACTGTGGGAGGCGCTGTGGACTTCGTCTACATGTTAAAAAACATGTAGTGACCCTGGGGGGaacttccaa                                                                                                     >  2:682853/1‑242 (MQ=255)
                                                                                                                                                    gtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATa                                                                                                  <  1:777766/250‑1 (MQ=255)
                                                                                                                                                          tCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATACCGGCa                                                                                            <  1:428789/250‑1 (MQ=255)
                                                                                                                                                                       gAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTT‑TCCCCGCGCCATGCCGCAGCAGCGGTACAACA‑GGCCCGCGGAAAAGTGGTCTGGGGGGCCTGTTCCCTCAACCGTGGTGGCCCCTGTGCACTTTGGCTACATTTTTAAAATTATTGAGAGA‑CCCTGGGGGAAAATGACAAAACCGGCAGCGAttactaagg                                                                               >  2:249024/1‑242 (MQ=255)
                                                                                                                                                                          tggtggGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGTTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCaa                                                                            <  1:767451/250‑1 (MQ=255)
                                                                                                                                                                            gtggGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCCTTT‑AACCTGCGCCATGCAGCAGCAGCGGGACAACA‑AGCCCCCGGAAAAGAGGTGTGGGGGGCCCGTTCCGCCACCAGTGGGAGCCGCTGTGCACTTCGTCTACACGTTAAAGATAATGAAGTGA‑CCTGGGGGGGAACAGACACAAACCGCAGAGGCGAGTAAGGCAAc                                                                           >  2:587035/1‑249 (MQ=255)
                                                                                                                                                                                            gAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGACAACCATCA‑GGTACGCGCCTg                                                         <  1:255888/250‑1 (MQ=255)
                                                                                                                                                                                              acacTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGGAAAAGGGGTCTGGGGGGCCTTTTCCGTCAACTGTGGGAGCCGCTGTGCACCTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTTACAATACCGGCAGCGATGAGTACGGCAACCATCA‑GGTACGCGCCTGtt                                                       >  1:575652/1‑250 (MQ=255)
                                                                                                                                                                                                   cGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTT‑TCCCTGCGCCATGCAGCAGCACCGGTACCACA‑GGCCCGCGCAAAAAGGGGCCGGGGGGCCTGTTCCGTCAACTGTGGGAGCCGCTGTGCACCTCGTCTTCATGTTTAAGATTATGCAATGA‑CCTGGGGGGGAACCTGCAATACCGGCAGCGATGAGGAAGGCCACCCACAGGGAACGC‑CCTGTTTgccc                                                  >  2:77622/1‑248 (MQ=255)
                                                                                                                                                                                                    gTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCA‑GGTACGCGCCTGTTTgcgc                                                  <  2:354552/249‑1 (MQ=255)
                                                                                                                                                                                                       aaaTCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCA‑GGTACGCGCCTGTTTGCGCGGt                                               <  2:229494/249‑1 (MQ=255)
                                                                                                                                                                                                        aaTCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCA‑GGTACGCGCCTGTTTGCGCGGTc                                              <  2:609711/249‑1 (MQ=255)
                                                                                                                                                                                                                gcgATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCA‑GGTACGCGCCTGTTTGCGCGGTCGCTCCATcc                                     <  1:398054/250‑1 (MQ=255)
                                                                                                                                                                                                                      ggTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCA‑GGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCgg                                <  1:697281/249‑1 (MQ=255)
                                                                                                                                                                                                                             tGGCGCGGGCTGAAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCA‑GGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTAtt                          <  2:296044/248‑1 (MQ=255)
                                                                                                                                                                                                                                 gctggctgCAGGCGGTTTTTCTTTTCCGTTT‑ACCCTGCGCAATGCAGCAGCAGGGGGAAAACC‑GGCCCGCGAAAAAGGGGTCTGGGGGGGCTGTTCCGGCAACTGTGGGAGCCCCCGTGCCCCTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGGAACTGTACATATCGGCAGCGGAGAGGAGGGCAACCCTCA‑GGGACGCCGCCGGTTGCGCGGGCGGTCCATCCGCCCGCGGATTAatcat                    >  2:322079/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  tttCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCA‑GGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTGAAtt    <  1:14209/249‑1 (MQ=255)
                                                                                                                                                                                                                                                   ttCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCA‑GGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTGAATTAc  >  1:164596/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                           
ATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTT‑ACCCTGCGCAATGCAGCAGCAGCGGTACAACA‑GGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGA‑CCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCA‑GGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTGAATTAC  >  NZ_CP009273/1652175‑1652667

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