Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,443,2770AC50.0% 67.6 / 17.3 54V274G (GTA→GGA) secYpreprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (24/3);  new base C (0/27);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 4.17e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAAT  >  NC_000913/3443132‑3443423
                                                                                                                                                 |                                                                                                                                                  
cGCATAGAGTAACACATAAAGCGGTTGCCCATGCTGCAAATGCAGCGAAATTGGTGTCAGCCAGTTCCAACCAGTACCGCCCCCGACTCATGAGGCGATGGTCGCCGGGAACAGAATAATTCTGGAAGCGAAGATTGCCGGGATTAccccc                                                                                                                                               <  2:753071/151‑1 (MQ=255)
  cATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGc                                                                                                                                             <  2:768359/151‑1 (MQ=255)
   aTAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGcc                                                                                                                                            >  1:612131/1‑151 (MQ=255)
   aTAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGcc                                                                                                                                            >  1:388404/1‑151 (MQ=255)
                 aagggggTTCCCCAGGCTGAAAATACAGGGAATTTTTTGCCACCTATTCCCAACCAGTCCCCCCCCCAACCCAGGACGGGTTGGCCGCCGGAAAAAAAATAATACTGGAGGCGAAATTGGCCGGAATTCCCCCCGCCATATTCACTTTCa                                                                                                                               <  2:343774/145‑1 (MQ=255)
                     cGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGt                                                                                                                          >  1:411557/1‑151 (MQ=255)
                      gtttGCCCAGGCTAAAAATACAGCAAATTTGTTGTCAGCCATTCCCAACCAGTCCCCCCCCCAAACCATGCCGCGATGGTCGCCGGAAACAAAATAATCCTGGAAGCAAAGATTGCGGGTATTCCCCCCGCCATATTCACTTTCAGCGGt                                                                                                                          <  1:539646/148‑1 (MQ=255)
                          gCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtg                                                                                                                     >  2:352792/1‑151 (MQ=255)
                              aGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgtgtg                                                                                                                 >  1:675906/1‑151 (MQ=255)
                                gCTGAAAATACAGCAAAATTGTTGTCAGCCATTTCCACCCAGTCCCGCCCCCGAACCATGACGCGATGGTCGCCGGAAACAAAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTATGct                                                                                                               <  1:761530/151‑1 (MQ=255)
                                  tgaaaaTACAGCAAATTTGTTGTCAGCCGGTCCCAACCAGTCCCGCCCCCGAACCATGCCGCGATGGCCGCCGGGAAAAAAATAATACTGGAAGCGAAGATTGCGGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                               <  1:314169/146‑1 (MQ=255)
                                  tGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctct                                                                                                             >  1:180294/1‑151 (MQ=255)
                                  tGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctct                                                                                                             >  1:47577/1‑151 (MQ=255)
                                     aaaaACGGCGAAGTTTTGGTCAGCAAGTACAACCAAGCCCCCCCCCCAACCCATGACGCGTTGGTGGCGGGGAACAGAATAATACGGGAAGGGAAGATTGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgt                                                                                                           <  2:808769/150‑1 (MQ=255)
                                      aatacgggaaattttttgtCGCCCATTTCCACCCATTCCCGCCCCCAAACCATGACGGGATGGTCCCCGGAACAAAAATATTACGGGAGCCGAAGTTTGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgc                                                                                                         <  2:148767/136‑1 (MQ=255)
                                            gCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                    >  1:918337/1‑150 (MQ=255)
                                             cGAAATGTTTGACACCCTGTTCCAACCGGACCCGCCCCCCGTCCACGACGCAAGGGCCGCCGGGAAAATAATAATACTGAAGGCAAAGATTGCCGGGAAACCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATa                                                                                                   <  2:545764/150‑1 (MQ=255)
                                              gAAATGGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGAAAATGTGCGCTCTGTgcagca                                                                                                     >  2:915125/1‑147 (MQ=255)
                                                     ttgtCAGCCAGTTCCACCCAGTCCCCCCCCCGAACCAGGCCGCGATGGTCGCCGGAAAAAAAATAATACTGGAAGCGAAGATGGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacga                                                                                            <  1:312677/149‑1 (MQ=255)
                                                     ctgtCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGCGCTCTGTGCAGCATAGACacgacg                                                                                          >  2:649066/2‑151 (MQ=255)
                                                         cAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCt                                                                                      >  1:664273/1‑151 (MQ=255)
                                                                    acccAGTCCCCCCCCCAAACCATGAGGGGATGGTCGCGGGAAAAAAAAAAATCCTGGAAGCAAAGTTTGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGt                                                                             <  2:899171/147‑1 (MQ=255)
                                                                        tgccccGCCCCCGACCCGTGACGCGAGGGCCGAGGGAAAAAAAATAATACTGGAAGGAAAAATGGCGGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcg                                                                         <  1:410774/145‑1 (MQ=255)
                                                                          tACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt                                                                      >  1:72549/1‑150 (MQ=255)
                                                                           aCCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt                                                                      >  1:870473/1‑149 (MQ=255)
                                                                              gCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGttt                                                                 >  2:667854/1‑151 (MQ=255)
                                                                                    gAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcac                                                            >  2:871241/1‑150 (MQ=255)
                                                                                     aaCCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcaca                                                           >  1:89899/1‑150 (MQ=255)
                                                                                       ccATGACGGGATGGTCGCCGGAAAAAAAATAATCCTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATg                                                        <  1:273375/151‑1 (MQ=255)
                                                                                       ccAGGACGCGATGGTCGCCGGGAAAAGAATAATACTGAAAGCGAAGTTTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATg                                                        <  1:257744/151‑1 (MQ=255)
                                                                                            aCGCGATGGTCGCCGGAAACAAAAAAATACTGAAAGCGAATATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcg                                                      <  1:957790/148‑1 (MQ=255)
                                                                                             cgcgtgggTCGCCGGGAACAAAATAATACTGGAAGCGAAGATTGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggc                                                    <  1:848175/143‑1 (MQ=255)
                                                                                                gATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTg                                                >  2:526510/1‑150 (MQ=255)
                                                                                                 aTGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGAc                                              >  1:23621/1‑151 (MQ=255)
                                                                                                      cGCCGGAAACAAAATAATAGTGGAAGCAAAGTTGGCCGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACg                                          <  1:139066/150‑1 (MQ=255)
                                                                                                        ccgaaaaCAAAATAATCCTGAAGGCGAATATTGCCGTGTTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCCCGc                                         <  2:592075/144‑1 (MQ=255)
                                                                                                                aGAATAATCCGGGAAGCGAAGTTTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                               <  2:634024/151‑1 (MQ=255)
                                                                                                                    taataCTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaac                           >  1:733370/1‑151 (MQ=255)
                                                                                                                       tacggaaaGCGAAGATTGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                         <  1:369134/144‑1 (MQ=255)
                                                                                                                        aCTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                         >  1:594943/1‑149 (MQ=255)
                                                                                                                           ggaaggaaAGAGTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGaa                     <  2:706197/150‑1 (MQ=255)
                                                                                                                           ggAAGCGAAGTTTGCGGGCATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGa                      <  2:951040/149‑1 (MQ=255)
                                                                                                                           ggAAGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAg                       <  1:150498/148‑1 (MQ=255)
                                                                                                                           gaaaGCAAAGTTTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACaaa                        <  1:403224/145‑1 (MQ=255)
                                                                                                                             aggCGAAGATGGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACg                   <  2:55019/148‑1 (MQ=255)
                                                                                                                             aaGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGt                  >  1:629514/1‑151 (MQ=255)
                                                                                                                                  aaaaTTGCCGGGTTGCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTg             <  2:388224/148‑1 (MQ=255)
                                                                                                                                  aaGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTg             >  2:611936/1‑151 (MQ=255)
                                                                                                                                     aTTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCa           >  2:839090/1‑150 (MQ=255)
                                                                                                                                      ttcccGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCa           <  2:503959/146‑1 (MQ=255)
                                                                                                                                        gCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCa           <  2:787494/147‑1 (MQ=255)
                                                                                                                                           ggtaTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatac      <  1:267976/146‑1 (MQ=255)
                                                                                                                                           ggcaTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaat        <  2:702589/144‑1 (MQ=255)
                                                                                                                                              aTTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGACGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaat  >  1:568503/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAAT  >  NC_000913/3443132‑3443423

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: