Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,722,0120AC74.3% 63.5 / 20.8 35G240G (GGT→GGGxylBxylulokinase
Reads supporting (aligned to +/- strand):  ref base A (6/3);  new base C (16/10);  total (22/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGA  >  NZ_CP009273/3721879‑3722149
                                                                                                                                     |                                                                                                                                         
aaCGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACCCCAAc                                                                                                                                      >  2:441168/1‑139 (MQ=255)
  cGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACACCAACTg                                                                                                                                    >  1:28396/1‑139 (MQ=255)
             cgcATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACCCCAACTGCACCAGCTGCa                                                                                                                         >  1:371394/1‑139 (MQ=255)
                 tGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACCCCAACTGCACCAGCTGCATTGt                                                                                                                     >  1:588651/1‑139 (MQ=255)
                  ggCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACCCCAACTGCACCAGCTGCATTGTc                                                                                                                    >  2:498770/1‑139 (MQ=255)
                  ggCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACCCCAACTGCACCAGCTGCATTGTc                                                                                                                    >  2:289427/1‑139 (MQ=255)
                      aaaGCTATGTACGGCGCTTTCTGGCTTGCTTAAGCACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACCCCAACTGCACCAGCTGCATTGTCGCCa                                                                                                                <  2:669961/139‑1 (MQ=255)
                      aaaGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACCCCAACTGCACCAGCTGCATTGTCGCCa                                                                                                                <  1:441168/139‑1 (MQ=255)
                          cTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACCCCAACTGCACCAGCTGCATTGTCGCCACCGc                                                                                                            >  2:256071/1‑139 (MQ=255)
                             tGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTg                                                                                                         >  1:827916/1‑139 (MQ=255)
                               tACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACCCCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCg                                                                                                       <  2:529109/139‑1 (MQ=255)
                                             gcttgcttAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACCCCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCAc                                                                                            <  1:498770/136‑1 (MQ=255)
                                             gcttgcttAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCg                                                                                          <  2:669138/138‑1 (MQ=255)
                                                  cttAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACCCCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGcc                                                                                     <  1:289427/138‑1 (MQ=255)
                                                        aaCCCTTCGCTGACAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATAccc                                                                               >  1:406479/1‑138 (MQ=255)
                                                             ttCGCTGACAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACCCCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACa                                                                                                    >  2:127997/1‑112 (MQ=255)
                                                             ttCGCTGACAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACCCCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACa                                                                                                    <  1:127997/112‑1 (MQ=255)
                                                                                  ccccGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCg                                                                                                       <  2:367517/88‑1 (MQ=255)
                                                                                  ccccGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCg                                                                                                       >  1:367517/1‑88 (MQ=255)
                                                                                   cccGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAg                                                   <  2:406479/139‑1 (MQ=255)
                                                                                     cGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACCCCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCa                                                                                    <  2:507676/104‑1 (MQ=255)
                                                                                     cGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACCCCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCa                                                                                    >  1:507676/1‑104 (MQ=255)
                                                                                       aCGCCCCCAGCAATAACATTGCCTGATTAGCATCAACCATTCCCACCCCAACTGCCCCAGCTGCATTGCCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCAcc                                               <  2:574420/139‑1 (MQ=255)
                                                                                                 cGATAACATTGCCTGATTAGCATCAACCATTCCCACCCCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCg                                     >  2:213949/1‑139 (MQ=255)
                                                                                                    tAACATTGCCTGATTAGCATCAACCATTCCCACCCCAACTGCACAAGCTGCATTGTCGCCACCGCCGCCGCCAACTGCCACCGCAGCCATACCCCACGCTTTCCCAACTTCAGGTAACAAACCACTCATAATTGCGCTg                                  >  2:459355/1‑139 (MQ=255)
                                                                                                      aCATTGCCTGATTAGCATCAACCATTCCCACCCCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGcc                                >  1:225836/1‑139 (MQ=255)
                                                                                                           gCCTGATTAGCATCAACCATTCCCACCCCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCAc                                                                                            >  2:401488/1‑74 (MQ=255)
                                                                                                           gCCTGATTAGCATCAACCATTCCCACCCCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCAc                                                                                            <  1:401488/74‑1 (MQ=255)
                                                                                                           gCCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCg                            >  2:21543/1‑138 (MQ=255)
                                                                                                                 ttAGCATCAACCATTCCCACCCCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATg                     >  1:165901/1‑139 (MQ=255)
                                                                                                                          aCCATTCCCACCCCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCt                               <  2:126304/120‑1 (MQ=255)
                                                                                                                          aCCATTCCCACCCCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCt                               >  1:126304/1‑120 (MQ=255)
                                                                                                                              ttCCCACCCCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGt         >  2:344991/1‑138 (MQ=255)
                                                                                                                              ttCCCACCCCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGt         >  1:439393/1‑138 (MQ=255)
                                                                                                                              ttCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGt         >  2:178670/1‑138 (MQ=255)
                                                                                                                                     acCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACgaga  >  2:440399/1‑138 (MQ=255)
                                                                                                                                     acCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACgaga  >  2:809488/1‑138 (MQ=255)
                                                                                                                                     |                                                                                                                                         
AACGTTGCGGTAGCGCATGGCAAAAGCTATGTACGGCGCTTTCTGGCTTGCTTAAGAACCCTTCGCTGACAGCAAAATAGACCCCCGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGA  >  NZ_CP009273/3721879‑3722149

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