Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,492,5810AC35.3% 131.0 / 12.5 68L229R (CTG→CGG) insGIS4‑like element IS4 family transposase
Reads supporting (aligned to +/- strand):  ref base A (6/38);  new base C (24/0);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 3.95e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCCCCCATTTCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGAAGGCTGCATCGTTCTCTGGTGTATCC  >  NZ_CP009273/4492335‑4492829
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          
tCCCCCATTTCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGc                                                                                                                                                                                                                                                         <  1:175698/250‑1 (MQ=255)
 cccccATTTCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCt                                                                                                                                                                                                                                                        <  2:257833/250‑1 (MQ=255)
     cATTCCCCCCCCGGGGAAGCGCATGGCGTCCGTCATCGCCGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAg                                                                                                                                                                                                                                                    <  2:177531/250‑1 (MQ=255)
             ctccGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCAtt                                                                                                                                                                                                                                              <  2:708523/247‑1 (MQ=255)
               ccGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACa                                                                                                                                                                                                                                         <  1:257551/250‑1 (MQ=255)
                 ggggAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGt                                                                                                                                                                                                                                       <  1:227304/250‑1 (MQ=255)
                          cATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGGTTTTTCCCCCCCGGCTCCAGGGACTTTAACAGCCCAGTggg                                                                                                                                                                                                                              >  2:250931/1‑248 (MQ=255)
                           aTGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTTCCGCCCGGCCCCCGGCCCTTAACAGTcccccggggt                                                                                                                                                                                                                             >  2:359738/1‑243 (MQ=255)
                             ggCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCCCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCCGGCTCCAGGCATTTAACAGTCCCCGGGAGgaa                                                                                                                                                                                                                           >  1:23839/1‑250 (MQ=255)
                               cgtccgtcATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGtaata                                                                                                                                                                                                                         <  2:105714/250‑1 (MQ=255)
                                      cATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTa                                                                                                                                                                                                                  <  2:46883/250‑1 (MQ=255)
                                        tCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCCCCAGGCTCCCGGCATTTAACAGTCCCCGTGAGTAATAACCttttttc                                                                                                                                                                                                               >  1:107614/1‑246 (MQ=255)
                                                  cagATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGAGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTc                                                                                                                                                                                                      <  2:558424/250‑1 (MQ=255)
                                                   agATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCa                                                                                                                                                                                                     <  1:443987/250‑1 (MQ=255)
                                                   agATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCa                                                                                                                                                                                                     <  1:474448/250‑1 (MQ=255)
                                                           gACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTa                                                                                                                                                                                             <  2:293214/250‑1 (MQ=255)
                                                                 tCCTTTGCGCGTCCCGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGcc                                                                                                                                                                                        <  1:156913/249‑1 (MQ=255)
                                                                  ccTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTCCCCCCAGGCTCCCGGCATTTTACAGTCCCAGTGAGTGATAAACCTTATTCATTAACGTCAGAGTGTTATCGCCgg                                                                                                                                                                                      >  1:697749/1‑250 (MQ=255)
                                                                   cTTTGCGCGTCACGATCAGCAGGCGTGCAGTCTCTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCgg                                                                                                                                                                                      <  1:565391/249‑1 (MQ=255)
                                                                      agcgcGTCGCGTACACCAGGCGGGCAGGCACTTCATTTCCCAGTCCCGGCCCCTTTTATCGTGCCTGCGGGCTGGTGTTCAGCTTCACCAGATGATCGCCTTTACCCAGAGTTCTGATCTCTTCAAATTGCGCTCCCTTTCTGAAAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGtt                                                                                                                                                                                    <  2:612980/247‑1 (MQ=255)
                                                                            tCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCtat                                                                                                                                                                            <  2:359300/250‑1 (MQ=255)
                                                                              aCGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCtata                                                                                                                                                                           <  1:624776/249‑1 (MQ=255)
                                                                                     gcagGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTc                                                                                                                                                                  <  1:780560/250‑1 (MQ=255)
                                                                                       agGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGATTGTTCTATAAGTTGCTCAg                                                                                                                                                                <  1:774314/250‑1 (MQ=255)
                                                                                           gggCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCCCCCGGCCCCCGGCCTTTTACAGTCCCCGTGAGGAATAACCTTTATCCCTTAACGCCAGAGTGTTATTGCCCGGTTGTTCTTTTAGTTTCTCAGCAAg                                                                                                                                                            >  1:62947/1‑250 (MQ=255)
                                                                                               aGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTCCCGCCCCGCTCCCGGCCTTTTACAGGCCCAGTGAGTAATAACCTTTTTCCATTTACGTCAGAGTGTTATTGCCGGCTTGTTCTTTTAGTTGTTCAGCAAGCCCa                                                                                                                                                        >  2:125137/1‑250 (MQ=255)
                                                                                                gTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCAt                                                                                                                                                       <  1:305970/250‑1 (MQ=255)
                                                                                                    cTTCATTTCCCAGTCCCGGCCCCTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTc                                                                                                                                                   <  2:783113/250‑1 (MQ=255)
                                                                                                    cTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTc                                                                                                                                                   <  1:390155/250‑1 (MQ=255)
                                                                                                     ttCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTCCCGCCCGGCTCCCGGCCTTTTACAATCCCCGTGGGTAATAACCTTTTTCCATTTAACTCAGAGTGTTTTCTCCCGGTTTTTCTATAAATTTCTCCGCCAGCGCATTTTct                                                                                                                                                  >  2:152587/1‑249 (MQ=255)
                                                                                                       cATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCg                                                                                                                                                  <  2:23979/248‑1 (MQ=255)
                                                                                                       cATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCCCCCGGCTCCCGGCATTTTACAGTCCCCGTGAGGAATAAACTTTATTCCTTAAAGTCCGGGGGTTTTTGCCGGGTTGTTCTTTAAGTTGCTCAACCAGCTCAttttttct                                                                                                                                                >  2:296046/1‑246 (MQ=255)
                                                                                                        aTTTCCCGGTCCCGGCCACTTTTTTCGTGCTTGCGGGCTGGTTTTCAGCTTCCCCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTGTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCt                                                                                                                                                <  2:376391/249‑1 (MQ=255)
                                                                                                          ttCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCCGGCTCCCGGCATTTAACAGTCCCAGTGAGGAATAAACTTTATCCATTAACGTCAGAGGGTTTTCGCCGGTTTTTTCTTTAAGTTGCTCAGCAAGCCCCTTTTTGCTGtt                                                                                                                                             >  1:462229/1‑250 (MQ=255)
                                                                                                              cAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCCGGCTCCCGGCATTTAAAAGTCCCCGTGGGGAATAAACCTTTTCCATTAACGTCAGAGGGTTTTTGCCGGGTTTTTTTTTAAATTGCTCAGCAAACTCATTTTTGCTGttcttt                                                                                                                                         >  1:365905/1‑249 (MQ=255)
                                                                                                                gTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCCGGCTCCAGGCGTTTTACAGTCCCAGTGAGGAATAAACCTTTTCCATTTACGTCAGAGTGTTATCGCCGGTTTGTTCTTTAAATTGCTCAGCCAGCTCCTTTTTGCTGttcttcct                                                                                                                                       >  2:383812/1‑248 (MQ=255)
                                                                                                                          cTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAg                                                                                                                             <  2:174608/250‑1 (MQ=255)
                                                                                                                              tttCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGGTCTCCCCCCCAGCTCCAGGCATTTTACAAGCCCCGGGAGTAATAAACCTTATCCATTAACGTCAGAGGGGTTTTGCCGGGTTTTTTTTTAAATTGCTCCAGAAAGCCCTTTTTTCTGTTTTTCATCAGGCCgcaatgcg                                                                                                                         >  2:349403/1‑243 (MQ=255)
                                                                                                                               ttCTTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTg                                                                                                                         <  2:448599/249‑1 (MQ=255)
                                                                                                                                  gTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCCGGCTCCAGGCATTTAACAGTCCCCGTGAGGAATAAACTTTATCCATTAACGTCAGAGTGTTATCGCCGGGTTGTTCTTTAAGTTGCTCCGCAAGCTCCTTTTCGCTGTTTTTCCTCGTGCCGAAGGCTGccac                                                                                                                     >  2:310035/1‑247 (MQ=255)
                                                                                                                                    gCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTCCCGCCCCGCTCCAGGCCTTTTAAAGTCCCCGTGAGTAATAAACTTTATTCCTTAACGTCAGAGGTTTATCCACGGTTTGTTCTTTTAATTGCTCCGCAAACTCATTTTCGCTGTTCTTCATCGTTCCGAAAGCTGCAGCCg                                                                                                                   >  1:705875/1‑250 (MQ=255)
                                                                                                                                         cGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTcagc                                                                                                              <  1:383812/250‑1 (MQ=255)
                                                                                                                                                   ttCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTTCCGCCCGGCTCCAGGCATTTAACAATCCCCGTGAGTAATAAACCTTTTCCATTAACCTCAGAGGGGTATCGCCGGGTTGTTTTTTTAGTTTCTCAGCAAGCACATTTTCGCTGTTTTTCATTGTGCCCACGGCTGCAGCCGTCAGCAGATGGCTgg                                                                                                    >  2:57432/1‑250 (MQ=255)
                                                                                                                                                            aCCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCa                                                                                           <  2:111294/250‑1 (MQ=255)
                                                                                                                                                              cAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTCTTTTCCCGCCAGGCTCCCGGGAATTAAAAGTCCCAGTGAGGAATAAACTTTTTCCATTAAAATCAGAGGGTTATAGCCCGTTTTTTCTATTAACTGCTCCGCAAGATCATTTTTGCTGTTCTTCTTCGTGCCGAAAGCTGCAGCCCCCAGAAAAGTGCTGGTCAGTTTCacc                                                                                         >  2:573108/1‑248 (MQ=255)
                                                                                                                                                                        ccTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCCGGCTCCAGGCATTTTACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGGTTGTTCTTTAAATTGCTCAACAAGCCCATTTTTGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCCGTTTCATCTGGCAGAccc                                                                               >  2:543051/1‑249 (MQ=255)
                                                                                                                                                                            tACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCAtt                                                                             <  1:607295/248‑1 (MQ=255)
                                                                                                                                                                              cccAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTCCCCCCCGGCTCCAGGCATTTAACAGTCCCAGTGAGGAATAAACCTTATCCATTAACGTCAGAGGGTTATCGCCGGTTTGTTTTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCGGCCACCGCCAGCAGATGGCTGGTCAGTTTCCTCTGGCAGACCAttttta                                                                         >  2:642664/1‑248 (MQ=255)
                                                                                                                                                                              cccAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTg                                                                          <  2:687898/249‑1 (MQ=255)
                                                                                                                                                                                cAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGAcc                                                                       <  1:31363/250‑1 (MQ=255)
                                                                                                                                                                                 atgtttCAGATCTCTTCATATTGCGCTCCCGTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCt                                                                      <  2:740417/247‑1 (MQ=255)
                                                                                                                                                                                 aggtttCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCt                                                                      <  2:187428/247‑1 (MQ=255)
                                                                                                                                                                                     ttCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTCCCCCCCGGCTCCAGGCATTTAACAGTCCCCGTGAGTAATAACCTTTTTCCATTAAAGTCAGAGTGTTATCGCCGGGTTGTTTTATAAGTTGCTCAGCCAGCTCATTTTTGCTGTTTTTTATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCCTCTGGCAGACCATTTTTACCTGCgg                                                                  >  1:101685/1‑250 (MQ=255)
                                                                                                                                                                                         gATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTTCCGCCCGGCCCCCGGCATTTAACAGTCCCCGTGAGGAATAAACTTTTTCCATTAAAGTCAGAGTGGTATTGCCGGGTTGTTTTTTTAGTTTCTCAGGAAGCTCATTTTTTCTGTTCTTCATCGTGCCGCAAGCTGCAGCCGTCACCAGATGGCTGGGCCGGTCCATCTGGCAGCACCTTTTTACCTGCCGGTag                                                              >  2:568074/1‑250 (MQ=255)
                                                                                                                                                                                             tctTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTTCCCCCCCGCTCCCCGCATTTAAAAATCCCCGTGAGGAATAACCTTTTTTCCTTATCGTCCGAGGGTTATTGCCCGTTTTTTCTTTAAATTTTTCCGCCAGCTCATTTTTTCTGGTCTTCATTCTGCCGAAGGCCGACGCCGTCAGCAGTTGGCTGGTACGTTGCCTCCTGCAGGACATTTTTGCCTGTggggtggggg                                                          >  2:432629/1‑243 (MQ=255)
                                                                                                                                                                                             tctTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTCCCGCCCCGCTCCCGGCATTTAACAGTCCCAGTGAGGAATAACCTTTATCCCTTAACGTCCGAGTGTTATCGCCGGTTTGTTTTATAAGTTTCTCAGCAAACTCATTTTTGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCCTCTTGCCGACCATTTTTACCTTCGGGTGAAgcg                                                          >  2:396214/1‑250 (MQ=255)
                                                                                                                                                                                                    atTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCCGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAACTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCCGCAGATGGCTGGTCAGTTCCATCTGGCAGACCCTTTTGACCTGCGGGTAGAGCGCCGGGtt                                                   >  1:482988/1‑250 (MQ=255)
                                                                                                                                                                                                    atTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGGGTTTTTCCGCCCGGCTCCCGGCATTTTAAAGTCCCCGGGAGTAATAACCTTTATCCATTAACGTCAGAGGGTTTTCGCCGGGTTTTTCTTTAAGTTGGTCAGCAAGCTCATTTTTGCTGTTTTTCATCGTGCCGCAAGCTGCAGCCGTCAGCAGATGGGTGGGTAGTTTCATCTTGCAGACCCCTTTTACCTGCGGGGAGAGCGCCCGGtt                                                   >  2:614318/1‑250 (MQ=255)
                                                                                                                                                                                                     ttGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGtt                                                   <  1:731972/249‑1 (MQ=255)
                                                                                                                                                                                                       gcgcTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCCGGCTCCCGGCATTTAACAGTCCCAGTGGGGAATAACCTTTATCCATTAACGTCAGAGTGGTATCGCCGGTTTGTTTTTTAAGTTGCTCAGCAAGCCCATTTTTGCTGTTCTTCATCGTGCCGCAGGCTGCAGCCCTCCGCAGATGGCTGGTCAGTTTCCTCTGGCCGACCCTTTTGGCCCGCGGGGAGAGCGCCGGGTTccc                                                >  2:22287/1‑250 (MQ=255)
                                                                                                                                                                                                        cgcTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTCCCGCCAAGCTCCAGGCATTTAACAATCCCCGTGAGTAATAACCTTTATCCATTAACGTCCGAGTGTTATTGCCGGGTTGTTTTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGCAGGCTGCCGCCGTCCGCAGATGGCTGGTCAGTTCCCTCTGGCCGACCATTTTGGCCCTCGGGTAGAGCGCCGGGTTcccc                                               >  1:737470/1‑249 (MQ=255)
                                                                                                                                                                                                        cgcTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCg                                               <  1:311399/250‑1 (MQ=255)
                                                                                                                                                                                                        cgcTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCCCCAGGCTCCAGGGATTTTACACTCCCAGTGAGTAATAACCTTTATCCATTAAAGTCCGAGTGTTATCGCCGGTTTGTTCTTTAAGTTGCCCAGCAAGCTCATTTTTGCTGTTCTTTTTCCTGCCGAAGGCTGCAGCCGCCCGCAGAAGGCTGGTCAGTTTCATCTGGCACACCATTTTGACCCGCGGGGAGAGCGCCGGGTTCCCg                                               >  1:473789/1‑250 (MQ=255)
                                                                                                                                                                                                         gcTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCgg                                              <  2:19077/250‑1 (MQ=255)
                                                                                                                                                                                                                      gagGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCg                                 <  2:163295/250‑1 (MQ=255)
                                                                                                                                                                                                                             catcCAGTGGCGGTGTTCTCCCCCCCGGCCCCCGGCATTTTACAGTCCCCGTGAGTAATAACCTTTATCCCTTTACGTCAGAGTGTTATCGCCCGTTTTTTCTATAAGTTTCTCCGCAAGCTCCTTTTTTCTGTTCTTCCTCGTTCCGAAGGCTTCCGCCCTCCGCCGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTTACCTGCGGGTAGAGCGCCCGGGTCCCGGCCTGTGTCTGGCGGGGGGAg                          >  2:360077/1‑250 (MQ=255)
                                                                                                                                                                                                                                               tcCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGCAGGCTGCATCGGTCTCTggc        >  2:313858/1‑249 (MQ=255)
                                                                                                                                                                                                                                                     cAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGAAGGCTGCATCGTTCTCTGGTGTATcc  <  2:440033/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          
TCCCCCATTTCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGAAGGCTGCATCGTTCTCTGGTGTATCC  >  NZ_CP009273/4492335‑4492829

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