Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,303,9280AG69.2% 32.6 / 20.6 26L200S (TTG→TCG) kchvoltage‑gated potassium channel protein
Reads supporting (aligned to +/- strand):  ref base A (5/3);  new base G (8/10);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGCATTGTATGATTGTTTCCTTTTACAAGTTTGTTGAATCCCCCGCGGATAAGCGGACCAAAAATTGAGGTCATAGATGTGGCAAAAACGGTAATGCCGGAAATAATGACCGAAATAGTGAACAATCGTGCTGATTCAGAAACAGGGACAATATCGCCATAGCCGACGGTTGACATGGTTTCTATCGAAAAATAGAACGCGGTCATCAAACTTTCTATTCGCGGATTAAAACCTTCGCTTAAATAAAGCG  >  NZ_CP009273/1303805‑1304054
                                                                                                                           |                                                                                                                              
tGCATTGTATGATTGTTTCCTTTTACAAGTTTGTTGAATCCCCCGCGGATAAGCGGACCAAAAATTGAGGTCATAGATGTGGCAAAAACGGTAATGCCGGAAATAATGACCGAAATAGTGAACGATCGTg                                                                                                                          >  2:784964/1‑130 (MQ=255)
tGCATTGTATGATTGTTTCCTTTTACAAGTTTGTTGAATCCCCCGCGGATAAGCGGACCAAAAATTGAGGTCATAGATGTGGCAAAAACGGTAATGCCGGAAATAATGACCGAAATAGTGAACGATCGTg                                                                                                                          <  1:784964/130‑1 (MQ=255)
tGCATTGTATGATTGTTTCCTTTTACAAGTTTGTTGAATCCCCCGCGGATAAGCGGACCAAAAATTGAGGTCATAGATGTGGCAAAAACGGTAATGCCGGAAATAATGACCGAAATAGTGAACGATCGTGCTGAtt                                                                                                                    >  1:658379/1‑136 (MQ=255)
  cATTGTATGATTGTTTCCTTTTACAAGTTTGTTGAATCCCCCGCGGATAAGCGGACCAAAAATTGAGGTCATAGATGTGGCAAAAACGGTAATGCCGGAAATAATGACCGAAATAGTGAACGATCGTGCTGATTCAGa                                                                                                                >  2:264293/1‑138 (MQ=255)
          gATTGTTTCCTTTTACAAGTTTGTTGAATCCCCCGCGGATAAGCGGACCAAAAATTGAGGTCATAGATGTGGCAAAAACGGTAATGCCGGAAATAATGACCGAAATAGTGAACAATCGTGCTGATTCAGAAACAGGGa                                                                                                        >  1:118800/1‑138 (MQ=255)
               tttCCTTTTACAAGTTTGTTGAATCCCCCGCGGATAAGCGGACCAAAAATTGAGGTCATAGATGTGGCAAAAACGGTAATGCCGGAAATAATGACCGAAATAGTGAACGATCGTGCTGAtt                                                                                                                    <  2:851217/121‑1 (MQ=255)
               tttCCTTTTACAAGTTTGTTGAATCCCCCGCGGATAAGCGGACCAAAAATTGAGGTCATAGATGTGGCAAAAACGGTAATGCCGGAAATAATGACCGAAATAGTGAACGATCGTGCTGAtt                                                                                                                    >  1:851217/1‑121 (MQ=255)
               tttCCTTTTACAAGTTTGTTGAATCCCCCGCGGATAAGCGGACCAAAAATTGAGGTCATAGATGTGGCAAAAACGGTAATGCCGGAAATAATGACCGAAATAGTGAACAATCGTGCTGATTCAGAAACAGGGACAatat                                                                                                  >  1:689603/1‑139 (MQ=255)
                 tCCTTTTACAAGTTTGTTGAATCCCCCGCGGATAAGCGGACCAAAAATTGAGGTCATAGATGTGGCAAAAACGGTAATGCCGGAAATAATGACCGAAATAGTGAACGATCGTGCTGATTCAGAAACAGGGACAATATCg                                                                                                >  2:446239/1‑139 (MQ=255)
                      ttACAAGTTTGTTGAATCCCCCGCGGATAAGCGGACCAAAAATTGAGGTCATAGATGTGGCAAAAACGGTAATGCCGGAAATAATGACCGAAATAGTGAACGATCGTGCTGATTCAGAAACAGGGACAATATCGCCa                                                                                             <  1:665445/137‑1 (MQ=255)
                      ttACAAGTTTGTTGAATCCCCCGCGGATAAGCGGACCAAAAATTGAGGTCATAGATGTGGCAAAAACGGTAATGCCGGAAATAATGACCGAAATAGTGAACGATCGTGCTGATTCAGAAACAGGGACAATATCGCCa                                                                                             <  1:264293/137‑1 (MQ=255)
                                     aTCCCCCGCGGATAAGCGGACCAAAAATTGAGGTCATAGATGTGGCAAAAACGGTAATGCCGGAAATAATGACCGAAATAGTGAACAATCGTGCTGATTCAGAAACAGGGACAATATCGCCATAGCCGACGGTTGACAt                                                                            <  2:689603/139‑1 (MQ=255)
                                                        aCCAAAAATTGAGGTCATAGATGTGGCAAAAACGGTAATGCCGGAAATAATGACCGAAATAGTGAACGATCGTGCTGATTCAGAAACAGGGACAATATCGCCATAGCCGACGGTTGACATGGTTTc                                                                      <  1:107823/126‑1 (MQ=255)
                                                        aCCAAAAATTGAGGTCATAGATGTGGCAAAAACGGTAATGCCGGAAATAATGACCGAAATAGTGAACGATCGTGCTGATTCAGAAACAGGGACAATATCGCCATAGCCGACGGTTGACATGGTTTc                                                                      >  2:107823/1‑126 (MQ=255)
                                                        aCCAAAAATTGAGGTCATAGATGTGGCAAAAACGGTAATGCCGGAAATAATGACCGAAATAGTGAACGATCGTGCTGATTCAGAAACAGGGACAATATCGCCATAGCCGACGGTTGACATGGTTTCTATCGaaaaa                                                            >  1:730333/1‑136 (MQ=255)
                                                                               tgGCAAAAACGGTAATGCCGGAAATAATGACCGAAATAGTGAACGATCGTGCTGATTCAGAAACAGGGACAATATCGCCATAGCCGACGGTTGACATGGTTTCTATCGAAAAATAGAACGCGGTCATCAAACTTTCTAt                                  >  2:214036/1‑139 (MQ=255)
                                                                                ggCAAAAACGGTAATGCCGGAAATAATGACCGAAATAGTGAACAATCGTGCTGATTCAGAAACAGGGACAATATCGCCATAGCCGACGGTTGACATGGTTTCTATCGAAAAATAGAACGCGGTCATCAAACTTTCTAtt                                 >  1:217520/1‑139 (MQ=255)
                                                                                ggCAAAAACGGTAATGCCGGAAATAATGACCGAAATAGTGAACAATCGTGCTGATTCAGAAACAGGGACAATATCGCCATAGCCGACGGTTGACATGGTTTCTATCGAAAAATAGAACGCGGTCATCAAACTTTCTAtt                                 >  2:292587/1‑139 (MQ=255)
                                                                                              tGCCGGAAATAATGACCGAAATAGTGAACGATCGTGCTGATTCAGAAACAGGGACAATATCGCCATAGCCGACGGTTGACATGGTTTCTATCGAAAAATAGAACGCGGTCATCAAACTTTCTATTCGCGGATTAAAAcc                   <  2:730333/139‑1 (MQ=255)
                                                                                              tGCCGGAAATAATGACCGAAATAGTGAACGATCGTGCTGATTCAGAAACAGGGACAATATCGCCATAGCCGACGGTTGACATGGTTTCTATCGAAAAATAGAACGCGGTCATCAAACTTTCTATTCGCGGATTAAAAcc                   <  2:744705/139‑1 (MQ=255)
                                                                                                ccGGAAATAATGACCGAAATAGTGAACGATCGTGCTGATTCAGAAACAGGGACAATATCGCCATAGCCGACGGTTGACATGGTTTCTATCGAAAAATAGAACGCGGTCATCAAACTTTCTATTCGCGGATTAAAACCtt                 <  1:214036/139‑1 (MQ=255)
                                                                                                   gAAATAATGACCGAAATAGTGAACGATCGTGCTGATTCAGAAACAGGGACAATATCGCCATAGCCGACGGTTGACATGGTTTCTATCGAAAAATAGAACGCGGTCATCAAACTTTCTATTCGCGGATTAAAACCTTCGc              <  1:812492/139‑1 (MQ=255)
                                                                                                              cGAAATAGTGAACAATCGTGCTGATTCAGAAACAGGGACAATATCGCCATAGCCGACGGTTGACATGGTTTCTATCGAAAAATAGAACGCGGTCATCAAACTTTCTATTCGCGGATTAAAACCTTCGCTTAAATAAAgc   <  2:217520/139‑1 (MQ=255)
                                                                                                               gAAATAGTGAACGATCGTGCTGATTCAGAAACAGGGACAATATCGCCATAGCCGACGGTTGACATGGTTTCTATCGAAAAATAGAACGCGGTCATCAAACTTTCTATTCGCGGATTAAAACCTTCGCTTAAATAAAgcg  <  2:765360/139‑1 (MQ=255)
                                                                                                                      tGAACAATCGTGCTGATTCAGAAACAGGGACAATATCGCCATAGCCGACGGTTGACATGGTTTCTATCGAAAAATAGAACGCGGTCATCAAACTTTCTATTcg                               >  2:703457/1‑103 (MQ=255)
                                                                                                                      tGAACAATCGTGCTGATTCAGAAACAGGGACAATATCGCCATAGCCGACGGTTGACATGGTTTCTATCGAAAAATAGAACGCGGTCATCAAACTTTCTATTcg                               <  1:703457/103‑1 (MQ=255)
                                                                                                                           |                                                                                                                              
TGCATTGTATGATTGTTTCCTTTTACAAGTTTGTTGAATCCCCCGCGGATAAGCGGACCAAAAATTGAGGTCATAGATGTGGCAAAAACGGTAATGCCGGAAATAATGACCGAAATAGTGAACAATCGTGCTGATTCAGAAACAGGGACAATATCGCCATAGCCGACGGTTGACATGGTTTCTATCGAAAAATAGAACGCGGTCATCAAACTTTCTATTCGCGGATTAAAACCTTCGCTTAAATAAAGCG  >  NZ_CP009273/1303805‑1304054

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