Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,436,6360AC51.8% 37.9 / 46.8 27V274G (GTA→GGA) secYpreprotein translocase subunit SecY
Reads supporting (aligned to +/- strand):  ref base A (13/0);  new base C (0/14);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 4.99e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-06
Rejected as consensus: Frequency below/above cutoff threshold.

TAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCA  >  NZ_CP009273/3436500‑3436773
                                                                                                                                        |                                                                                                                                         
tAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAcc                                                                                                                                         >  1:132611/1‑139 (MQ=255)
              gTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCAc                                                                                                                           >  1:529773/1‑139 (MQ=255)
                tGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACtt                                                                                                                         >  1:165226/1‑139 (MQ=255)
                tGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACtt                                                                                                                         >  1:574404/1‑139 (MQ=255)
                  ccgcgccTGAAAATACAGAGAAATTTGTGTCAGCCTTTTCCACCCAGTCCCCCCCCCAACCCGTGGCGCGTTGGTCCCCGGAAACAAAAAAATTGGGGAAGGAAATTTTGCGGGTTTTCCCCCCGCCATATTCACTTTc                                                                                                                       <  2:22499/133‑1 (MQ=255)
                    cAGGCTGCAAATACAGAGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAGGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAg                                                                                                                     <  1:322149/139‑1 (MQ=255)
                           cAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTaaa                                                                                                              >  1:22499/1‑139 (MQ=255)
                            aaaTACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAt                                                                                                             >  1:34271/1‑139 (MQ=255)
                                   gcaaatttgttgCCAGCCAGTCCAACCCAGTCCCGCCCCCGACCCAGGACGCGATGGTCGCCGGAAACAGAATAATACGGAAAGGGAAGATGGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                      <  1:8153/133‑1 (MQ=255)
                                   gcaaaaTTTTTGTCAGCCAGTCCCACCCAGTCCCGCCCCCGAACCATGCCGCGATGGTCGCCGGAAACAAAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                      <  1:329108/136‑1 (MQ=255)
                                   gcaaaaTTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGAAACAAAAAAAACCTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                      <  1:549085/136‑1 (MQ=255)
                                   gaaaaaTTTTTGTCAGCCAGTCCCACCCAGTACCGCCCCCAAACCATGCCGCGATGGCCGCCGGAAACAGAATAATACTGAAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                      <  1:460397/136‑1 (MQ=255)
                                     gAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgtgctgtct                                                                                                    >  1:458171/1‑133 (MQ=255)
                                                          caacaaGTCCCCCCCCCAAACCATGCCGCGATGGCCGCCGGAAACAAAATAATACTGGAAGCGAAGATTGCCGGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACgacgac                                                                               <  1:349446/139‑1 (MQ=255)
                                                                      cccccGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGtt                                                                   >  1:38141/1‑139 (MQ=255)
                                                                      cccccGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGtt                                                                   >  1:43593/1‑139 (MQ=255)
                                                                           gacccATGCCGCGTGGGCCGCCGGAAAAAAAATAATACGGGAAGCAAAGTTTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgcg                                                              <  2:38141/136‑1 (MQ=255)
                                                                             cccATGCCGCGATGGCCGCGGGAAACAAAATAATTCTGGAAGGGAAGTTTGCCGGGATCCCCCCCGCCATATTCACTTTCAGCGGTAGATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTa                                                            <  2:505274/138‑1 (MQ=255)
                                                                                         tGGCCGCCGGAAACAGAATAATACGGAAAGCGAAGATGGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACGCGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAAt                                                <  1:236332/139‑1 (MQ=255)
                                                                                                        aaataatACTGGAAGCGAAGTTTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACg                                 <  2:90054/138‑1 (MQ=255)
                                                                                                         aataatTCTGGAAGAGAAGAGTGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGc                                <  2:303899/139‑1 (MQ=255)
                                                                                                              gcatGGAAGCGAAGTTGGCCGGGGTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaac                           <  2:111630/136‑1 (MQ=255)
                                                                                                              caCTGGAAGCGAAATTTGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCATAATGCGGCGTTGACCACGCTcaac                           <  2:237181/138‑1 (MQ=255)
                                                                                                                   gAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                      >  2:342808/1‑139 (MQ=255)
                                                                                                                     aGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATac                     >  1:536536/1‑138 (MQ=255)
                                                                                                                     aGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATac                     >  1:3398/1‑138 (MQ=255)
                                                                                                                                      ttACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACATCAAAGAACGTCACTg    >  1:53546/1‑138 (MQ=255)
                                                                                                                                        aCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCa  <  1:342808/138‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                         
TAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCA  >  NZ_CP009273/3436500‑3436773

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