Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,436,6360AC45.6% 78.5 / 22.5 57V274G (GTA→GGA) secYpreprotein translocase subunit SecY
Reads supporting (aligned to +/- strand):  ref base A (25/6);  new base C (0/26);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 7.52e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCA  >  NZ_CP009273/3436502‑3436773
                                                                                                                                       |                                                                                                                                         
acacATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAcccc                                                                                                                                       >  1:251101/1‑139 (MQ=255)
acacATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAcccc                                                                                                                                       >  2:691538/1‑139 (MQ=255)
acacATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAcccc                                                                                                                                       >  2:314006/1‑139 (MQ=255)
acacATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAcccc                                                                                                                                       >  1:435693/1‑139 (MQ=255)
      aaaGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCa                                                                                                                                  >  2:326491/1‑138 (MQ=255)
        aGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCatat                                                                                                                               >  1:599355/1‑139 (MQ=255)
            gTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCAc                                                                                                                           >  1:551540/1‑139 (MQ=255)
            gTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCAc                                                                                                                           >  2:372084/1‑139 (MQ=255)
                cccAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACttt                                                                                                                        >  1:408478/1‑138 (MQ=255)
                  acagCTGAAAATACAGCAAAATTGTTGTCAGCCAGTCCCAACCAGTCCCCCCCCCGACCCAGGCCGCGTTGGCCGCCGGGAACAAAATAATACTGGAAGCAAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAg                                                                                                                     <  2:903521/136‑1 (MQ=255)
                       ggaaaaTACAGCGAATTTGTTGACAGCCGTTCCAACCCAGTCCCCCCCCCGACCCATGCCGCGATGGCCGCCGGAAACAAAATAATACTGGAGGCAAAGATTGCCGGGATCCCCCCCGCCATATTCACTTTCAGCGGTa                                                                                                                <  2:245961/136‑1 (MQ=255)
                       cccAAATACGGCAAAATGGTTGTCAGCCAGTTCCACCAATTACCCCCCCCGAACCATGCCGCGTTGGTCGCCGGGAACAAAATAATATTGGAAGGGAATATTGCCGGAATTCCCCCCGCCATATTCACTTTCAGCGGTa                                                                                                                <  2:371090/137‑1 (MQ=255)
                         cAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTCCCGCCCCCGAACCATGACGCGATGGTCGCGGGAAACAAAATAATACGGGAAGCGAAGTTTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTa                                                                                                                <  1:64191/137‑1 (MQ=255)
                         cAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTaaa                                                                                                              >  1:106865/1‑139 (MQ=255)
                         cAAATACAGCAAAATTGTTTTCAGCCGGTTCCACCCATTCCCCCCCCCGAACCATGACGCGATGGTCGCCGGAAACAAAATAATCCTGGAGGCGAATATTGCCGTGTTTCCCCCCGCCATATTCACTTTCAGCGGTa                                                                                                                <  1:363019/137‑1 (MQ=255)
                         aaaaTCCAGCAAAATTTTTGTCAGCCAGTTCCACCCAGTCCCCCCCCCGAACCATGACGCGATGTCCGCCGGGAACAGAAAAATACTGAAAGCGAAGATTGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTa                                                                                                                <  1:355043/136‑1 (MQ=255)
                         aaaaTACAGAGAATTGGTTGTCAGCCAGTCCAAACAATTCCCCCCCCCGACCCATGCCGGGGTGGTCGCCGGAAAAAAAATAATACGGGAAGGGAAGATTGCCGGTTTTCCCCCCGCCATATTCACTTTCAGCGGTa                                                                                                                <  2:951230/136‑1 (MQ=255)
                           aaTACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtg                                                                                                            >  1:572673/1‑139 (MQ=255)
                               cAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgtgtg                                                                                                        >  1:371090/1‑139 (MQ=255)
                                 gcaataTTGTTGCCAGCCAGTTCCACCCACTCCCCCCCCCGAACAATGCGGCGATGGCCGCGGGGAAAAAAATAATTCTGGAAGCAAAGTTTGCCGGAATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                      <  2:106865/134‑1 (MQ=255)
                                 gcaaaaTTGTTGTCAGCCAGTCCCACCCAGTCCCGCCCCCGAACCATGACGCGAGGGTCGCCGGGAACAAAATAATACTGGAAGCGAAGATGGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                      <  1:327465/136‑1 (MQ=255)
                                 gcaaaaTTGTTGTCAGCCAGTCCCAACCAGTACCCCCCCCGAACCATGACGCGATGGTCGCCGGGAACAAAAAAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                      <  1:691538/136‑1 (MQ=255)
                                 gcaaaaTTGTTGCCAGCCAGTTCCAACCAGTCCCCCCCCCAAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                      <  1:858609/136‑1 (MQ=255)
                                 gCGAATTTGTTGTCACCCATTCCCACCCAGTCCCCCCCCCGAACCATGACGCGATGGTCGCCGGAAACAAAATAATACTGGAGGCGAAGATTGCCGGGTTTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                      <  1:778260/139‑1 (MQ=255)
                                 gCGAAATTGTTGTCATCCAGTCCCAACCAGTCCCCCCCCCAAACCATGACGCGTTGGTCGCCGGGAACAGAAAAATACGGGAAGGGAAGTTTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                      <  2:251101/139‑1 (MQ=255)
                                 gCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgtgctgt                                                                                                      >  1:903521/1‑135 (MQ=255)
                                 gCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                      >  1:296078/1‑139 (MQ=255)
                                 gCGAAATTGTTGCCAGCCAGTTCCAACCAGTCCCGCCCCCGAACCATGACGCGATGGGCGCCGGGAACAGAATAATCCTGGAAGCAAATAGTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                      <  1:922161/139‑1 (MQ=255)
                                 cccaaaTTGTTTACACCCATTACCAACCATTCCCGCCCCCAACCCATGACGCGATGGTCGCCGGAAACAGAAAAATCCTGAAAGCGAAGTTTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                      <  1:540358/136‑1 (MQ=255)
                                      ttttttgcCAGCCATTCCCAACCAGTACCGCCCCCAACCCATGCCGCGATGGCCGCGGGGAACAAAAAAATACTGAAAGCGAAGATGGCGGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg                                                                                                 <  1:741767/135‑1 (MQ=255)
                                      atggttgcCACCCAGTTCCACCCAGACCCCCCCCCAACCCAGGCCGCGATGGTCGCGGGAAACAAAATAATACTGAAAGCGAAGATTGCGGTGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg                                                                                                 <  2:600442/136‑1 (MQ=255)
                                        tgttgcAAGCCATTCCAACCCAGTACCGCCCCCGAACCAGGACGCGATGGCCGCCGGGAAAAAAATAATACGGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgc                                                                                                <  1:656712/138‑1 (MQ=255)
                                               aGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGa                                                                                        >  1:83780/1‑139 (MQ=255)
                                                         cacccAGTCCCGCCCCCGAACCAGGCCGCGATGTTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACAcgacga                                                                                <  2:269371/135‑1 (MQ=255)
                                                                          gacccATGACGCGATGGCCGCCGGGAACAAAATAATCCGGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgcg                                                              <  2:551540/136‑1 (MQ=255)
                                                                          gAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgcg                                                              <  1:219826/139‑1 (MQ=255)
                                                                          gAACCATGACGCGATGGTCGCCGGAAACAGAATAATACTGGAAGCAAAGATGGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgcg                                                              <  1:914106/139‑1 (MQ=255)
                                                                                tGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGtt                                                         >  1:825882/1‑138 (MQ=255)
                                                                                tGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGtt                                                         >  2:101812/1‑138 (MQ=255)
                                                                                tGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGtt                                                         >  2:483037/1‑138 (MQ=255)
                                                                                 gACGCGAGGGCCGCCGGGAACAAAATAATCCTGAAAGGGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTa                                                       <  1:314006/139‑1 (MQ=255)
                                                                                    gcgATGGTCGCGGGGAACAAAATAATCCTGAAAGCAAAGTTGGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACca                                                    <  2:408478/139‑1 (MQ=255)
                                                                                         gcgcGCCGGAAACAAAATAATACGGAAAGCGAAGAGTGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATg                                               <  2:52224/136‑1 (MQ=255)
                                                                                                        aataatACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGc                                >  1:912856/1‑139 (MQ=255)
                                                                                                             taCTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaa                            >  2:575707/1‑138 (MQ=255)
                                                                                                                  gAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                      >  2:33922/1‑59 (MQ=255)
                                                                                                                  gAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                      <  1:33922/59‑1 (MQ=255)
                                                                                                                  gAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                      >  2:767902/1‑139 (MQ=255)
                                                                                                                  gAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGAAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGAGTAGTTTACCACAATGCGGCGTTGACCACGCTAAACAAAt                       >  2:194262/1‑138 (MQ=255)
                                                                                                                   aaGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcg                                                                <  1:167360/96‑1 (MQ=255)
                                                                                                                   aaGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcg                                                                >  2:167360/1‑96 (MQ=255)
                                                                                                                    aGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATac                     >  2:796176/1‑138 (MQ=255)
                                                                                                                       gAAGATTGCCGGATTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaaca                 <  2:104271/139‑1 (MQ=255)
                                                                                                                            ttGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGaa            <  1:576940/139‑1 (MQ=255)
                                                                                                                                ggggATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTc        <  1:50061/138‑1 (MQ=255)
                                                                                                                                cGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTc        <  1:418162/139‑1 (MQ=255)
                                                                                                                                    aTTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTg    <  2:296078/139‑1 (MQ=255)
                                                                                                                                       aCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCa  <  1:575707/138‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                         
ACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCA  >  NZ_CP009273/3436502‑3436773

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