Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,340,9020GA34.9% 75.1 / 10.6 43S33F (TCT→TTT) pdeRcyclic di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (25/3);  new base A (0/15);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 5.38e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAG  >  NZ_CP009273/1340657‑1341149
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       
gCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGaa                                                                                                                                                                                                                                                       >  1:665822/1‑248 (MQ=255)
               gCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTAAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGGTCGGCAGATAACTGCAATGTTCGCACAGGGTCTGTGGTTTCGGTGGGAGAAAAGCGCAAGACactgct                                                                                                                                                                                                                                      >  2:156497/1‑245 (MQ=255)
                cccGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCgctg                                                                                                                                                                                                                                     >  1:17057/1‑250 (MQ=255)
                    tgtgATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCttt                                                                                                                                                                                                                                 >  2:746219/1‑250 (MQ=255)
                              cccGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCg                                                                                                                                                                                                                       <  1:81252/248‑1 (MQ=255)
                                         ccggtggtttCATCCCCTGCAACCGCGTGCCCACCCCTTCCCCTTTTTTTGTTTTTCCTTCCAACCAAAGGCACAGAAAGTTCGTTTCCAGCGCCGGTCAGACTCACCATCGGGCGGGGGGTGATGACCGTCATTCCCCGGATGCGAGCGGCCTGTGCGGCAGAAAACGGCAATGTACGATCAGGTTCTGTGTTTTCGGTGGTAAAAAAGCGAAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGAt                                                                                                                                                                                                            <  2:858773/241‑1 (MQ=255)
                                                                  cgTTGCCAGCCCATCCCCTTTTATTGTTTTTTCGCCCTACCAGATGCACAAAAAGATCGTTTTCATCGACGTCCAGACTCATCATCAGGCTGGAGGTGATGCCCTTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGTTTCTGTGTTTTGGGTGGTAAAAAAGCGAAAAACATCGCTGCTTCCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACaa                                                                                                                                                                                   <  1:856882/250‑1 (MQ=255)
                                                                     gcccACCCCTTTCCTGTTTATTTATTTTCCGCCCTACCATATACACAAAAAGTTCTCTTCAATCGACGGACAATCACACAATCGGGGGGGATGTGATGCCCTTCTTTTCCCTGATGCGAGCAGCCTGTTCGGCAAATAATGCCAATGTCCGATTGGGTTTTGTGGTTTGGGTGGTAAAAAAGCGCAAAACACCGCTGTTTTCCTTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAAc                                                                                                                                                                                  <  2:868240/246‑1 (MQ=255)
                                                                         accccATTCCCTTTATTGGATTTTCCGTCCTACCAGAGGCAGAAAAAGACCCTTTTCATCGCCGGTCAGCCTAACCATAAGGCGGGAGGTGATGACCGTCTTGTCCGTGATGCGAGCAGCCTTTTGGGAAGAAACCTGCAATGTACGACCAGTTTCGGTGGTTTGGGTGGAAAAAAAGCACAAACCATCGCTTTTTTCCGTAAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTtgtg                                                                                                                                                                            <  2:693841/248‑1 (MQ=255)
                                                                          ccccATTCCCTTTTATTGTTTTTCCGCCCTACAAGATGCAAAAAAAGATCGCTTTCATCGGCGGTCAGCCTCACCATCAGGTTGGGGGTGATGCCCGTCTTTCCCCGGATGCGAGCGGCTTGTTCGGAAAATAACTGCAATTTCCGATCGGGTTCTGTGGTTTCGGTGGTAAAAAAGCGAAAAACATCGCTGTTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTgg                                                                                                                                                                           <  2:144370/249‑1 (MQ=255)
                                                                          ccccATTCCCTTTTATTGATTTTTCGTCCTCCCAGATGAACAAAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGCCCGTCTTTCCCCTGATGCGAGCAGCCTGTTCGGAAGATAACTGCAATGTACGATCAGTTTCTGTGTTTTGGGTGGTAAAAAAGCGCAAACCATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTgg                                                                                                                                                                           <  1:312683/249‑1 (MQ=255)
                                                                            ccATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGAc                                                                                                                                                                         >  1:470169/1‑250 (MQ=255)
                                                                               ttCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCACTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAAAATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTcc                                                                                                                                                                      >  2:201311/1‑250 (MQ=255)
                                                                                tCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTccc                                                                                                                                                                     >  1:101602/1‑250 (MQ=255)
                                                                                                       cTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGCTGCatat                                                                                                                                               >  2:783416/1‑250 (MQ=255)
                                                                                                         cccAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCGGGAGGTGATGACCGTCTTTCCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGTTTCTGTGTTTTGGGTGGTAAAAAAGCGCAAACCATCGCTGCTTCCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATAtcc                                                                                                                                             <  1:80347/249‑1 (MQ=255)
                                                                                                                         gATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCg                                                                                                                               >  2:696968/1‑248 (MQ=255)
                                                                                                                               tttCATCGACGGTCAGACTCATCATCAGGCTGGCGGTGATGACCGTCATTTCCCTGATGAGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAAACGCCAGTAAGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCcgt                                                                                                                       >  1:660594/1‑249 (MQ=255)
                                                                                                                                     cGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTta                                                                                                                  >  1:273513/1‑249 (MQ=255)
                                                                                                                                     cGACGGCAAGACTCATCATGAGGCGGGGGGTGATGACCGTCATTCCCCGGAGGCGAGCACCCTGTTCGGAAGATAACTGAAATGTACGATCGGGTTCTGTGGTTGGGGGGGAAAAAAAGCGAAAAACATCGCTGTTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaa                                                                                                                 <  1:485955/250‑1 (MQ=255)
                                                                                                                                           tCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCTGTtaaataa                                                                                                             >  1:824390/1‑248 (MQ=255)
                                                                                                                                           tCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaataa                                                                                                             >  1:529630/1‑248 (MQ=255)
                                                                                                                                            cAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGTGGTTAAATAAAAg                                                                                                          >  1:806020/1‑250 (MQ=255)
                                                                                                                                            cAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAg                                                                                                          >  2:691207/1‑250 (MQ=255)
                                                                                                                                                       tcaGGCTGGAGGTGATGACCGTCTTTCCCCTGATGCGAGCAGCTTGTCGGGCAGATAACTGCAATGTACGATCAGGTTCTGTGTTTTGGGTGGAAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGa                                                                                               <  1:559771/250‑1 (MQ=255)
                                                                                                                                                             tGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACa                                                                                         >  1:340025/1‑250 (MQ=255)
                                                                                                                                                                 ggTGATGACCGTCATTCCCCTGTTGCGAGCCGCCTTTTCGGAAGAAAACTGCATTGAAGGATCAGTTTTTGGGTTTTCGGTGGTAAAAAAGCGCAAAACACCTCTCTTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGaa                                                                                     <  1:878131/250‑1 (MQ=255)
                                                                                                                                                                                  cccgtaggCGAGCAGCCGGTGCGGAAGATAACTACAATGTACGATCAGGTTTGTGTGTTGCGGTGAAAAAAAAGCGAAAACCATCGCTTCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAAc                                                                    <  2:529630/243‑1 (MQ=255)
                                                                                                                                                                                     tGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCa                                                                 >  1:206913/1‑250 (MQ=255)
                                                                                                                                                                                               gtcttGTTCGGAAAATACCTGAAATGTACGATGAGGTTTTGTGGTTTGGGTGGAAAAAAAGCGCAAAACTCCGTTTTTTCCCATCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGttttttt                                                       <  2:81723/246‑1 (MQ=255)
                                                                                                                                                                                               agcCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGttttttt                                                       >  1:115752/1‑250 (MQ=255)
                                                                                                                                                                                               agcCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAAAGCATTGttttttt                                                       >  2:2414/1‑250 (MQ=255)
                                                                                                                                                                                                 ccTGTTCGGCAGATACCTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAACGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGACCGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAACAGGAACAACATGACTTACATGAAATTAACGGCGGCTCAACGCATTGgtttttttt                                                     >  1:141264/1‑249 (MQ=255)
                                                                                                                                                                                                                    cAATTTACGATCAGGTTCTGTGGTTTGGGTGGTAAAAAAGCGCAAAACTCCCCTGTTTCCCGTAAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCa                                  <  1:248207/250‑1 (MQ=255)
                                                                                                                                                                                                                       tgatcatttcggTTTTGGGGGTTCGTGTGGAAAAAAACCGAAAAACCCCCCTGCTTCCCACAACCCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGt                                <  2:428827/240‑1 (MQ=255)
                                                                                                                                                                                                                       tGTACGTTCAGTTTTTGTGGTTTGGGTGGTAGAAAAGCGAAAAACATGTCTGTTTTCCGTAAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGtt                               <  1:563117/250‑1 (MQ=255)
                                                                                                                                                                                                                            gATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTTTTTATTTATTTGTAATCCAGTTTTgg                           >  2:378626/1‑249 (MQ=255)
                                                                                                                                                                                                                                  gTTTTGTGTGTTCGGGGATAAAAAAGCAAAACACGGCTCTTCTTTCCCACAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGTGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAAc                     <  2:801105/249‑1 (MQ=255)
                                                                                                                                                                                                                                       gtgGTTGGGGTGGTAAAAAAGCAAAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAg                <  1:56390/249‑1 (MQ=255)
                                                                                                                                                                                                                                               cGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCAAGTTTTTaa         >  1:194971/1‑249 (MQ=255)
                                                                                                                                                                                                                                                   tggtAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACg      >  1:647696/1‑249 (MQ=255)
                                                                                                                                                                                                                                                   tggtAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAAAGCCAGTTTTCAAACg      >  1:810597/1‑249 (MQ=255)
                                                                                                                                                                                                                                                      tAGAAAAGCGCAAAATATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCCTTTTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAg  >  2:497583/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       
GCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAG  >  NZ_CP009273/1340657‑1341149

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