Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,560,4950CT34.7% 118.9 / 23.1 72T301T (ACG→ACAdosCdiguanylate cyclase DosC
Reads supporting (aligned to +/- strand):  ref base C (7/40);  new base T (25/0);  total (32/40)
Fisher's exact test for biased strand distribution p-value = 2.21e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGT  >  NZ_CP009273/1560249‑1560733
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
aGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGt                                                                                                                                                                                                                                                 <  1:425619/249‑1 (MQ=255)
      aaCAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGCTAAGTAATTTTGTCatgaca                                                                                                                                                                                                                                          >  2:295210/1‑245 (MQ=255)
         aaTGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTCAAGTAATTTTGttcgcacctca                                                                                                                                                                                                                                       >  2:296923/1‑240 (MQ=255)
             aTAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATAc                                                                                                                                                                                                                                   <  1:121950/250‑1 (MQ=255)
               aaaTTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGGCTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCg                                                                                                                                                                                                                                 <  2:349268/250‑1 (MQ=255)
                  ttCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACt                                                                                                                                                                                                                              <  2:63477/250‑1 (MQ=255)
                          ccccGTAGCGGAAAACATAATCGCTACTGCGGACGATGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCACACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGtcg                                                                                                                                                                                                                      <  2:271159/250‑1 (MQ=255)
                                     aaaaCATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGGTAAATAATTTCTTCAGTAAATCCCTACCGACTTCCTGGCGCGAGActtttt                                                                                                                                                                                                           >  2:68316/1‑250 (MQ=255)
                                       aaCATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCa                                                                                                                                                                                                         <  1:133380/250‑1 (MQ=255)
                                                  tactGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAAATTTGTTAGTACAACCATAACGACCTCGTGGGGCGGTACTTTTTTAAACAAATTACg                                                                                                                                                                                              >  2:234962/1‑250 (MQ=255)
                                                       cGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACgcag                                                                                                                                                                                          <  1:478531/248‑1 (MQ=255)
                                                          aCGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGGTTAATAAATTTGTCAGTACATCCATACCGCCTTCGGGGGGCGGTAATTCTTTAAAAAAATAACACAACAAgg                                                                                                                                                                                    >  2:288894/1‑250 (MQ=255)
                                                                tCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGtaata                                                                                                                                                                               <  1:159933/249‑1 (MQ=255)
                                                                           cTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTTAGTAATTTTGTCCGTACATCCATACCCACTTCGTGGCGCGGTACTTCTTCAAACACATTACTCCGCAAGGCAATAATTTTTGATAct                                                                                                                                                                   >  2:697179/1‑249 (MQ=255)
                                                                           cTGAGAGAATTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGACACGGCGGCTTAGTAATTTCGTCGATAAATCCATATCGACATCGTGGGGGCGTACTTTTTAAAACAACTCACACACCAAGAGaataatatttgaaacg                                                                                                                                                                   >  2:683991/1‑240 (MQ=255)
                                                                            tGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTTCGTACATACATACCGACTTCGTGGGGCGATTATTTTTCTAAAAATTCACGCAGCAAGGTAATAATTTTCGATAcag                                                                                                                                                                  >  2:148805/1‑248 (MQ=255)
                                                                              agagaCTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTTAGTAATTTTGTTAGGACATCCATAACGACTTCGGGGCGCGATAATTTTTCAAACAATTCACGCAGCAAGGGAATAAttttggaaaaggta                                                                                                                                                                >  1:711151/1‑239 (MQ=255)
                                                                                agaCTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTAtt                                                                                                                                                             <  1:446013/250‑1 (MQ=255)
                                                                                    tttACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTTAGTAATTTTTTCAGTAACTCCCTACCGACTTCGGGGCGCGGAACTTCTTCTAAACATTACCGACGCAATGTAATAATTTTTGATACCGTAtttttt                                                                                                                                                         >  2:24325/1‑250 (MQ=255)
                                                                                           aGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGta                                                                                                                                                  <  1:492362/250‑1 (MQ=255)
                                                                                            gAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTCCGTACATCCATAACGACTTCGTGGGGCGATACTTTTTCCAACAATTCACGCAGCAAGGGAAAAATTTGCGATACGGTATTTTTTATTTGtt                                                                                                                                                  >  1:740760/1‑248 (MQ=255)
                                                                                            gAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCAGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGtaa                                                                                                                                                 <  2:30320/250‑1 (MQ=255)
                                                                                             aaTTTCATCACCAGAGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGtaat                                                                                                                                                <  1:424424/250‑1 (MQ=255)
                                                                                                   atcaCCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATaaaaa                                                                                                                                           <  2:102080/249‑1 (MQ=255)
                                                                                                      aCCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTc                                                                                                                                       <  1:389080/250‑1 (MQ=255)
                                                                                                       ccAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTTCGGACATCCAAAACGACCTTGTTGGGCGATACTTCTTTAAACAATTTACCCAGCAAGGGTAAAATTTTCGGTACGGGTTTTCTTTTCTGTAATAAAAAattta                                                                                                                                      >  1:147428/1‑248 (MQ=255)
                                                                                                          gtgtTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAAATTTGTCAGTAAAACCATACCGACTTTGTGGCGCGATACTTCTTCAAAAAAATTACGCAGCAAGGGAATAATTTTCGATACGGGATTTTTTATCTGTAATAAAAAATTTAccc                                                                                                                                   >  1:342255/1‑250 (MQ=255)
                                                                                                             ttATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCa                                                                                                                                <  1:595086/250‑1 (MQ=255)
                                                                                                                         gTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTTAGTAATTTTGTCAGTACATCCATACCAACTTCGTGGCGCGATAATTCTTCCAACAATTCACGCAGCAAGGGAATAATTTTCGATACGGGATTTTTTTTCTTTAATAAAAATTTTACCCCCCGACTTTTGtttt                                                                                                                    >  1:118135/1‑248 (MQ=255)
                                                                                                                            tCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACAGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTc                                                                                                                  <  2:18352/249‑1 (MQ=255)
                                                                                                                                    ctctTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTCTGGACATCCATACCGACTTCGTGGCGCGATAATTCTTTAAACAATTTACCCAAGCAGGGTATAATTTGTGATCACGTATTTTTTATTTTTTATAAAAATTTTATCCACAGACCTCTCTTTTTTCAATTTCTg                                                                                                         >  2:580178/1‑250 (MQ=255)
                                                                                                                                       tttGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTg                                                                                                         <  1:286992/247‑1 (MQ=255)
                                                                                                                                       tttGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTg                                                                                                         <  2:212381/247‑1 (MQ=255)
                                                                                                                                       tttGAAGTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTc                                                                                                           <  1:526078/245‑1 (MQ=255)
                                                                                                                                         tGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAAATTTGTTAGTACATCCATACCCACTTCGTGGCGCGGGACTTTTTTCAACAAATTACGCAGCAAGGGAATAATTTGCGATACGGTATTTCTTATCTGTTATAAAAATTTTACCCCCCGACTTCTGTTTTTTAAATTTTTGGTGtt                                                                                                    >  2:11036/1‑250 (MQ=255)
                                                                                                                                         tGAATTTATCAACGTCAAGAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATTCTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGgtg                                                                                                      <  1:743336/248‑1 (MQ=255)
                                                                                                                                                aTCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGATAAAAAATTTCGTAAGTATCTCCATCCCGACTTCCTGGCGCGGTAATTTTTCAAACAATTCACTCAGCAAGGTAATAATTTTCGATACGGTATTTCTTATCTTTTATAAAAAATTCTACCGCAGCATTCTGTTATTTTAATTTTTTGTGTTATGGATg                                                                                             >  2:414039/1‑250 (MQ=255)
                                                                                                                                                   aaCGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTTAGTAAATTTGTCCGGACATACACCCCGACTTCCTGGCGCGGTAATTCTTCAAACAATTAACGCAGCAAGGGAATAATTTTCGATATGGGTTTTTTTTTTTTTAATAAAAAATTTTCACACAGAAATATTTTTTTTCAATTTTTTGTGTTATCCATGGtt                                                                                          >  2:337864/1‑250 (MQ=255)
                                                                                                                                                                aGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTTAAGAATTTTTTCAGTACATCCATAACGACCTCGTGGCGGGGAAATTTTTTAAACAAATCACCCAGAAAAGGACAAAATTTGTATATGGTTATTCTTTTTTTTTAATAAAAATTTAACCCCACGACTTTTGTATTTAAATTTTTTGTGTTTCGTATGTGTTGATTTAAAAtta                                                                             >  2:409908/1‑248 (MQ=255)
                                                                                                                                                                    cTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTTAATAAATTTGTCAGTACATTCAAACACAAATCGGGGGGCGATAATTCTTTCAAAAAATCACCCAGCAAGGGAAAAATTTTGGAAACGGGATTTCTTTTTTTTAATAAAAATTTCACCCGCACACTTCTTTTATTCAAATTTCTGGTGTTATGCATGGGTTGTTTGAAAATACAGTc                                                                         >  2:437948/1‑250 (MQ=255)
                                                                                                                                                                         aGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAAt                                                                    <  1:295210/250‑1 (MQ=255)
                                                                                                                                                                          gCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATc                                                                   <  2:211829/250‑1 (MQ=255)
                                                                                                                                                                          gCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGATTTGATTGAAAATACCGTCGAAATc                                                                   <  2:355509/250‑1 (MQ=255)
                                                                                                                                                                                tACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAAt                                                             <  2:140751/250‑1 (MQ=255)
                                                                                                                                                                                    ggtccggtTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGa                                                         <  1:296923/250‑1 (MQ=255)
                                                                                                                                                                                         ggtTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGga                                                    <  1:589744/250‑1 (MQ=255)
                                                                                                                                                                                         ggtTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGga                                                    <  1:476382/250‑1 (MQ=255)
                                                                                                                                                                                           ttGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGGTAAGTAATTTTGTCAGTACAACCATACCGACTTCGTGGCGCGGGACTTCTTTAAAAAAATCACGCAGCAAGGGAATAAATTGCGATACGGGTTTTCTTTTATGTAATAAAAATTTCAACCCCCGACATCTGTTGTTTTAATTTCTTGTGTTACGCATGGTTTGATTGAAAATAACGTCTAGATCCTGCAACAGAAGGgagg                                                  >  2:73253/1‑249 (MQ=255)
                                                                                                                                                                                              gCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGatat                                               <  1:50914/250‑1 (MQ=255)
                                                                                                                                                                                                aTGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTCAGTACATCCATAACGACTTCGTGGCGCGATACTTTTTTAAACAAATTACGCAGCAAGGGAATAATTTGCGATAAGGGATTTCTTTTCTGTAATAAAAATTTCAACCGCAGACTTCTGTTATTCAAATTTTTGGTGTTACGCATGGGTTGGTTTAAAATACCGTCGAAATCCTGAATCAGAAGGGGGATATgg                                             >  1:42864/1‑250 (MQ=255)
                                                                                                                                                                                                    gCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGGAAATTCGTCAGTACATCCATACCGACTTCGTGGCGCGGTACTTCTTCAAACAAATCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTa                                         >  1:717806/1‑250 (MQ=255)
                                                                                                                                                                                                    gCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAATTATTTTGTTAGTAAAACCATACCGACTTCGTGGGGCGGTAATTCTTTAAACAATTCACGCAACAAGGGAATAATTTTCGATACGGGATTTTTTTTCTTTAATAAAAAATTTACCCCCCAACTTCTTTTTTTTTAATTTTTTGTGTTTTGCCTTGTTTGATTTGAAATAACGTTGTAATCCTTAATTAGAAGGGAGatattggctt                                         >  2:163405/1‑244 (MQ=255)
                                                                                                                                                                                                       aTTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACtt                                      <  2:495999/250‑1 (MQ=255)
                                                                                                                                                                                                            gcgTTTGAAGATAGTCGGTAGGAAACGGCGGTTTAGAAAATTTGTCTGTACATCCATACCGACCTCGTGGGGCGGTACTTCTTTAAACAAATCACCCAGCAAGGGAATAAATTTCGGTACGGTATTTTTTTTCTTTAATAAAAAATTTAACCCCAGACTTTTGTTTTTTAAATTTTTTGGGTTTCGCATGGTTTGCTTGAAAATATCGTCCAAATACTGAATCAGACGGGGGATATGGGCCACTTCggg                                  >  2:626732/1‑248 (MQ=255)
                                                                                                                                                                                                            gcgTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAATAATTTTTCCAGTAAAACCATACCCCATTCCTTGGGCGATACTTCTTTTAACAATTTACGCAGCAAGGTAATAAATTGCGATACGGAATTTCTTTTCTGTAATAAAAATATCAACCTCACACCTCTGTTATTCAAATTTTTGGTGGTTCACATGGTTTGAATGAAAATAACGTCTAAATACTGAAACAAACGGGCGCTTTTGCCTACTTCggg                                  >  1:635926/1‑248 (MQ=255)
                                                                                                                                                                                                              gTTTGAAGATAGTCGGTAGGAAACGGCGGTTTAGGAATTTTTTTAGTAAATCCATAACGACTTAGTGGCGCGGGACTTCTTTAAACAATTTACGCAGCAAGGTAATAATTTGCGTTAAGGGTTTTTTTATTTGTTATAAAAATTTCACCCTCAGACTTCTGTTTTTTTAATTTTTTGTGTTAAGCCTGGTTTTATTGAAAATAACATCGAAATTCTGAATTAGACGGGAGAAATGGCTTACTTCTGCaa                                >  2:310065/1‑249 (MQ=255)
                                                                                                                                                                                                              gTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAAt                               <  1:201902/250‑1 (MQ=255)
                                                                                                                                                                                                              gTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAAt                               <  1:237263/250‑1 (MQ=255)
                                                                                                                                                                                                               tttGAAGATAGTCGGTAGGAAACGGCGGTTAAATAAATTTTTCAGTACAACCATAACGACTTCGGGGGGCGGTAATTCTTTAAACAATTTACCCAGCAAAGGAATAATTTTCGATAAGGTATTTTTTTTTTGTAATAAAAATTTTACCCGCAGACTTCTGTTATTCTAATTTTTTGTGGTATGCCTGGTTTTATTTAAAATACCGTCGAAATTCTTGATCAAACGGGGGATATTGGCTAATTCCGCCATa                              >  2:230908/1‑250 (MQ=255)
                                                                                                                                                                                                                       tAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTaaa                      <  1:37745/250‑1 (MQ=255)
                                                                                                                                                                                                                          tCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAaata                   <  1:441768/250‑1 (MQ=255)
                                                                                                                                                                                                                            ggTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAaataa                  <  1:335787/249‑1 (MQ=255)
                                                                                                                                                                                                                                ggAAACGGCGGTTAAGTAAATTTGTCAGTACATCCATACCGACTTCGGGGCGCGAAACTTTTTCAAACAATTTACGCAGCAAGGTAATAATTTTCGATACGGTATTTTTTTTCTTTAATAAAAATTTCACCCGCAGACTTCTGTTTTTTAAATTTCTTGTGTTAAGCATGGTTTGTTTTAAAATACCGTCGAAATTCTGAATCAAACGGGAGATATGGCCTACTTTGGCCATACCACTAAAAAAATGTc              >  2:57897/1‑249 (MQ=255)
                                                                                                                                                                                                                                  aaaCGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGAc           >  2:544004/1‑250 (MQ=255)
                                                                                                                                                                                                                                  aaaCGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAAGAATGTCGAc           <  2:571597/250‑1 (MQ=255)
                                                                                                                                                                                                                                     cggcggTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCtt        <  2:94535/250‑1 (MQ=255)
                                                                                                                                                                                                                                      ggcggTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCttt       <  1:30907/250‑1 (MQ=255)
                                                                                                                                                                                                                                      ggcggTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCttt       <  2:128502/250‑1 (MQ=255)
                                                                                                                                                                                                                                         ggTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTAt                                                                                                                    >  1:588738/1‑138 (MQ=255)
                                                                                                                                                                                                                                         ggTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATg    >  2:293226/1‑250 (MQ=255)
                                                                                                                                                                                                                                          gTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTAtt                                                                                                                   <  2:588738/138‑1 (MQ=255)
                                                                                                                                                                                                                                           ttAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGt  <  2:541997/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
AGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGT  >  NZ_CP009273/1560249‑1560733

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: