Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,652,4240GT33.3% 114.3 / 16.6 72L33F (TTG→TTTynfEDmsA/YnfE/YnfF family dimethyl sulfoxide reductase
Reads supporting (aligned to +/- strand):  ref base G (19/29);  new base T (24/0);  total (43/29)
Fisher's exact test for biased strand distribution p-value = 1.08e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTGAATTACCCA  >  NZ_CP009273/1652176‑1652670
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        
tACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGAGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGTCTGCAGGCGGTTTTTTTTTGc                                                                                                                                                                                                                                                         >  2:435390/1‑250 (MQ=255)
          tCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATGTTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCt                                                                                                                                                                                                                                               <  1:382341/250‑1 (MQ=255)
                   attaCTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTACCCTTCGCAATg                                                                                                                                                                                                                                       >  1:356164/1‑249 (MQ=255)
                   attaCTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTAACCTGCGCCATg                                                                                                                                                                                                                                       >  2:593509/1‑249 (MQ=255)
                              ccTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTACCCTGCGCAATGCAGCAGCAGCgg                                                                                                                                                                                                                           >  1:222288/1‑250 (MQ=255)
                              ccTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTCTTTTTTTGCCGTTTACCCTTCGCCATGCAGCAGACGCgg                                                                                                                                                                                                                           >  2:68772/1‑250 (MQ=255)
                                  cTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTac                                                                                                                                                                                                                        <  2:173710/249‑1 (MQ=255)
                                  cTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTa                                                                                                                                                                                                                         <  1:724181/248‑1 (MQ=255)
                                     ttGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTATTTTCCGTTTTACCCGCGCCATgcagcagcacgcgtccaac                                                                                                                                                                                                                     >  2:402936/1‑239 (MQ=255)
                                           tGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCcg                                                                                                                                                                                                              <  1:564969/250‑1 (MQ=255)
                                                taataaTTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTTCCCTTCGCAATGCAGCAGCCGCGGTGCAACAAGCCcgcgga                                                                                                                                                                                                         >  2:196166/1‑247 (MQ=255)
                                                        cTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTACCCTGCGCAATGCCACAGCCGCGGGAAAACCGGCCCGCGCAAAAgggggc                                                                                                                                                                                                >  2:651258/1‑248 (MQ=255)
                                                             cccAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTGCCGTTTACCCTGCGCAATTCAGCCGCAGCGGGACAACAGGCCCGCCGaaaaaagggcttggg                                                                                                                                                                                            >  1:410685/1‑240 (MQ=255)
                                                                       gTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTACCCTGCGCAATTCCGCAGCAGCGGGACAACAGGCCCGCGGAAAAGGGGTCTGGGGGGCCTGTTcc                                                                                                                                                                                 >  1:701667/1‑250 (MQ=255)
                                                                        tAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTGCCGTTTTCCCTGCGCAATGCAGCAGGCGCGGTACCACAGGCCCGCGGAAAAGTGGGCTGGGGGGCCTGTTTCg                                                                                                                                                                                >  1:3057/1‑250 (MQ=255)
                                                                                     ttGTTAGCGCGCTAAATATTCAATATATAAACTTTAATAAAACGATAAAGAACGGGAGTTGAGTTATGTCCAAAAAGGAAGGATTGGTGGGCATCGGCCGCAGACCACTCTTTAAATCTCCCGGGATGGGTTCTCGGGCGCTGGCGGCGGGCGGTTTTTCTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAg                                                                                                                                                                   <  2:603112/249‑1 (MQ=255)
                                                                                             gcgcTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTTACCTGCGCAATGCAGCAGCAGCGGGACAACAGGCCCGCGAAAAAGTGGTCTGGGGGGCCTGTTTCGTCAAATGTGGTAGCCGCtgtg                                                                                                                                                           >  1:717572/1‑250 (MQ=255)
                                                                                             gcgcTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCtgtg                                                                                                                                                           <  1:448961/250‑1 (MQ=255)
                                                                                               gcTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGGACAACAGGCCCGCGGCAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCCGGGCa                                                                                                                                                         >  1:153883/1‑250 (MQ=255)
                                                                                                 tAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTTCCCTGCGCACAGCAGCACCACGGGCACAACCGGCCCGCGCCAAAGAGGTCTGGGGGGCCTGTTCCTGCAACTTGGGTATCCGCCGTGccca                                                                                                                                                       >  2:612614/1‑247 (MQ=255)
                                                                                                     tatTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTTCCCTGCGCAAAGACGCAGCAGCGGGACACAAGGCCCGCGGAAAAGGGGGGTGGGGGGCCCGTTCCGTCCACTGTGGTAGCCGCTGTGGCCTTCGt                                                                                                                                                   >  2:409433/1‑250 (MQ=255)
                                                                                                       ttCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCCGTTACCCTGCGCAATGCAGCAGCAGCGGGATAAAAGGCCCGCGGAAAAGAGGGCCGGGGTGCCTGTTTCGTCCACCGGGGTAGCCGCCGGGGACTTCGCCt                                                                                                                                                 >  2:389887/1‑250 (MQ=255)
                                                                                                       ttCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGGGCCTGTTCCGTCAACCGTGGGAGCCGCTGTGCACTTCGTCt                                                                                                                                                 >  1:259437/1‑250 (MQ=255)
                                                                                                         cAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTa                                                                                                                                                <  2:239250/249‑1 (MQ=255)
                                                                                                          aaTATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCCTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGCAAAAGAGGTCTGGGGGGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGGCTACa                                                                                                                                              >  2:12849/1‑250 (MQ=255)
                                                                                                            tatataAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATg                                                                                                                                            <  2:251341/250‑1 (MQ=255)
                                                                                                             atataAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTACCCTGCGCAATGGAGCAGCCGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGGCTGTTCCCTCAACCGTGGGAGACGCTGTTCACGTCGTCTACATGt                                                                                                                                           >  2:241176/1‑250 (MQ=255)
                                                                                                                 aaaCTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTTCCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTa                                                                                                                                         >  1:221754/1‑248 (MQ=255)
                                                                                                                  aaCTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTaaa                                                                                                                                       <  1:257196/249‑1 (MQ=255)
                                                                                                                  aaCTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAg                                                                                                                                      <  1:494461/250‑1 (MQ=255)
                                                                                                                        tatataACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAAt                                                                                                                                 <  2:104328/249‑1 (MQ=255)
                                                                                                                           ataACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGGAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTGGCCGCTGTGGCATTCGTCTACATGTTAAAGATTATGAAg                                                                                                                             >  1:717148/1‑250 (MQ=255)
                                                                                                                                 aTAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGCTTTCCCTTCGCCCTGCCGCCACAGCGGGACAACAGGCCCCCGGAAAAGTGGGCTGGGGTGCCTGTTCCGTCAACTGTGGTTGCCGCTGTGCACTTCGGCTACATGTTTAAGGTAATGAAGTGACCt                                                                                                                       >  2:297918/1‑250 (MQ=255)
                                                                                                                                    aaGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACGGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTggg                                                                                                                    <  2:460440/250‑1 (MQ=255)
                                                                                                                                     aGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGt                                                                                                                   <  1:536526/250‑1 (MQ=255)
                                                                                                                                      gAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTg                                                                                                                  <  2:200710/250‑1 (MQ=255)
                                                                                                                                         cAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGGTTACCCTGCGCAATGCAGCAGCAGCGGGACAACAGGCCCGCGAAAAAGTGGGCTGGGGTGCCTGTTCCGTCCACTGTGGGAGCCGCTGTGCACTTCGTCCACATGGTAAAGATAATGAAGTGACCTGGGTggg                                                                                                                >  2:166197/1‑248 (MQ=255)
                                                                                                                                          aGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGGGGTCTGGGGGGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCCCTTCGTCTACATGTTTAAGATAATGAAGTGACCTGGGGGGaaa                                                                                                              >  1:654555/1‑250 (MQ=255)
                                                                                                                                                 tgagtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGCAAAAGTGGTCTGGGGGGCCCGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGAc                                                                                                        >  1:613127/1‑249 (MQ=255)
                                                                                                                                                 tgagtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTTCCCTGCGCAATGCAGCCGCAGCCGGACAAAACGCCCGCGGAAAAGTGGTTTGGGGGGCCCGTTTCCTCAACTGTGGGTACCCCTGTGCACCTTCTCTCCCAGTTTAAAGTAATGAAGTGACCTGGGGGGgaaacgac                                                                                                         >  2:249823/1‑244 (MQ=255)
                                                                                                                                                  gagtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTACCCTGCGCAATGCAGCAGCAGCGGGACAAAAAGCCCCCGAAAAAGAGGTCTGGGGGGCCCGGTCCGTCAACTGTGGGAGACGCCGTGAACCTCGTCTTCATGTTAAAGATAATGAAGTGACACGGGGGGAAACTGACa                                                                                                       >  2:164052/1‑249 (MQ=255)
                                                                                                                                                   agtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACaa                                                                                                      <  1:547818/249‑1 (MQ=255)
                                                                                                                                                    gtTAGGTCAAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATa                                                                                                    <  2:379647/250‑1 (MQ=255)
                                                                                                                                                          tCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTTCCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGGGCCTGTTTCCTCAACTGTGGTAGCCGCTGGGCACTTCGTCTACCTGTTAAAGATAATGAAGGGACCTGGGGGGGAACTGACAAAACCgg                                                                                                >  1:238835/1‑248 (MQ=255)
                                                                                                                                                            cAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGc                                                                                            <  1:588746/250‑1 (MQ=255)
                                                                                                                                                            cAAAAATGAACGAAGGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAg                                                                                             <  2:434675/249‑1 (MQ=255)
                                                                                                                                                             aaaaaTGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTACCCTGCGCAATGCAGCAGCAGCGGGACAACAGGCCCGCGAAAAAGTGGTCTGGGGGGCCTGTTCCGTCAACTGTGGGAGCCCCCGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGGCCTGGGGGGGAACTGACAAAACCGGCAACg                                                                                           >  2:711912/1‑250 (MQ=255)
                                                                                                                                                                aaTGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGGACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGAt                                                                                         >  1:179784/1‑249 (MQ=255)
                                                                                                                                                                 aggagagAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCAGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGGCTGGGGTGTGTGTTCCGTCAACTGTGGTAGCGGCTGTGCACTTCGTCTGCATGTTAAAGATAAGGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGGGATGa                                                                                       <  2:197397/244‑1 (MQ=255)
                                                                                                                                                                     aCGAATGGTGGGCATCAGCCGCAGAACGCTCGTTAAATCTGCCGCGATAGGTTCTCTGGGGCTGGCTGCAGGCGGGTTTTCTTGGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGGCTGGGGTGCCGGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTAc                                                                                   <  2:130967/250‑1 (MQ=255)
                                                                                                                                                                        aaTGGTTGACATCGGCGGCATAACACGCGTTAAATCTACCGCGATAGGTGCTCTGGCGCTTGCGGCAGGCGGGTTTTCTTTGCCGTTTACCCTGCGCAATGCATCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGCCAATACCGGCAGCGATGAGTACGGc                                                                                <  2:606950/250‑1 (MQ=255)
                                                                                                                                                                         aTGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAAGCCCGCGAAAAAGGGGTCTGGGGGGCCTGTTCCGTCCACTGTGGTAGCCGCTGTGCACCTCGTCTACATGTTAAAGATAATGGAGTGGACTGGGGGGGAACTGACAAAACCGGCAGCGGTGAGTACGGc                                                                                >  2:681013/1‑249 (MQ=255)
                                                                                                                                                                         aTGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGCTTACCCCGCGCAATGCAGCCGCCGCGGTACAACAGGCCCCCGGAAAAGCGGTCCGGGGGGCCCGTTCCGTCAACCGTGGGAGCCGCCCTGCCCCTCGTCTCCCTGTTAAAGATAATGAAGTGACCCGGGTGGGAACTGACCATACCCGCCGCGGTGAGGACGGcc                                                                               >  2:205760/1‑249 (MQ=255)
                                                                                                                                                                                  aTCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTACCCTGCGCAATGCAGCAGCAGCGGGACAAAAGGCCCGCGAAAAAGTGGTCTGGGGGGCCTGTTCCGTCAACTGTGGGAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAgg                                                                      >  1:394567/1‑250 (MQ=255)
                                                                                                                                                                                      gCCGCAGAACACGCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTAcg                                                                  <  1:457842/250‑1 (MQ=255)
                                                                                                                                                                                                 cTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTGCCGTTTACCCTGCGCAATTCAGCAGCAGCGGTACAACAGGGCCGCGGAAAAGTGGGCTGGGGGGCCTGTTTCGTCAACTGTGGTAGCCGCCGTGCACCTCGTCTACCTGTTAAAGATAATGAAGGGACCTGGGGGGAAAATGACCATACCGGCCGCGATGAGGACGGCCACCACTAGGGACGCGCCCGttttc                                                       >  2:360827/1‑248 (MQ=255)
                                                                                                                                                                                                   cGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTAACCTGCGCCATGCAGCAGCAGCGGGACAAAAGGCCCGCGGAAAAGTGGTCTGGGGGGCCCGTTCCGTCAACTGGGGGAGCCGCTGTGCAATTCGGCTACATGTTAAAAATAATGAAAGGGACTGGGGGGGAAAAGACAAAACCGGCAGCGGTGAGTACGGCAACCATCAGGGAAGCGCCCGTTTgcgc                                                     >  2:536051/1‑250 (MQ=255)
                                                                                                                                                                                                    gTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTgcgc                                                     <  1:313957/249‑1 (MQ=255)
                                                                                                                                                                                                       aaaTCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTTCGCTTTCCCTTCGCAATGCAGCAGCAGCGGGACAACAGGCCCGCGGAAAAGTGGGCCGGGGTGCCTGGTCCGTCAACTGTGGTAGCCGCTGTGGACCTCCTCTTCCTTTTTAAGATAATGGAGGGGCCTGGGGGGGAACTGGCCAAACCGGCAGCGACTAGGACCGGCACCACCACGGACGCCCCCGGTTGCGCgggc                                                 >  2:726581/1‑248 (MQ=255)
                                                                                                                                                                                                       aaaTCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTACCCTGCGCAATTCAGCAGCAGCGGGACAACAGGCCCGCGGAAAAAAGGGCTGGGGGGCCTGTTCCGTCAAATGTGGTAGCCGCTGTGCCCCTCGTCTACCTGTTTAAGATAATGAAGTGACCTGGGGGGGAACTGAACAAACCCGGAGGCGATGGGACGGCAACCAACAGGGACCGGCCTGTTTTGGCgggc                                                 >  2:484437/1‑248 (MQ=255)
                                                                                                                                                                                                       aaaTCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGGGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGGACGGCAACCATCAGGGACGCGCCTGTTTTCGCGGTc                                                 >  1:64947/1‑250 (MQ=255)
                                                                                                                                                                                                        aaTCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCCATGCAGCAGCAGCGGGACAACAGGCCCGCGAAAAAGGGGGCTGGGGGGCCCGTTCCGTCAAACGTGGGAGCCGCTGTGCACTTTGTCTAAATGTTAAAGATAAAGAAGTGACCTGGGTGGGAACTGACAATAACGGCAGCGATGGGTACGGCAACCATCAGGGACGCGGCTGTTTGCGCGGGCg                                                >  2:422350/1‑250 (MQ=255)
                                                                                                                                                                                                                  gATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATccgc                                      <  1:645003/250‑1 (MQ=255)
                                                                                                                                                                                                                       gTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTACCCTGCGCAATGCAGCAGCAGCGGGACAACAGGCCCGCGAAAAAGTGGGCTGGGGGGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGGGGGAACTGACAATACCGGCAGCGGTGAGTAAGGCAAACATCAGGGACGCGCCTGTTTGCGCGGGCGCCCCCTCCGCCGGCg                                 >  2:273842/1‑250 (MQ=255)
                                                                                                                                                                                                                               gcgcTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTaa                           <  1:121212/248‑1 (MQ=255)
                                                                                                                                                                                                                                cgcTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAatca                        <  2:356164/250‑1 (MQ=255)
                                                                                                                                                                                                                                         cAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATc                <  1:523850/249‑1 (MQ=255)
                                                                                                                                                                                                                                         cAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATc                <  1:143661/249‑1 (MQ=255)
                                                                                                                                                                                                                                             cGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCtt            >  1:456474/1‑249 (MQ=255)
                                                                                                                                                                                                                                               gTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACAGGCAGCGATGAGTACGGCAACCATCAGGGACGCGCCTGTTTGCGAGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTGaa         >  2:518678/1‑250 (MQ=255)
                                                                                                                                                                                                                                                tttttCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCAACCGCCGGCGTATTAATCATCCCGATCGCTTGAAt        <  1:703320/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 ttttCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTGAAt        >  2:254280/1‑249 (MQ=255)
                                                                                                                                                                                                                                                        ggCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCACCCCGATCGCTTGAATTACCCa  <  1:258953/247‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        
TACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTGAATTACCCA  >  NZ_CP009273/1652176‑1652670

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