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

CATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTGAATTACCCAAT  >  NZ_CP009273/1652181‑1652672
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                            
caTAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGgtt                                                                                                                                                                                                                                                        >  2:439013/1‑250 (MQ=255)
 aTAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTGCCGGTTa                                                                                                                                                                                                                                                      >  1:805101/1‑250 (MQ=255)
 aTAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGGTTa                                                                                                                                                                                                                                                      >  1:321818/1‑250 (MQ=255)
            ttattaCTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGCTTTTTTTTGGCGTTTaacgctctgcag                                                                                                                                                                                                                                           >  2:382956/1‑239 (MQ=255)
                taCTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATgcag                                                                                                                                                                                                                                       >  1:231651/1‑250 (MQ=255)
                  cTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGcc                                                                                                                                                                                                                                                            >  1:51979/1‑228 (MQ=255)
                   tCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCg                                                                                                                                                                                                                                                           <  2:51979/228‑1 (MQ=255)
                     aTACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATgcagcag                                                                                                                                                                                                                                    <  2:883122/248‑1 (MQ=255)
                                   atatGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTTCGTTTTTCCTGCGCAATGCAGGAGCTGCGGTATGACAGGc                                                                                                                                                                                                                   >  2:644773/1‑250 (MQ=255)
                                   atatGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTTCCCTTCGCCATGACGCAGCCGCGGTCCAACAGGc                                                                                                                                                                                                                   >  2:592011/1‑250 (MQ=255)
                                      tGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGccc                                                                                                                                                                                                                 <  2:400581/249‑1 (MQ=255)
                                        gATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTTCCCTTCGCAAATCATCAGCCGCGGTtcacccagcacccc                                                                                                                                                                                                               >  2:404436/1‑239 (MQ=255)
                                         aTTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGa                                                                                                                                                                                                             <  1:344765/250‑1 (MQ=255)
                                             ataatTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTGCCGTTTTCCCTGCGCAATGCAGCCGCAGCGGGACAACAGGCCCGCGaaa                                                                                                                                                                                                           >  1:695606/1‑248 (MQ=255)
                                             ataatTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTTACCTGCGCCAAGCCGCAGCACCGGTACAACCAGCCcgcgcaca                                                                                                                                                                                                          >  2:723809/1‑245 (MQ=255)
                                               aatTCTTATCCCAAAATGTGTAGACTCCCCCTATGGTAAGCGCGCTAAAGATTCAATATATAAACTTTAATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAGCGATGGGTGGGCATCAGCCGCAGAACACTCGTTAAAGCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAgg                                                                                                                                                                                                       <  2:528276/250‑2 (MQ=255)
                                                atTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCCGGCGGTTTTTTTTTGCCGTTTACTCTGCCCAATGCACCAAGCGCCGTGAAAAAGGCCCGCCCAAAAATg                                                                                                                                                                                                      >  1:797433/1‑250 (MQ=255)
                                                     tAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTg                                                                                                                                                                                                 >  1:553759/1‑250 (MQ=255)
                                                            aaaTGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGGACAACAAGCCCGCGAAAAAGAGGTCTGGGGGGcc                                                                                                                                                                                          >  1:414323/1‑250 (MQ=255)
                                                                   tAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCg                                                                                                                                                                                   <  1:2901/250‑1 (MQ=255)
                                                                      aCTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTACCCTGCGCAATGCAGCAGCACCGGGACACCAAGCCCGCGCAAAAGAGGGTTGGGGGGCCCGTTcccgc                                                                                                                                                                                 >  2:398054/1‑246 (MQ=255)
                                                                       cTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTACCCTGCGCAATGCAGCCGCAGCGGGACAACAGGCCCGCGGAAAAGTGGGCTGGGGGGCCTGTTCCGGCaa                                                                                                                                                                               >  1:770760/1‑250 (MQ=255)
                                                                           ctctATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGGACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTAACCCTGCGCAATGAGCAGCAGCGGGTCAAAAGGGCCCGGGAAAAGGGGGTTTGGGGTGCCTTTCCggtaaccgt                                                                                                                                                                           >  2:787060/1‑242 (MQ=255)
                                                                                ttGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGt                                                                                                                                                                        <  2:666608/248‑1 (MQ=255)
                                                                                      gcgcgcTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCtg                                                                                                                                                                <  1:787618/250‑1 (MQ=255)
                                                                                      gcgcgcTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCt                                                                                                                                                                 <  2:684001/249‑1 (MQ=255)
                                                                                             aaaTATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTACCCTACGCAATGCAGCAGCAGCGGGACAACAGGCCCGCGGAAAAGCGGGCTGGGGGGACTGTTCCGGCAACTGTGGGAGCCGCCTTGCACtt                                                                                                                                                         >  1:77358/1‑250 (MQ=255)
                                                                                               atatTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCg                                                                                                                                                       <  1:827151/250‑1 (MQ=255)
                                                                                               atatTCAATATATAAACTTTTATATAACGATAAAGAAAAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCg                                                                                                                                                       <  2:630818/250‑1 (MQ=255)
                                                                                                     aaTATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTACCCTGCGCAATTCAGCCGCCGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGGGCCCGTTCCCGCCACTGTGGTAGCCGCTGTGGCCTTCGTATAcc                                                                                                                                                 >  2:883225/1‑249 (MQ=255)
                                                                                                         tataAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGt                                                                                                                                              <  1:486378/249‑1 (MQ=255)
                                                                                                             aaCTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAg                                                                                                                                         <  2:543983/250‑1 (MQ=255)
                                                                                                              aCTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGa                                                                                                                                        <  2:93494/250‑1 (MQ=255)
                                                                                                               cTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGAt                                                                                                                                       <  1:33129/250‑1 (MQ=255)
                                                                                                                  ttATATAACGATTAAGAACAGGTAGTGAGTTAAGTCCAAAAATGAACGAATGGTATCCATCAGCTACCGAACACGCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCGGACGGTTTTTCTTAACCGTTTACCCAGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGGCTGGGGTGCATGTTGCGTCAACTGTGGTAGCCGCGATGAACTTCGTCTACATGATAAAGATaa                                                                                                                                     <  2:355604/249‑1 (MQ=255)
                                                                                                                   tatataACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTTACCCGCGCAATGCAGCAGCCGCGGGACCACCGGCCCGCGAAAAAGGGGGCTTGGGGGGCTGTTCCGTCCACCGTGGGTGCCGCTGGGCACCTCGGCTACCTGTTTAAGATAAt                                                                                                                                    >  2:221821/1‑249 (MQ=255)
                                                                                                                     tataACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGaa                                                                                                                                 <  1:44503/250‑1 (MQ=255)
                                                                                                                        aaCGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTGCCGTTTAACCCGCGCAATGCAGCAGCAGCGGGACAACAGGCCCGCGAAAAAGTGGTCTGGGGGGGCTGTTCCGGCAACTGTGGTAGcc                                                                                                                                                                    >  2:477686/1‑212 (MQ=255)
                                                                                                                         aCGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCg                                                                                                                                                                   <  1:477686/212‑1 (MQ=255)
                                                                                                                                  aaCAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTGCCGTTTTACCTGCGCAATGCAGCAGCAGGGGGACAAAAGGCCCGCGGAAAAGGGGTCTGGGGGGCCTGTTCCGGCAACTGTGGTAGCCGCCGTGCACATCGTATACCTGTTTAAAGTAATGAAGAGAGCAgggggg                                                                                                                   >  2:307790/1‑250 (MQ=255)
                                                                                                                                     aGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTGCCGTTTACCCTGTGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGGCTGGGGGGCCTGTTCCGTCAACTGTGGTAGGCGCTGTGCCCTTTGGCTACATGTTAAAGATAATGAAGTGAACTggggggga                                                                                                                 >  2:484564/1‑247 (MQ=255)
                                                                                                                                             agtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTGCCGTTTATCCTGCGCAATGCAGCAGCAGCGGGACAACAAGGCCCCGAAAAAGTGGTCTGGGGGGCCTGTTCCCTCCACTGTGGGAGGCGCTGTGGACTTCGTCTACATGTTAAAAAACATGTAGTGACCCTGGGGGGaacttccaa                                                                                                          >  2:682853/1‑242 (MQ=255)
                                                                                                                                              gtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATa                                                                                                       <  1:777766/250‑1 (MQ=255)
                                                                                                                                                    tCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCa                                                                                                 <  1:428789/250‑1 (MQ=255)
                                                                                                                                                                 gAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTTCCCCGCGCCATGCCGCAGCAGCGGTACAACAGGCCCGCGGAAAAGTGGTCTGGGGGGCCTGTTCCCTCAACCGTGGTGGCCCCTGTGCACTTTGGCTACATTTTTAAAATTATTGAGAGACCCTGGGGGAAAATGACAAAACCGGCAGCGAttactaagg                                                                                    >  2:249024/1‑242 (MQ=255)
                                                                                                                                                                    tggtggGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGTTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCaa                                                                                 <  1:767451/250‑1 (MQ=255)
                                                                                                                                                                      gtggGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCCTTTAACCTGCGCCATGCAGCAGCAGCGGGACAACAAGCCCCCGGAAAAGAGGTGTGGGGGGCCCGTTCCGCCACCAGTGGGAGCCGCTGTGCACTTCGTCTACACGTTAAAGATAATGAAGTGACCTGGGGGGGAACAGACACAAACCGCAGAGGCGAGTAAGGCAAc                                                                                >  2:587035/1‑249 (MQ=255)
                                                                                                                                                                                      gAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGACAACCATCAGGTACGCGCCTg                                                              <  1:255888/250‑1 (MQ=255)
                                                                                                                                                                                        acacTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGGAAAAGGGGTCTGGGGGGCCTTTTCCGTCAACTGTGGGAGCCGCTGTGCACCTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTTACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGtt                                                            >  1:575652/1‑250 (MQ=255)
                                                                                                                                                                                             cGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTTCCCTGCGCCATGCAGCAGCACCGGTACCACAGGCCCGCGCAAAAAGGGGCCGGGGGGCCTGTTCCGTCAACTGTGGGAGCCGCTGTGCACCTCGTCTTCATGTTTAAGATTATGCAATGACCTGGGGGGGAACCTGCAATACCGGCAGCGATGAGGAAGGCCACCCACAGGGAACGCCCTGTTTgccc                                                       >  2:77622/1‑248 (MQ=255)
                                                                                                                                                                                              gTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTgcgc                                                       <  2:354552/249‑1 (MQ=255)
                                                                                                                                                                                                 aaaTCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGt                                                    <  2:229494/249‑1 (MQ=255)
                                                                                                                                                                                                  aaTCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTc                                                   <  2:609711/249‑1 (MQ=255)
                                                                                                                                                                                                          gcgATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATcc                                          <  1:398054/250‑1 (MQ=255)
                                                                                                                                                                                                                ggTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCgg                                     <  1:697281/249‑1 (MQ=255)
                                                                                                                                                                                                                       tGGCGCGGGCTGAAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTAtt                               <  2:296044/248‑1 (MQ=255)
                                                                                                                                                                                                                           gctggctgCAGGCGGTTTTTCTTTTCCGTTTACCCTGCGCAATGCAGCAGCAGGGGGAAAACCGGCCCGCGAAAAAGGGGTCTGGGGGGGCTGTTCCGGCAACTGTGGGAGCCCCCGTGCCCCTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGGAACTGTACATATCGGCAGCGGAGAGGAGGGCAACCCTCAGGGACGCCGCCGGTTGCGCGGGCGGTCCATCCGCCCGCGGATTAatcat                         >  2:322079/1‑250 (MQ=255)
                                                                                                                                                                                                                                            tttCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTGAAtt         <  1:14209/249‑1 (MQ=255)
                                                                                                                                                                                                                                             ttCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTGAATTAc       >  1:164596/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   gCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTGAATTACCCAAt  >  1:498140/1‑249 (MQ=255)
                                                                                                                                                                                                                                                   gCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTGAATTACCCAAt  >  1:212684/1‑249 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                            
CATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTGAATTACCCAAT  >  NZ_CP009273/1652181‑1652672

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: