Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,367,341 T→A 8.0% E191V (GAA→GTA)  pspF ← psp operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,367,3410TA 8.0% 141.0 / 3.9 50E191V (GAA→GTA) pspFpsp operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (21/25);  new base A (2/2);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

GTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACC  >  NC_000913/1367193‑1367481
                                                                                                                                                    |                                                                                                                                            
gTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTAtt                                                                                                                                              <  1:128453/149‑1 (MQ=255)
gTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTAtt                                                                                                                                              <  8:55747/149‑1 (MQ=255)
gTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTAtt                                                                                                                                              <  2:252405/149‑1 (MQ=255)
gTGGCGATACACTGAACGGGCCACCACGTTTTTCAATTCACTAATTTTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTACACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTAtt                                                                                                                                              <  8:237267/149‑1 (MQ=255)
 tGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTc                                                                                                                                             <  1:157358/149‑1 (MQ=255)
 tGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCATGATTTCCCGACACATCTGGATGGCAAAGTATTc                                                                                                                                             >  3:132258/1‑149 (MQ=255)
   gCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATTGCAAAGTATTCt                                                                                                                                            <  4:262846/148‑1 (MQ=255)
   gCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTg                                                                                                                                           <  8:92821/149‑1 (MQ=255)
   gCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTCGATTTCCCGACACATCTGGATGGCAAAGTATTCTg                                                                                                                                           <  6:157258/149‑1 (MQ=255)
               aCGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTa                                                                                                                               <  1:174458/149‑1 (MQ=255)
                       ccacGTTTTTCAATTCACGAATATTTCCCGGCCACCGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCACGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCtt                                                                                                                       >  3:241172/1‑149 (MQ=255)
                            tttttCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGACGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCTTCAACATAATGTCGCTTtcgcg                                                                                                                  >  1:121844/1‑149 (MQ=255)
                                tCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGAAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTtcgcgctcg                                                                                                              >  4:242533/1‑149 (MQ=255)
                                    ttCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTtcgcgctcgcg                                                                                                            >  3:253725/1‑147 (MQ=255)
                                    ttCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCa                                                                                                          <  6:97104/149‑1 (MQ=255)
                                     tCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAg                                                                                                         >  5:239779/1‑149 (MQ=255)
                                       aCGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCTGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAgt                                                                                                        <  6:1007/148‑1 (MQ=255)
                                       aCGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAgtg                                                                                                       <  4:253725/149‑1 (MQ=255)
                                              tttCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCa                                                                                                 <  5:284542/148‑1 (MQ=255)
                                                  cgggCCAACGATAATTCAGCAATGTTTCTCTGGCGCGGGCCGTAAACCCCGGGTTCAGAGGCAGCTTGATTTCCCGCCAGATCTGGATGGCAATGTATTCTGCCATCAAGATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGt                                                                                            <  1:160272/147‑1 (MQ=255)
                                                   cGGCAAACGATAATTCAGAAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTGCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTa                                                                                           <  4:241172/149‑1 (MQ=255)
                                                       cAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTacaac                                                                                       <  3:121873/149‑1 (MQ=255)
                                                             tAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCaaa                                                                                 <  7:18297/149‑1 (MQ=255)
                                                                    gCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAgt                                                                          >  1:148756/1‑149 (MQ=255)
                                                                          tttCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAgtcggtc                                                                     >  1:297026/1‑148 (MQ=255)
                                                                             ctctGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCt                                                                                                                                            >  3:161288/1‑74 (MQ=255)
                                                                             ctctGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCt                                                                                                                                            <  4:161288/74‑1 (MQ=255)
                                                                                tGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAgg                                                              >  2:211997/1‑149 (MQ=255)
                                                                                tGGCGCCGTCCGTAATCCCCGGGAACAGACGCAGCCTGATTACCCCCCACATCTGGCATGCAAATTATTCTGCCAACAACAATATGTCGCTTTCGCGCTCGCGCAGTGGGGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGc                                                                 <  2:337151/146‑1 (MQ=255)
                                                                                  gcgcgcTCCGTAAACCCCGGGAACAGAGGCAGCATGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGt                                                             <  4:132258/148‑1 (MQ=255)
                                                                                  gcgcgcTCCGTAAACCCCGGAAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAg                                                               <  3:242533/146‑1 (MQ=255)
                                                                                        tCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATCGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCTCGCTCGCGCAGTGGTGGAAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGCTCAgcgc                                                       >  6:165108/1‑148 (MQ=255)
                                                                                          cGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATACTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGa                                                     >  5:170855/1‑148 (MQ=255)
                                                                                          cGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATACTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGa                                                     >  5:174976/1‑148 (MQ=255)
                                                                                                  ccGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTg                                               >  1:333545/1‑146 (MQ=255)
                                                                                                  ccGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTg                                               <  2:333545/146‑1 (MQ=255)
                                                                                                          agagGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGAc                                    >  7:143388/1‑149 (MQ=255)
                                                                                                          agagGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGAc                                    >  2:297399/1‑149 (MQ=255)
                                                                                                                     gATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCgg                           <  6:239779/147‑1 (MQ=255)
                                                                                                                     gATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCgg                           >  8:119646/1‑147 (MQ=255)
                                                                                                                     gATTTCCCGACACATCTGGAGGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCgg                           <  5:302503/147‑1 (MQ=255)
                                                                                                                           ccGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCg                    <  2:297026/148‑1 (MQ=255)
                                                                                                                                  aTCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCgg                           >  8:159734/1‑134 (MQ=255)
                                                                                                                                  aTCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATtcgt            >  8:10732/1‑149 (MQ=255)
                                                                                                                                       gATGGCAGAATACTTTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGATGTACAACATCAAAAGCCACTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGACATCGGCATTCGTcgcg        >  8:282741/1‑148 (MQ=255)
                                                                                                                                       gATGGCAAAGTATACTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATtcgtc           <  6:170855/145‑1 (MQ=255)
                                                                                                                                       gATGGCAAAGTATACTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATtcgtc           <  6:174976/145‑1 (MQ=255)
                                                                                                                                           gCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGga                         <  2:143112/127‑1 (MQ=255)
                                                                                                                                           gCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGga                         >  1:143112/1‑127 (MQ=255)
                                                                                                                                           gCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATa    <  1:297399/148‑1 (MQ=255)
                                                                                                                                           gCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATa    <  1:211997/148‑1 (MQ=255)
                                                                                                                                             aaaGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGCCCATCGCCGGGAGATCGGCATTCGTCGCGCATAcc  >  2:197757/1‑148 (MQ=255)
                                                                                                                                              aaGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATc                     <  4:208473/128‑1 (MQ=255)
                                                                                                                                              aaGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATc                     >  3:208473/1‑128 (MQ=255)
                                                                                                                                                   ttCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGAc                                    >  8:152605/1‑108 (MQ=255)
                                                                                                                                                   ttCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGAc                                    <  7:152605/108‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
GTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACC  >  NC_000913/1367193‑1367481

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