Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,436,6360AC41.3% 77.3 / 11.1 46V274G (GTA→GGA) secYpreprotein translocase subunit SecY
Reads supporting (aligned to +/- strand):  ref base A (22/5);  new base C (0/19);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 1.06e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTC  >  NZ_CP009273/3436503‑3436767
                                                                                                                                     |                                                                                                                                   
cacaAAAAGCGGTCCCCCGGACAGAAAAAACAGCAAATTTGTTGTCCCCCATTCCCACCCAGTCCCGCCCCCAACCCAGGCCGCGATGGTCGCCGGGAACAAAATAATACTGGAAGCAAAGATTGCCGGGATTcccccc                                                                                                                                <  2:671673/139‑1 (MQ=255)
     aaaGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCa                                                                                                                            >  2:369190/1‑138 (MQ=255)
       aGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCatat                                                                                                                         >  2:370196/1‑139 (MQ=255)
           gTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCAc                                                                                                                     >  2:551338/1‑139 (MQ=255)
            ttGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACt                                                                                                                    >  1:300440/1‑139 (MQ=255)
                 cAGGCTGCAAATACAGCGAAATTGTTGTCAGCCATTTCCAACCAGTCCCGCCCCCGAACCATGACGCGAGGGTCGCGGGGAAAAAAATAATACGGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAg                                                                                                               <  1:348454/139‑1 (MQ=255)
                      cgCAAATACAGCAAAATTGTTGCCGGCCAGTTCCAACCAGTCCCCCCCCCAAACCATGACGCGATGGTCGCCGGAAACAAAATAATACTGGAAGCGAAAATTGCCGGATTTACCCCCGCCATATTCACTTTCAGCGGTa                                                                                                          <  1:518374/138‑1 (MQ=255)
                        aaaaTACAGCAAAATTGTGGTCAGCCAGTTCCACCCAGTCCCCCCCCCGAACCATGACGCGATGGTCGCCGGAAACAAAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTa                                                                                                          <  1:253325/136‑1 (MQ=255)
                           aTACAGCGAATTTGTTGCCAGCCAGTTCCAACCAGTCCCGCCCCCGAACCATGCCGCGATGGTCGCCGGGAAAAAAATAATACTGGAAGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAAtgt                                                                                                     <  2:300440/139‑1 (MQ=255)
                               aGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGCAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCAATATTCACTTTCAGCGGGAAATGTGAGc                                                                                                 >  1:381167/1‑139 (MQ=255)
                                gcaaatttgttgTCAGCCAGTTCCAACAAGTCCCCCCCCCGAACAATGACGCGATGGTCGCCGGAAACAGAATAATACTGAAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                <  1:367926/133‑1 (MQ=255)
                                gcaaatttgttgCCAGCCATTCCCACCCAGTCCCGCCCCCGAACCATGCCGCGATGGTCCCCGGAAACAGAATAATCCGGAAAGCGAAAATTGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                <  1:369190/133‑1 (MQ=255)
                                gCGAATTTTTTTCCACCCATTCCCACCCAGTCCCCCCCCCGACCCATGACGCGATGTTCGCCGGAAAAAAAATAATACTGGAAGCGAAGATTGCCGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                <  2:92275/139‑1 (MQ=255)
                                     ttttttgtCGCCCATTCCCACCCAGCCCCCCCCCCAACCCAGGGCGGGATGGCCGCCGGAAAAAAAATAATACTGGAAGCGAAGATTGCCGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg                                                                                           <  1:199877/135‑1 (MQ=255)
                                     tttgttgTCAGCCAGTTCCAACCAGTCCCGCCCCCGAACCAGGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg                                                                                           <  1:903866/138‑1 (MQ=255)
                                     aTTGTTGTCAGCCAGTTCCAACCAGTCCCGCCCCCGACCCATGACGCGATGGTCGCCGGAAACAAAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg                                                                                           <  1:370196/139‑1 (MQ=255)
                                     aTTGTGGCCAGCCAGTTCCAACCTGTCCCCCCCCCAAACCAGGCCGCAATGGTCCCCGGAAACAAAATAATACTGGAAGGGAAGATGGCGGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg                                                                                           <  1:551338/139‑1 (MQ=255)
                                              aGCCATTCCCACCCAGCCCCCCCCCCGAACCAGGACGCGATGGTCGCCGGAAACAAAATAATACTGGAGGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGa                                                                                  <  2:502478/139‑1 (MQ=255)
                                                          ccAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                      >  1:616349/1‑139 (MQ=255)
                                                              tACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCt                                                                   >  1:500404/1‑138 (MQ=255)
                                                              tACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCt                                                                   >  1:348938/1‑138 (MQ=255)
                                                              tACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCt                                                                   >  1:767851/1‑138 (MQ=255)
                                                              tACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCt                                                                   >  1:395566/1‑138 (MQ=255)
                                                                        gacccATGCCGCGATGGTCGCCGGGAACAAAATAATACTGGAAGCGAAGATTGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgcg                                                        <  1:765966/136‑1 (MQ=255)
                                                                        gAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgcg                                                        >  1:282906/1‑139 (MQ=255)
                                                                        gAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgcg                                                        >  2:478218/1‑139 (MQ=255)
                                                                              tGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGtt                                                   >  1:770661/1‑138 (MQ=255)
                                                                                  gcgtgggTCGCCGGGAACAGAATAATCCTGAAAGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACca                                              <  1:820375/134‑1 (MQ=255)
                                                                                      tgttCGCCGGGAACAGAATAATACTGGAAGCTAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGAAAATGTGTGCTCTGTGCAGCATAGCCACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAAt                                          >  2:445572/4‑139 (MQ=255)
                                                                                      ggtccGCCGGAAAAAAAATAATACTGAAAGCAAATATGGCGGGGATCACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAAt                                          <  2:380164/135‑1 (MQ=255)
                                                                                       ggcggCCGGGAAAAGAAAAATACGGGAGGCGAAGGTGGCGGGAATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATg                                         <  2:767851/135‑1 (MQ=255)
                                                                                               ggAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGtt                                  >  1:294322/1‑138 (MQ=255)
                                                                                                aaaCAAAATAATACTGGAAGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGa                                <  2:69978/138‑1 (MQ=255)
                                                                                                  aCAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGAc                               >  1:709384/1‑138 (MQ=255)
                                                                                                      aaaaataCTGGAGGCGAAGGTGTCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGc                          <  2:294322/136‑1 (MQ=255)
                                                                                                           taCTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggc                                     <  1:104359/123‑1 (MQ=255)
                                                                                                           taCTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggc                                     >  2:104359/1‑123 (MQ=255)
                                                                                                           taCTGGAAGCAAAGTTTGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaac                     <  2:256288/139‑1 (MQ=255)
                                                                                                                gAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAAtgt                                                                                                     <  1:400697/54‑1 (MQ=255)
                                                                                                                gAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgt                                                                                                     >  2:400697/1‑54 (MQ=255)
                                                                                                                gAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGtt                                                   >  1:872250/1‑104 (MQ=255)
                                                                                                                gAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                >  1:527661/1‑139 (MQ=255)
                                                                                                                gAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                >  2:615214/1‑139 (MQ=255)
                                                                                                                gAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                >  1:256288/1‑139 (MQ=255)
                                                                                                                gAAGCGAAGATTGCCGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGtt                                                   <  2:872250/104‑1 (MQ=255)
                                                                                                                              ggttaTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTc  <  2:493817/135‑1 (MQ=255)
                                                                                                                              cGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTc  <  2:348938/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                   
CACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTC  >  NZ_CP009273/3436503‑3436767

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