Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,001,8540GT28.3% 140.3 / 13.2 61G33C (GGC→TGC) yedEselenium metabolism membrane protein YedE/FdhT
Reads supporting (aligned to +/- strand):  ref base G (16/27);  new base T (17/0);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 4.15e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATTAAAGCCCTGATGTCGCCAGCCTTTCAAGGTTCAGAGAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGG  >  NZ_CP009273/2001608‑2002100
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         
aTTAAAGCCCTGATGTCGCCAGCCTTTCAAGGTTCAGAGAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTGTTCTCTCTACCTACTATTTTTGcc                                                                                                                                                                                                                                                        >  2:149402/1‑249 (MQ=255)
    aaGCCCTGATGTCGCCAGCCTTTCAAGGTTCAGAGAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTa                                                                                                                                                                                                                                                     <  2:133686/249‑1 (MQ=255)
    aaGCCCTGATGTCGCCAGCCTTTCAAGGTTCAGAGAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTAc                                                                                                                                                                                                                                                    <  1:559887/250‑1 (MQ=255)
               tCGCCAGCCTTTCAAGGTTCAGAGAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCtttt                                                                                                                                                                                                                                         <  2:520689/250‑1 (MQ=255)
                      ccTTTCAAGGTTCAGAGAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCTCCCATCCCTGCGGTCATCGCGGCGGGTAGTTTCTTTTCTTACCATCTTTGCCTTCATGGCCGCTTTTGTGgtgt                                                                                                                                                                                                                                 >  2:338682/1‑250 (MQ=255)
                      ccTTTCAAGGTTCAGAGAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGTTTCTCTCTACTTAATATTTTTTGATTTCTTGCCACTTTTGGGCCGt                                                                                                                                                                                                                                 >  1:343479/1‑250 (MQ=255)
                         ttCAAGGTTCAGAGAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTTTCTCTCCTTACTTTTTTTGCATTACTTGCACCTTTTTGGCTGTccc                                                                                                                                                                                                                              >  1:477280/1‑248 (MQ=255)
                                   agagagAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTTTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAAttt                                                                                                                                                                                                                    >  1:522270/1‑250 (MQ=255)
                                     agagaAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGGGAAttttc                                                                                                                                                                                                                  >  1:465351/1‑248 (MQ=255)
                                                aTCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGAGTTTCTCTACTTACTATTTTTGCCTTTCTGGCCCCCTTTGTGTTGTCTCCGGTGAATTTTCCCGCTgtgggg                                                                                                                                                                                                       >  2:407564/1‑248 (MQ=255)
                                                   tAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTACTTTCTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTAACCGTTGGGGCGGcc                                                                                                                                                                                                    >  1:210561/1‑250 (MQ=255)
                                                     aCTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGggcgg                                                                                                                                                                                                      <  2:400890/246‑1 (MQ=255)
                                                           attTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTTGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTg                                                                                                                                                                                           >  1:255251/1‑250 (MQ=255)
                                                             ttttGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTCTCTTCTTTCCATTTTTTGTTTTATTGGACCTTTTTGGGTGTTCTGGGGTGATTTTCCCCTTGGGGCGGCCCGCCCCTGc                                                                                                                                                                                          >  1:156604/1‑249 (MQ=255)
                                                              tttGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCa                                                                                                                                                                                         <  1:398897/249‑1 (MQ=255)
                                                                     attattACCCCGAAAATAAATTTTCATGGTTTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCg                                                                                                                                                                                 <  2:274941/250‑1 (MQ=255)
                                                                      ttattaCCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTTTAATTTCTATTTTTGCATTACTGGCCCCTTTTTGGCTGTCCCGGGTGAATTTTTCCGTTTGGGCGGGCCGCTCCCGCAACTGtttcg                                                                                                                                                                                >  2:384300/1‑247 (MQ=255)
                                                                        atttCCCCGAAAATAAATTTGCATGGTGTAGGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCGTGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGc                                                                                                                                                                               <  2:545109/249‑1 (MQ=255)
                                                                         ttaCCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCg                                                                                                                                                                              <  1:171468/249‑1 (MQ=255)
                                                                          taCCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTACTTACTATTTTGGCATTTCTGGCACCTTTTTGGCTGTCACGGGGGGATTTTCCCGTTTGGGGGGCCCGCTCCTCCACCTGTTCGGCggc                                                                                                                                                                            >  2:229561/1‑248 (MQ=255)
                                                                               cGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTTGCATTACTGGCACCCTTTTGGCTGTCACGGGGGAATTTACCCGTTGGGGCGGCCAGCCCCTGCAACTGTTCGGCGTCCATGc                                                                                                                                                                       >  1:350137/1‑250 (MQ=255)
                                                                                gAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCt                                                                                                                                                                      <  1:224674/250‑1 (MQ=255)
                                                                                gAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATGCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCt                                                                                                                                                                      <  1:273033/250‑1 (MQ=255)
                                                                                 aaaaTAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTg                                                                                                                                                                     <  2:616251/250‑1 (MQ=255)
                                                                                 aaaaTAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGGGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTTGGCGTCCATGCt                                                                                                                                                                      >  1:38470/1‑249 (MQ=255)
                                                                                  aaataaatTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGa                                                                                                                                                                    <  2:83472/250‑1 (MQ=255)
                                                                                   aataaatTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACTTTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTTTCTACTTTCTATTTTGGCATTACTGGCACCTTTTTGGCTGTCACGGGTGAATTTTCCCGTTGGGGCGGCCAGCTCCTGCAACTGTTTGGCGTCCATTCTg                                                                                                                                                                     >  1:55051/1‑248 (MQ=255)
                                                                                     taaatTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGAATCTCTCTTCTTTCTATTTTTGCCTTTCTTGCCCCCTTTTGGCTTTTCCGGGTGAATTTTTCCGCTTGGGCGGCGCGCTCCTGCAACTGTTCGGCGTCCATGCttac                                                                                                                                                                   >  2:574417/1‑245 (MQ=255)
                                                                                       aatTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAagag                                                                                                                                                                <  1:221400/249‑1 (MQ=255)
                                                                                             cATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGtt                                                                                                                                                         <  2:133911/249‑1 (MQ=255)
                                                                                             atgggTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGGCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGtt                                                                                                                                                         <  2:288984/247‑1 (MQ=255)
                                                                                                 gtgtATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTCTCTCCTTTCTATTTTTGCATTTCTGGCCCCCTTTGGGCTGTCACGGGTGGAATTACCCGTTTGGGCGGGCCGCCCCTGCCACTGGTTCGCGTCCCTGCTGCAGCGGGGGGTCACtt                                                                                                                                                     >  2:261645/1‑250 (MQ=255)
                                                                                                  tgtATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACttt                                                                                                                                                    <  1:600725/250‑1 (MQ=255)
                                                                                                   gtATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTCTCTTCCTACTATTTTTGCATTACTGGCACCTTTTTGGCTTTCACGGGGGAATTTTACCGTTGGGGCGGCCAGCTCCTGCCACTGTTTGGCGTCCATGCTGAAGAGTGGGGTTTCttt                                                                                                                                                    >  1:553989/1‑249 (MQ=255)
                                                                                                        aTTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAAtt                                                                                                                                              <  1:634836/250‑1 (MQ=255)
                                                                                                           cgcgACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTCCTTACTATTTTGTCCTTACTGTCACCCTTTGGGCTGTCTCGGGTGAAATTACCCCGTGGGGCGGCCCTCCCCTCCATCTGTTCTGCTTCCCTTCTGCAGAGTGTGGTGACTTTtaaagtttc                                                                                                                                           >  2:152242/1‑245 (MQ=255)
                                                                                                                ctctTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTTTCTCTACTTACTATTTTTGCATTACTGGCACCTTTTGGGCTTTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTTAAATTATCCATCt                                                                                                                                      >  1:179797/1‑250 (MQ=255)
                                                                                                                 tctTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTCTCTACTTACTATTTTGGCCTTACTGGCACCCTTTTGGCTTTCACGGGGTAATTTTCCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGGCCCTTCTTGAGAGGGGGGGTACTTTTAAATTATCCCTCTg                                                                                                                                     >  2:187315/1‑250 (MQ=255)
                                                                                                                   ttGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGTTTTTCTCTACTTACTCTTTTTGCATTACTGGCACCCTTTTGGCTGTCTCGGGTGGATTTACCCGTTTGGGCGGCCAGCTCCTGCCACTGTTTGGGGTCCACGCTGGAGAGTGGGGTTACTTTTAAATTTTCCATCTgg                                                                                                                                    >  2:500598/1‑249 (MQ=255)
                                                                                                                      tttCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGTATATTTCTTCTTTCCCTTTTTGCGTTACTTGCACCCCTTGGGGTTGCCCGGGGGGATTTACCCCTTGGGGCGGCCAGCCCCTCCCACCTGTTCGCGGCCCTGCTGATGGGGGGGGGTTCCTTACAAATATTCCTCTGGGAgg                                                                                                                                >  2:384072/1‑250 (MQ=255)
                                                                                                                         cAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGGTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGAt                                                                                                                              <  2:501985/249‑1 (MQ=255)
                                                                                                                          aaaCCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTTTCTCTACTTCCTCTTTTTGCATTACTGGCACCTTTTTGGCTGTCACGGGGGAATTTATCCGTTGGGGCGGCCAGCTCCTGCAACCGTTCGGCGTCCCTGCTGCAGAGTGGGGGTACTTTTAAATTCTCCATCTTGAAGGATCg                                                                                                                            >  1:469676/1‑250 (MQ=255)
                                                                                                                             ccTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCa                                                                                                                         <  2:30042/250‑1 (MQ=255)
                                                                                                                                         agagGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTTCTTACTATTTTTGCATTACTGGCACCCTTTTGGCTGTCACGGGTGGATTTTCCCGTTTGGGCGGCCAGCTCCCTCCAATGTTCGGCGGCCATGCTGAAGAGTGGGGGTTCTTTTAAATTTTCCCTCTTGAAGGGTCGCCATTTAAccccctc                                                                                                             >  2:60240/1‑246 (MQ=255)
                                                                                                                                                         tcatgGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATc                                                                                             <  1:626235/249‑1 (MQ=255)
                                                                                                                                                                       cAAACCCGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACATCCTATTTTCGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGCATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTgg                                                                               <  1:178947/249‑1 (MQ=255)
                                                                                                                                                                         aaCACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTggcgg                                                                            <  1:147978/250‑1 (MQ=255)
                                                                                                                                                                            acGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTCTCTACTTACCATTTTTGCATTACTGGCACCCTTTTGGGTGTCACGGGGGAATTTTCCCGTTGGGGGGGCCAGCTCCTGCAACTGTTCGGCGTCCATGGTGAAGAGGGGGGGTACTTTTAAATTATCCATCTGGAAGGATCGCCATTAACCCGCCTCGCAGGGAGGTTGATCCCCGGGATGTTTGGGGgctg                                                                         >  2:108048/1‑250 (MQ=255)
                                                                                                                                                                                  ggTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGGCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGGTGATGATCCTCGGTATGTTTGGCGGCTGCTTTGc                                                                   >  1:280706/1‑250 (MQ=255)
                                                                                                                                                                                                        cATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTAATTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGGGAATTTACCCGTTGGGGGGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGGTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAAACCGCATCGACGGGGTGATGATCCTCGGGATGTTTGGCGGCTGCTTTGCCGCCGCGCTGGGGGCCAACAAt                                             >  1:559426/1‑250 (MQ=255)
                                                                                                                                                                                                                            cggGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGTTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCcgcg                         >  1:378056/1‑250 (MQ=255)
                                                                                                                                                                                                                             gggTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCgcgc                        <  2:412062/250‑1 (MQ=255)
                                                                                                                                                                                                                                  ttCTCTCTACTTACTATTTTGGCATTACTGGCTCCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGcc                    <  1:151814/249‑1 (MQ=255)
                                                                                                                                                                                                                                    ctctctACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTa                 >  2:556270/1‑250 (MQ=255)
                                                                                                                                                                                                                                          acttactATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTCTCCGCa           >  2:211382/1‑250 (MQ=255)
                                                                                                                                                                                                                                            ttactATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCAt          >  1:167986/1‑249 (MQ=255)
                                                                                                                                                                                                                                               ctATTTTGGCATTACTGGCACCTGTTGGGCTGTCACGGGTGAATTTTCCCGTTGGGGCGGCCAGCTCCGGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGTTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATg      <  1:192826/250‑1 (MQ=255)
                                                                                                                                                                                                                                               ctATTTGGGCATTACTGGCTCCTTTTGGGCTGTCCCGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGTGGTTACTTTAAAATTATCCATCGGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTAt       <  2:562481/249‑1 (MQ=255)
                                                                                                                                                                                                                                                   ttttGGCATTACTGGCGCCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCAAATGCCGCGCAGCCGTATCCGCATTATGc     <  1:229561/248‑1 (MQ=255)
                                                                                                                                                                                                                                                   ttttGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCa    <  2:553989/249‑1 (MQ=255)
                                                                                                                                                                                                                                                    tttGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAgg  >  2:37153/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         
ATTAAAGCCCTGATGTCGCCAGCCTTTCAAGGTTCAGAGAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGG  >  NZ_CP009273/2001608‑2002100

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