Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,784,808 C→T 7.0% G769R (GGG→AGG)  ppsA ← phosphoenolpyruvate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,784,8080CT 7.0% 157.1 / 3.0 57G769R (GGG→AGG) ppsAphosphoenolpyruvate synthase
Reads supporting (aligned to +/- strand):  ref base C (18/35);  new base T (2/2);  total (20/37)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

TGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGCCATCGACAGCAGTGCTTTCACCGCATCGTTGCGCTCATCGAACAATTCA  >  NC_000913/1784662‑1784953
                                                                                                                                                  |                                                                                                                                                 
tgAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCct                                                                                                                                                 <  2:156303/149‑1 (MQ=255)
tgAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTATTCAGTTCAGCAAGGCTTAACCAGGATTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCct                                                                                                                                                 >  2:334630/1‑149 (MQ=255)
   agagagTGAATAAGTGGCATCTTCGGGGAGCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCCGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCtt                                                                                                                                              <  1:692474/149‑1 (MQ=255)
   agagATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACACGCTATCGATCCCCTCtt                                                                                                                                              >  2:826171/1‑149 (MQ=255)
           aaTAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAAc                                                                                                                                      <  2:501346/149‑1 (MQ=255)
           aaTAACGTTCATCTTCGGGGCTCAAATAACCCCTGCGACTAAACGCCGCCGGGTAGTTCTATTATTTCTTCAGTTCAGCCATGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAAc                                                                                                                                      <  2:669369/149‑1 (MQ=255)
              aaGTTTCATCTTCGGGGATCACATACGCCCGGAGACTCAACGACGCCGGGGATTTATTTTATTTCTGCAGTTCAGCGAGGCGTAACCAGGTAGGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCAGCAACCAt                                                                                                                                   <  1:666723/149‑1 (MQ=255)
                 tttCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCg                                                                                                                                >  1:1022453/1‑149 (MQ=255)
                      tCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCaaa                                                                                                                           <  1:998425/149‑1 (MQ=255)
                         tCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTc                                                                                                                        <  1:413494/149‑1 (MQ=255)
                          ctgggCTCACATATCCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTGCAGCCAGGCTTAACCAGGTTTGCTCCACGGTGTCCGGGGTCAGAGACACGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCt                                                                                                                       <  2:781530/147‑1 (MQ=255)
                          cGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCt                                                                                                                       <  1:552904/149‑1 (MQ=255)
                            gggATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTc                                                                                                                     <  1:189153/149‑1 (MQ=255)
                            gggATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTc                                                                                                                     >  2:423407/1‑149 (MQ=255)
                                 cacaTAACCCCGGCGACTAATCGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGt                                                                                                                <  1:1016640/149‑1 (MQ=255)
                                     tAACCCCGGCGACTAAACGCCGCCGTGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTcgga                                                                                                            <  2:500826/149‑1 (MQ=255)
                                           cGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGAcc                                                                                                      >  2:925574/1‑149 (MQ=255)
                                            ggCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTACCCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGAccc                                                                                                     <  1:667194/149‑1 (MQ=255)
                                             gCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCt                                                                                                    <  2:327657/149‑1 (MQ=255)
                                             gCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGGTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCt                                                                                                    <  1:124704/149‑1 (MQ=255)
                                                 cTAATCGCCGCCGGGGATTTATTATAATTCTTCAGATCAGGCTGGCTTAACCAGGTTTGCACGACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGActgaccctgacc                                                                                                <  1:886984/149‑1 (MQ=255)
                                                 cTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCCCGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGAcc                                                                                                <  1:784953/149‑1 (MQ=255)
                                                 cTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGAcc                                                                                                <  1:891190/149‑1 (MQ=255)
                                                    aaCGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTCACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGCCTTTCGATCCCCTCTTGCATCAACGATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCgaa                                                                                             >  1:947644/1‑147 (MQ=255)
                                                         cgccgGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATcc                                                                                        <  2:267907/149‑1 (MQ=255)
                                                         cgccgGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATcc                                                                                        <  1:322170/149‑1 (MQ=255)
                                                          gccgGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTGGCACCACGGTGTCCGGGTTCAGATACAGGCTATCGATCCCCTCTTCCTTCAACCATGCGGAAAAGTCTTAGTGTTCGGACGGACACTGCCCGCAAATccc                                                                                       >  1:561503/1‑149 (MQ=255)
                                                                     attttatttCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGc                                                                            >  1:295323/1‑149 (MQ=255)
                                                                      ttttatttCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCTCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCGGACCGCAACTCCCGACATATTTGcc                                                                           <  1:699271/149‑1 (MQ=255)
                                                                      ttttatttCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGCCCGCAAATCCCGACATATTTGcc                                                                           >  2:611656/1‑149 (MQ=255)
                                                                      ttttatttCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCTCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGcc                                                                           <  1:697306/149‑1 (MQ=255)
                                                                        ttatttCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGCCAGGCTATCGATCCTCTCTTCGCTCAACCATGCGGCAATGTCTTCGTGGTCGGACGGACACTGCCCGCCAATCCCGCCAAATTTGccca                                                                         >  2:699271/1‑148 (MQ=255)
                                                                        ttatttCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCTCTCTTCCATCAACCATTCGGCAAAGTCTTAGTCGTCTGACGGACCATTACCGCATATCACTACATATTTGCCCt                                                                         >  2:697306/1‑149 (MQ=255)
                                                                          atttCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTTTACGGGTTCAGAGACAGGCTATCGATACCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGTTCGGACTGACCCTGCCCGAACATCCCTACATATTTGCACTGt                                                                       >  1:697620/1‑149 (MQ=255)
                                                                             tcttcAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTttc                                                                    <  2:45176/149‑1 (MQ=255)
                                                                                 cagttcagCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCCTGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTcg                                                                <  2:750387/149‑1 (MQ=255)
                                                                                  atttcagCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTcgc                                                               <  1:315117/147‑1 (MQ=255)
                                                                                   gttcagCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACTGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTcgcg                                                              >  2:813041/1‑149 (MQ=255)
                                                                                    tgcagCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCgg                                                             <  2:996809/147‑1 (MQ=255)
                                                                                      cagCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCACAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCa                                                           <  1:901644/149‑1 (MQ=255)
                                                                                         ccAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGTAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACgg                                                        <  1:841703/149‑1 (MQ=255)
                                                                                           aGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGAt                                                      >  2:544064/1‑149 (MQ=255)
                                                                                            ggCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATa                                                     <  1:826171/149‑1 (MQ=255)
                                                                                               ttAACCAGGTTTGCACCACGGTGTCCGGGTGCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGcc                                                  <  1:822720/149‑1 (MQ=255)
                                                                                                  aCCCGGATTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCACCCATGCGGCAAAGTCTTCGTGGTCCGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGCCATc                                               <  2:591008/149‑1 (MQ=255)
                                                                                                       gTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGCCATCGAcag                                          >  2:1027202/1‑149 (MQ=255)
                                                                                                                   gtgtCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGCCATCGACAGCAGTGCTTTCAc                              <  1:726713/149‑1 (MQ=255)
                                                                                                                   gtgtCCGGGTTCAGAGACAGGCTAGCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGCCATCGACAGCAGTGCTTTCAc                              <  1:240036/149‑1 (MQ=255)
                                                                                                                            ttCAGAGACAGGCTATCGATCCCCTCTTCCATCACCCATGCGGCAAAGTCTTCGTGCTCGGACGGACCCTGACCGCAAATCCCGACATATATGCCCTGTTTCTTCGCGGCACGGAAAGCCAACGACAGCAGTGCTTTCACCGCATCGtt                     >  1:559710/1‑149 (MQ=255)
                                                                                                                            ttCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGGTTCTTCGCGGCACGGATAGCCATCGACAGCAGTGCTTTCACCGCATCGtt                     >  1:41719/1‑149 (MQ=255)
                                                                                                                               agagACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGCCATCGACAGCAGTGCTTTCACCGCATCGTTgcg                  <  1:127802/149‑1 (MQ=255)
                                                                                                                               agagACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGCCATCGACAGCAGTGCTTTCACCGCATCGTTgcg                  >  1:845628/1‑149 (MQ=255)
                                                                                                                                        gtATCGATCCCCTCTTCCATCAACCATGCGGCATTGTCTTCGTGGTCGGACGGACCCTGACGGCAAATCCCGTCATATTTGCCCTGTTTCTTCGCGGCACGGATAGCCATCGACAGCAGTGCTTTCACCGCATCGTTGCGCTCATCGaa         <  2:12667/148‑1 (MQ=255)
                                                                                                                                        ctttCGATCCCCTCTTCCATCAACCATGCGGCTAAGTCTTAGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGCCATCTACAGCAGTGATTTCACCGCATCGTTGCGCGCTTCGaa         >  1:692025/4‑149 (MQ=255)
                                                                                                                                              aTCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGCCATCGACAGCAGTGCTTTCACCGCATCGTTGCGCTCATCGAACAAtt    <  2:491174/148‑1 (MQ=255)
                                                                                                                                              aTCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGCCATCGACAGCAGTGCTTTCACCGCATCGTTGCGCTCATCGAACAAtt    >  1:491174/1‑148 (MQ=255)
                                                                                                                                               tCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGCCATCGACAGCAGTGCTTTCACCGCATCGTTGCGCTCATCGAACAATTCa  <  2:553615/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGCCATCGACAGCAGTGCTTTCACCGCATCGTTGCGCTCATCGAACAATTCA  >  NC_000913/1784662‑1784953

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