Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,240,6900GT33.3% 123.9 / 18.4 63L25F (TTG→TTT) BW25113_RS06225protein YcgY
Reads supporting (aligned to +/- strand):  ref base G (7/35);  new base T (21/0);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 4.29e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCGGCGTA  >  NZ_CP009273/1240445‑1240937
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        
cATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGATCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGTAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAACTGGGGCGGATAGTGAGGCGGATTTTGCTg                                                                                                                                                                                                                                                      >  1:197332/1‑250 (MQ=255)
 aTAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTgc                                                                                                                                                                                                                                                     <  1:446902/250‑1 (MQ=255)
           tAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGaa                                                                                                                                                                                                                                          <  1:662895/250‑1 (MQ=255)
                     gATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAAt                                                                                                                                                                                                                                <  2:395009/250‑1 (MQ=255)
                       tGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTg                                                                                                                                                                                                                              <  2:252853/250‑1 (MQ=255)
                          cGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGGTGG‑GCTTTTTTAAAATTTAGAATTTAt                                                                                                                                                                                                                            >  1:611927/1‑249 (MQ=255)
                               tCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATc                                                                                                                                                                                                                       <  1:473799/249‑1 (MQ=255)
                                     ttCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGC‑GCTTTTTGAACAATTTGAATTGTTCAAACCTTCttg                                                                                                                                                                                                                >  1:476488/1‑248 (MQ=255)
                                        acaTCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCa                                                                                                                                                                                                              <  2:262473/249‑1 (MQ=255)
                                         caTCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGa                                                                                                                                                                                                            <  2:209755/250‑1 (MQ=255)
                                             aGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGGTGG‑GGCTATTTAACAATTATAATTGATCAAACCtttttggcgatta                                                                                                                                                                                                         >  2:70209/1‑238 (MQ=255)
                                                 cGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGcc                                                                                                                                                                                                     >  1:202972/1‑249 (MQ=255)
                                                  gTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCaa                                                                                                                                                                                                   <  1:65870/250‑1 (MQ=255)
                                                         cGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGC‑GCTTTTTTTAAATTTTTAATTTATTAAACCTTTTGTCCGAATAAACCAAAATAacc                                                                                                                                                                                            >  1:610724/1‑250 (MQ=255)
                                                          gACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATAccc                                                                                                                                                                                            <  2:616622/249‑1 (MQ=255)
                                                              ccTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCg                                                                                                                                                                                       <  2:679614/250‑1 (MQ=255)
                                                                tGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGC‑GCTTATTTAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCgg                                                                                                                                                                                      >  1:339908/1‑249 (MQ=255)
                                                                 gCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGC‑GCTTATTTGACATTTAGAATTGATCAATCCTTCTGGCCGATTAAGCCAAAATAACCGTTTGGc                                                                                                                                                                                     >  1:46435/1‑249 (MQ=255)
                                                                 gCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGC‑GCTTTTTTAACATTTTGAATTTATCAATCCTTCTGGCCGATTAAACCAAAATAACCGTTCGGc                                                                                                                                                                                     >  1:621369/1‑249 (MQ=255)
                                                                  cTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGc                                                                                                                                                                                     <  1:559026/248‑1 (MQ=255)
                                                                         ggTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTTCTGG‑GCTTATTTAACATTTTGAAATTATCAATTCTTCTGGTCGATTATGTCGAAATACCCCTTTTGCAATGCtt                                                                                                                                                                              >  1:105635/1‑248 (MQ=255)
                                                                          gTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGGTTTTTTTAACTTTTGTAATT‑ATTCATTCTTCTTGCCGATTAAGACAAAATAAACACTTGGCAATGCTTGAg                                                                                                                                                                           >  2:388805/1‑250 (MQ=255)
                                                                             aaaaaaaTCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGC‑GCTTTTTGAACATTTAGAATTGATCAATCCTTTTGGCAGATTAAGACAAAATAACCGGTCGGCAATGCTTTAGTTa                                                                                                                                                                        >  1:663579/1‑250 (MQ=255)
                                                                                  aaTCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGC‑GCTTTTTTAACATTTAGAATTTATTAATCCTTCTGGCAGATTTAGCCAAAATACCCGTTCGGGAAAGCTTCAGTTAcgcg                                                                                                                                                                    >  1:312659/1‑249 (MQ=255)
                                                                                    tCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAag                                                                                                                                                                 <  2:520190/250‑1 (MQ=255)
                                                                                       atatAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAgg                                                                                                                                                              <  1:715885/250‑1 (MQ=255)
                                                                                          taGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAAt                                                                                                                                                           <  2:635983/250‑1 (MQ=255)
                                                                                            gCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATgg                                                                                                                                                         <  1:554493/250‑1 (MQ=255)
                                                                                            gCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATgg                                                                                                                                                         <  1:662479/250‑1 (MQ=255)
                                                                                                     tatCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTc                                                                                                                                                <  1:343802/250‑1 (MQ=255)
                                                                                                            tAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCttttttt                                                                                                                                         <  2:576588/250‑1 (MQ=255)
                                                                                                                     ttCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTaa                                                                                                                                <  2:55407/250‑1 (MQ=255)
                                                                                                                           gTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATtctc                                                                                                                          <  2:211196/250‑1 (MQ=255)
                                                                                                                              acaacaacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGG‑GCTTTTTTAACATTTTGAATTGATCAATACTTCTTGCAAATTAAGCCAAAAAACCCCTTCGGCAATGCTTTCTTTTCGCCGAGAGGAAAGGTTTGAATTCTTTTTTTTGTTGCTAAATTCTTTgt                                                                                                                       >  2:460484/1‑249 (MQ=255)
                                                                                                                                aacaacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGG‑GCTTTTTTAACATTTTGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGGCTGAATTCTTTTTTTTGTTTCTAAATTCTCTgggg                                                                                                                     >  2:137140/1‑248 (MQ=255)
                                                                                                                                   aacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTTTTGAACAATTATAAATGATACATCCTTCTGGCAGATTAAACCAAAATACCCGTTCGGCAAAGCGTCAGTTACGCGAAGAGGGATGGTCTGGATTCTTTTTTTTGTTGTTAAATTCTCTGGAgtgt                                                                                                                  >  1:744905/1‑250 (MQ=255)
                                                                                                                                         aTCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGC‑GCTTATTGAACATTTAAAATTGATCAATCCTTTTTGCAGATTTAGCCAAAATAACCGTTCGGCAATGCTTTCGTTTCGCGAAGAGGAATGGTCTTAATTCTTTTTTTTGGTTCTAAATTTTCTTGAGTGTCTGGAt                                                                                                            >  1:285845/1‑250 (MQ=255)
                                                                                                                                                   aaTTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGatt                                                                                                   <  1:5011/249‑1 (MQ=255)
                                                                                                                                                    attattGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGatta                                                                                                  <  1:420574/249‑1 (MQ=255)
                                                                                                                                                                   tttttGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAg                                                                                   <  1:665864/249‑1 (MQ=255)
                                                                                                                                                                      ttgttgTGATTCAGAAAATAAATATTGTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGaaa                                                                                <  2:331690/249‑1 (MQ=255)
                                                                                                                                                                        ggTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGTAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGc                                                                              <  1:431855/249‑1 (MQ=255)
                                                                                                                                                                          tgtgATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCaaa                                                                           <  1:61961/250‑1 (MQ=255)
                                                                                                                                                                              aTTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGC‑GCTTATTTAACATTTTGAATTTATCAATCCTTCTGGCAGATTAAGCCAAAAAAACCGTTCGGCAATGCTTTAGTTACGCGAAGAGGGATGGTCTGAATTCTTTTTTTTGTTGCTTAATTTTCTTGAGGGTCTGGGTTATGTGATTTTTTATTTTAATCCAGAAAGCaaaaaa                                                                        >  1:256112/1‑249 (MQ=255)
                                                                                                                                                                               ttCAGAAAATAAATAGTTTAGAGAGTAGTCAGTCAGGCGTTAAAAATGGGGCAGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAAGCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCGGGAGTGTCTGGATTATGTGAGTATTAATTTAACTCCAGAAAGCAAAAAAAc                                                                      <  2:455096/250‑1 (MQ=255)
                                                                                                                                                                                tCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGC‑GCTTATTTAACATTTTGAATTGATTAATCCTTTTGGCAGATTAAGCCAAAAAACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTTTTTTTTTTGTTGCTAAATTTTCTTGAGTGTCTGGGTTATGTGTTTTTTTATTTAACTCCAGAAAACAAAAAAACg                                                                     >  1:249599/1‑250 (MQ=255)
                                                                                                                                                                                tCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGC‑GCCTATTTAAAATTTAGAATTTGTCAATCCTTCTTGCAGATTAAGCCAAAAAAACCGTTCGGCCATGCTTCAGTTTCGCGAAGAGGAAAGGTTTGAATTCTTTTTTTTGGTGCTAAATTTTCTTGAGGTTTTGGGTTTTTTGGTTATTAATTTTAATCCAGAAAACaaaaaaaag                                                                     >  2:634603/1‑248 (MQ=255)
                                                                                                                                                                                  aGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTTAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAACGtt                                                                   <  1:1957/250‑1 (MQ=255)
                                                                                                                                                                                   gAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACAATTAAAAATGATCAATCCTTCTGGGAGATTAAGGCAAAATAACCGTTTGGCAATGCTTCAGTTAAGCGAAAAGGGATGGTCTGAATTTTTTTTTTGGTTGCTAAATTTTCTGGAGTGTCTGGGTTATGTGATTATTTATTTAAATCCAGAAAGCaaaaaaaaggtg                                                                  >  1:687112/1‑245 (MQ=255)
                                                                                                                                                                                         aaataTTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAAATGATCAATCCTTCTGGCAGATTAAACCAAAATAACCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGGCTGAATTTCTTTTTTTGTTGCTAAAATCTCTGGAGTGTCTGGGTTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAAGGTGGTGAGt                                                            >  1:68823/1‑250 (MQ=255)
                                                                                                                                                                                              ttttAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGC‑GCTTATTTTACATTTAGAATTGGTCAATCCTTCTTGCAGATTAAGCCCAAAAAACCGGTCGGCCATGCTTCCGTTACGCGGAGAGGAATGGTCTGGATTCTTTTTTTTGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTTACTCCAGAAAGCAAAAAAAAGTTGATGAGTGAaca                                                       >  1:319909/1‑250 (MQ=255)
                                                                                                                                                                                              ttttAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAACGTTGATGAGTGAac                                                        <  1:462799/249‑1 (MQ=255)
                                                                                                                                                                                               tttAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGC‑GCTTTTTTAACATTTAGAATTTATCAATCCTTCTGGCAGATTAAGCCAAAAAACCCGTTTGGCAATGCTTTAGTTAAGCGAAGAGGGATGGGCTGAATTCTTTTTTTGGTTGCTAAATTTTCTGGAGTGTGTGGGTTTTGTGATTTTTTATTTTAATCCAGAAAGCAAAAAAAAGTTGTTGAGTGAacac                                                      >  1:360985/1‑250 (MQ=255)
                                                                                                                                                                                                   gACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGC‑AAAAAACGTTGATGAGTGAACACCGTa                                                  <  1:605571/249‑1 (MQ=255)
                                                                                                                                                                                                         aGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTCTTTGCTGC‑GCTTATTGTAAAATTAGAATTGATCAATCCTTCTGGTAGGTTAAGACAAAATACCCGCTCGGGCATGGTTCTGTTTCGCGCAGAGGAATGGTTTGAATTTTTTTTTTTGTTGCTAAATTTTCTTGAGTGTCTGGATTTTGTGATTTTTATTTTTATTCCAGAAAGCCAAAAAAAGTTTATGAGTTAATACCCTAACAata                                            >  2:518398/1‑250 (MQ=255)
                                                                                                                                                                                                             aGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGC‑GCTTATTTAACATTTAGAATTTATCAATCCTTCTGGCAGATTAAACCCAAATACCCCTTCGGCAATGCTTCAGTTACGCGGAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTTAATTCTCTGGAGTGTCTGGGTTTTGTGATTATTTATTTTACTTCCGAAAGCAAAAAAACGTTGATGAGGTAAAACCGTAACCATAtt                                          >  1:123769/1‑248 (MQ=255)
                                                                                                                                                                                                                    cGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGC‑GCTTATTTAAAATTTTGAATTGATCAATCCTTTTTGCAGATTAAACCCAAATAACCCTTCGGCAATGCTTCAGTTTTGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTTTTTGGAGGGTCTGGATTTTGTGATTATTAATTTTAATCCAGAAAAAAAAAAAACGTTGATGAGTGAACACCGTAAAAATATTCAGGGAGc                                  >  1:619979/1‑249 (MQ=255)
                                                                                                                                                                                                                     gTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCaa                                <  1:354382/250‑1 (MQ=255)
                                                                                                                                                                                                                        aaaaaTGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAAtt                              <  1:272906/249‑1 (MQ=255)
                                                                                                                                                                                                                                ggCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCTTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTta                     <  2:476488/250‑1 (MQ=255)
                                                                                                                                                                                                                                 gCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTtat                    <  1:288342/250‑1 (MQ=255)
                                                                                                                                                                                                                                   ggATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTtatat                  <  2:611927/250‑1 (MQ=255)
                                                                                                                                                                                                                                                    tGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGGATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCGGCGTa  >  1:636347/1‑249 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        
CATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGC‑GCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCGGCGTA  >  NZ_CP009273/1240445‑1240937

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