Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,217,6890GA43.1% 65.9 / 18.3 51pseudogene (2228/2243 nt)BW25113_RS25990autotransporter outer membrane beta‑barrel domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (18/11);  new base A (0/22);  total (18/33)
Fisher's exact test for biased strand distribution p-value = 2.35e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGGAGAAAAACAGCGGTCCGCACGGTGACGACTTTGAACCTTTTGTTGAAATGAACTGGATCCATAACAGTAAAGATTTTGCTGTCTCAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGCATTACCCATTACCG  >  NZ_CP009273/1217443‑1217927
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                              
aGGAGAAAAACAGCGGGCCGCACGGTGACGACTTTGAACCTTTTGTTGAAATGAACTGGATCCATAACAGTAAAGATTTTGCTGGCTCAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGaaa                                                                                                                                                                                                                                             <  2:495037/250‑1 (MQ=255)
   agaAAAACAGCCGTCCGCACGGTGACGACTTTGAACCTTTTGTTGAAATGAACTGGATCCATAACAGTAAAGATTTTGCTGTCTCAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAatat                                                                                                                                                                                                                                          <  1:182201/250‑1 (MQ=255)
    gaAAAACAGCGGTCCGCACGGTGACGACTTTGAACCTTTTGTTGAAATGAACTGGATCCATAACAGTAAAGATTTTGCTGTCTCAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAtata                                                                                                                                                                                                                                         <  2:410741/250‑1 (MQ=255)
    gaAAAACAGCGGTCCGCACGGTGACGACTTTGAACCTTTTGTTGAAATGAACTGGATCCATAACAGTAAAGATTTTGCTGTCTCAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAtata                                                                                                                                                                                                                                         <  2:609543/250‑1 (MQ=255)
       aaaCAGCGGTCCGCACGGTGACGACTTTGAACCTTTTGTTGAAATGAACTGGATCCATAACAGTAAAGATTTTGCTGTCTCAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGt                                                                                                                                                                                                                                      <  2:417723/250‑1 (MQ=255)
          cAGCGGTCCGCACGGTGACGACTTTGAACCTTTTGTTGAAATGAACTGGATCCATAACAGTAAAGATTTTGCTGTCTCAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCt                                                                                                                                                                                                                                   <  2:31017/250‑1 (MQ=255)
              ggTCCGCACGGTGACGACTTTGAACCTTTTGTTGAAATGAACTGGATCCATAACAGTAAAGATTTTGCTGTCTCAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATc                                                                                                                                                                                                                               <  1:433689/250‑1 (MQ=255)
                   gCACGGTGACGACTTTGAACCTTTTGTTGAAATGAACTGGATCCATAACAGTAAAGATTTTGCTGTCTCAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCccgcc                                                                                                                                                                                                                          <  1:246628/250‑1 (MQ=255)
                             gacTTTGAACCTTTTGTTGAAATGAACTGGATCCATAACAGTAAAGATTTTGCTGTCTCAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTaa                                                                                                                                                                                                                 >  2:185423/1‑249 (MQ=255)
                                        ttttGTTGAAATGAACTGGATCCATAACAGTAAAGATTTTGCTGTCTCAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcg                                                                                                                                                                                                     >  2:629504/1‑250 (MQ=255)
                                           tgttgAAATGAACTGGATCCAAAAAAGTAAAGATTTTGCTGTCTCAATGAAGGGTGTAAAAGTCAAAAAAGAGGGGGCGAGTAATTGGGGGGAAATAAAACTGGGCGTAAATGGCAACCTGAACCCAGCGGCCAGCGTCTGGGGAATGTGGGGGTGCAGCTGGGTGATAATGGCTACAATGCCACCGCAGTGATGGGGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTa                                                                                                                                                                                                   <  1:374717/249‑1 (MQ=255)
                                              tgAAATGAACTGGATCCATAACAGTAAAGATTTTGCTGTCTCAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCa                                                                                                                                                                                               >  1:435009/1‑250 (MQ=255)
                                              taaaaTGAACTGAACCCATAACAGAAAAGTTTTTGCTGTCACAATGAACGGTGTAAAAGTCGAAAAAGAGGGGGCGAGTAATTGGGGGAAAATTAAACTTGGCGTAAATGGCACCCTGAACCCAGCGGCCAGCGTCGGGGGAATGTGGGCGTGCAGCGGGGTGATAATGGCTACAATGACCCCGCAGTGATGGTGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc                                                                                                                                                                                                <  1:169513/247‑1 (MQ=255)
                                                aaaTGAACGGGACCCATAAAAGAAAAGATTTTGCTGTCTCAATAAAGGGTGTGAAAGTCAAACAAGATGGGGCGAGTAATGTGGGGGAAATAAAACTGGGCGTAAAGGCCACCCTGAACCCGGCGCCCAGCGTCTGGGGAATGTGGGCGTGCAGCGGGGTAATAATGGCTACAATCCCCCCTAGGTGATGGTGGGCCTAAAATATAAGTCCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                                                                                                                             <  2:213779/250‑1 (MQ=255)
                                                   tGAACTGGATCCATAACAGTAAAGATTTTGCTGTCTCAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAc                                                                                                                                                                                          >  1:560781/1‑250 (MQ=255)
                                                          tatacATAACAAAAAAGTTTGTGCGGTCTAAAGGAGGGGTGGAAAAAGAGAAGAAGGGGGGGGGGGTAATTGGGGGAAAATAAAATTTGGGAAAAAGGGCACCCGAAACCGGGCGGCCGGCGTGTGGGGAAGGGGGGGGGGGGGCGGGGGAATAATGGCTACAATGACCCCGCGGTGAGGGTGGGCGTGAAATATAAGTTGTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAc                                                                                                                                                                                   <  2:223700/246‑1 (MQ=255)
                                                           acccAAAACAGAAAAGATTTTGCTGTCTAAATGAACGGTGTAAAAGTCGAAAAAGAGGGGGCGAGTAATTGGGGGAAAATTAAACTGGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCGGGGGAATGGGGGCGTGCAGCTGGGTGATAATGGCTACAATGACCCCGCAGTGATGGGGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACa                                                                                                                                                                                  <  1:57105/248‑1 (MQ=255)
                                                                   cAGTAAAGATTTTGCTGTCTCAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaaa                                                                                                                                                                           >  1:127577/1‑249 (MQ=255)
                                                                          gATTTTGCTGTCTCAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttt                                                                                                                                                                   >  1:316815/1‑250 (MQ=255)
                                                                             tttGCTGTCAAAATGAAGGGTGTAAAAGTCAAACAAGAGGGGGCGAGTAATTGGGGGAAAATTAATCTTGGCGTAAATGGCACCCTAAATCCAGCGGCCAGCGTCTGGGGAATGGGGGGGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttttt                                                                                                                                                                 <  1:442323/249‑1 (MQ=255)
                                                                              ttGCTGTCAAAATAAGCGGTGTAAAAGTCAAAAAAGAGGGGGCGAGTAATTGGGGGAAAATTAAACTTGGCGAAAAGGGCAACCTGAACCCGGCGGCCGGCGTCGGGGGATTGGGGGCGTGCAGTTGGGTGATAATGGCTACAATGACCCCGCAGTGAGGGTGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttttt                                                                                                                                                                <  1:173662/249‑1 (MQ=255)
                                                                               tGCTGTCTCAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGc                                                                                                                                                              >  1:574588/1‑250 (MQ=255)
                                                                                gttgcataaaggAACGGGGTAAAAGACGAAAAAGGGGGGGGGAGTAATTGGGGGAAAATTAACCTGGGGGAAAAGGGCAACCTGAACCCGGCGGCCAGCGTGTGGGGAGGGTGGGGGGGCAGCTGGGTAATAATGGCAAAAATGACCCCGAAGTGATGGGGGGCTTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGc                                                                                                                                                              <  2:21896/238‑1 (MQ=255)
                                                                                      tcAATGAAGGGTGTAAAGGAGAAAAAAGAGGGGGCGAGTAATTGGGGGAAAATAAAACTGGGCGAAAATGGCAACCTGAACCCAGCGGCCGGGGTCGGGGGAATGGGGGCGTGCAGTTGGGTAATAATGGCAAAAATGCCCCCGAAGTGATGGGGGGCTTAAAATATATGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAg                                                                                                                                                        <  2:628740/249‑1 (MQ=255)
                                                                                            aaCGGTGTAAAAGTCAAACAAGGGGGGGCGAGTAATTTGGGGAAAATTAAACTTGGCGAAAATGGCACCCTAAACCCAGCGGCCGGCGTGGGGGAAATGTGGGGGGGCAGCGGGGTGATAATGGCTACAATGACCCCGCAGTGAGGGTGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCAt                                                                                                                                                 <  1:681342/250‑1 (MQ=255)
                                                                                                                      ggagattaatgggggggAAAATAAACTTGGCGTAAAGGGCAACCTGAATCCGGCGGCCGTCGTCTGGGAAATGGGGGGGGGCAGTGGGGTGATAATGGCTAAAAGGACACCGCAGTGATGGGGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACa                                                                                                                       <  1:571658/238‑1 (MQ=255)
                                                                                                                         gAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCa                                                                                                                    >  1:498207/1‑250 (MQ=255)
                                                                                                                            tAATTGGGGGAAAATAAAACTGGGCGTAAAGGGCAACCTGAACCCAGCGGCCAGCGTCGGGGGAATGGGGGCGTGCAGTTGGGTGATAATGGCTAAAATGACCCCGCAGTGAGGGGGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAttt                                                                                                                 <  1:367704/250‑1 (MQ=255)
                                                                                                                              gtttGGGGGAAAGTTAACTTTGGTGTAATGGCATCCTCGAATCCAGCGGCCGGGTTGTGGGGAAGGGGGGTGGGGAGCTGGGTGATAATGCCAACAATGCCACCGGAGGGATGGTGGGCGTAAAATATAAGTTTTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTg                                                                                                                <  2:162950/249‑1 (MQ=255)
                                                                                                                              aTTTGGGGGAAAATTAAATTGGGCAAAAAGGGCACCCTGAATCCGGCGGCCGGGTTCGGGGGAATGGGGGCGGGCAGGGGGTTAAAAAGGGCAAAAATGACCCCGCAGTGAGGGTGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAttt                                                                                                                 <  1:251204/249‑1 (MQ=255)
                                                                                                                                  ggggaAAATAAACCGTGGCGAAAAGGGAACCCTAAACCCGGCGGCCGGCGTCGGGGGAATGGGGGGGTGCGGTGGGGTAATAATGGCTAAAATGACCCCGCAGTGATGGGGGGCCTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGt                                                                                                           <  1:116751/250‑1 (MQ=255)
                                                                                                                                               aaCTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCa                                                                                                 >  1:25947/1‑249 (MQ=255)
                                                                                                                                                    gtcGTAAATGCAACCCGAAGTCCGGCGGCCGGCGTCGGGGGATGGTGGGGGTGCAGGTGGGTGAAAATGGTAACAATGCCCCCGAGGTGGTGGGGGGCCTAAAATAAAATTTCGGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTg                                                                                           <  2:32720/248‑1 (MQ=255)
                                                                                                                                                    ggCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGAGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCt                                                                                            >  1:666704/1‑249 (MQ=255)
                                                                                                                                                    ggCGAAAAGGGCACCCTGAACCCAGCGGCCAGCGTCTGGGAAATGGGGGCGGGCAGCTGGGTAATAATGGCTAAAATGACACCGCAGTGATGGTGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAAc                                                                                              <  1:554300/247‑1 (MQ=255)
                                                                                                                                                        tAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGAtt                                                                                        >  1:488192/1‑249 (MQ=255)
                                                                                                                                                         aaaTGAAAACCTGAACCCAGGGGCCAGCGTCGGGGGAATGTGGGCGTGCAGGTGGGTGATAATGGCCACAATGACCCCGCAGTGATGGTGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCa                                                                                      <  1:655762/250‑1 (MQ=255)
                                                                                                                                                                      gAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCa                                                                          >  1:691631/1‑250 (MQ=255)
                                                                                                                                                                            cAGCGGCCAGCGTCGGGGGAAGGTGGGCGTGCAGTTGGGTGATAATGGCTACAATGACACCGCGGTGAGGGTGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTg                                                                     <  1:62840/250‑1 (MQ=255)
                                                                                                                                                                             aGCGGCCAGGGTCTGGGGAGTGTGGGCGTGCAGCGGGGTGAAAATGGCTACAATGACACCGCAGTGAGGGGGGGCCTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTgg                                                                    <  1:668239/250‑1 (MQ=255)
                                                                                                                                                                                 gcgaGCGTGTGGGGAGTGTGGGCGTGGAGCTGAGTAATAATGCCTACAAGGACACCGCAGTGATGGTGGGCCTAAAATATATGTTTGGCCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCt                                                                <  2:239875/247‑1 (MQ=255)
                                                                                                                                                                                              gaaggggggCGTGCAGGTGGGTGATAATGGCTAAAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGAt                                                   <  1:253768/244‑1 (MQ=255)
                                                                                                                                                                                                          gCAGCTGGGTAATAATGACTAAAATGACCCCGCAGTGTTGGTGGGCCAAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAAttttt                                       <  1:302492/250‑1 (MQ=255)
                                                                                                                                                                                                                 ggTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACCGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCAc                                >  2:540391/1‑250 (MQ=255)
                                                                                                                                                                                                                        aTGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTgg                          >  1:415039/1‑249 (MQ=255)
                                                                                                                                                                                                                             aaCAATGACCCCGCAGTGATGGTGGCCCTAAAATAAAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGc                    <  1:406014/249‑1 (MQ=255)
                                                                                                                                                                                                                                         gCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGcattaccca        >  1:537693/1‑250 (MQ=255)
                                                                                                                                                                                                                                         gCAGTATTGGGGGGCATAAATTATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGcattaccca        <  1:223949/250‑1 (MQ=255)
                                                                                                                                                                                                                                            gTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGcattacccatta     >  1:219408/1‑250 (MQ=255)
                                                                                                                                                                                                                                             tGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGcattacccattac    >  1:29934/1‑250 (MQ=255)
                                                                                                                                                                                                                                                tggtggGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGCATTACCCATTACCg  >  1:169525/1‑249 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                              
AGGAGAAAAACAGCGGTCCGCACGGTGACGACTTTGAACCTTTTGTTGAAATGAACTGGATCCATAACAGTAAAGATTTTGCTGTCTCAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGCATTACCCATTACCG  >  NZ_CP009273/1217443‑1217927

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