Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,367,353 A→G 7.7% M187T (ATG→ACG)  pspF ← DNA‑binding transcriptional dual regulator PspF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,367,3530AG 7.7% 136.7 / 3.3 53M187T (ATG→ACG) pspFDNA‑binding transcriptional dual regulator PspF
Reads supporting (aligned to +/- strand):  ref base A (22/26);  new base G (2/2);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01

CCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATT  >  NC_000913/1367216‑1367492
                                                                                                                                         |                                                                                                                                           
ccacGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATg                                                                                                                                        >  1:884069/1‑143 (MQ=255)
 cacGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGt                                                                                                                                       >  1:1014535/1‑143 (MQ=255)
 cacGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGt                                                                                                                                       >  1:1061773/1‑143 (MQ=255)
 cacGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGt                                                                                                                                       >  2:51569/1‑143 (MQ=255)
  acGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTc                                                                                                                                      >  1:1007492/1‑143 (MQ=255)
        ttCAATTCACTAATATTGCCCGGCCAACGGTAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTTATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTtc                                                                                                                                <  2:688013/143‑1 (MQ=255)
           aaTTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGAAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCCGCCACCAACAGTCTGCCGCTTtcgcg                                                                                                                             >  2:852899/1‑143 (MQ=255)
              tCACGAATATTTCCCGGCCAACGGTACTGCAGCAATGTTTATCTGGCGCGCTCCGCAAACCCCGGGAACCGATGCATTTGCATTTCCCTACACATCTGGTTGCCCTAGTATTCTGCCATCAACATGAGTTCTCTTtcgccccc                                                                                                                          >  2:897189/1‑139 (MQ=255)
                  gtgtatTGCCCGGCCAACGACAGGTCAGCAATGCTTCTCTGGCGGGCTCCGTAAACCCCGGGAACAGAGGCAGCTGGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTtcgcgctcgcgc                                                                                                                      <  1:1072097/140‑1 (MQ=255)
                  gAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTtcgcgcccgcgc                                                                                                                      >  1:867804/1‑143 (MQ=255)
                  gAATATGTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTtcgcgctcgcgc                                                                                                                      <  1:177458/143‑1 (MQ=255)
                   aaTATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCa                                                                                                                     >  1:450084/1‑143 (MQ=255)
                     tatTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAgt                                                                                                                   <  2:786713/143‑1 (MQ=255)
                     tatTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAgt                                                                                                                   <  2:329696/143‑1 (MQ=255)
                             ggCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAg                                                                                                           <  1:1105181/143‑1 (MQ=255)
                               ccAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGtt                                                                                                         <  2:466779/143‑1 (MQ=255)
                               ccAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGtt                                                                                                         <  1:887896/143‑1 (MQ=255)
                                 aaCGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGACTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATCATGTCGCCTTCGCGCTCGCGCCGCGCCGCCAGTTGt                                                                                                       >  2:255235/1‑143 (MQ=255)
                                      gaaTGCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTacaac                                                                                                  <  1:107190/142‑1 (MQ=255)
                                             gaaaTGTTACTCTTGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCaaaa                                                                                           <  1:405013/141‑1 (MQ=255)
                                               aaTGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCGTTCGCGCTCGCGCAGTGGTGGCAGATGAACAACCTCacacgc                                                                                         >  2:94378/1‑138 (MQ=255)
                                                  gTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAg                                                                                      <  1:623478/143‑1 (MQ=255)
                                                         tGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAgtcggtcg                                                                               >  2:205262/1‑143 (MQ=255)
                                                         cggCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACGTTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAgtcggtcg                                                                               >  2:205267/2‑143 (MQ=255)
                                                         cggCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACGTTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAgtcggtcg                                                                               >  2:205263/2‑143 (MQ=255)
                                                          ggCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGGTGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGa                                                                              >  2:1056232/1‑143 (MQ=255)
                                                                cTCCGTAAACCCCGGGAACAGACGCAGCTTGCTTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGt                                                                        <  2:860291/143‑1 (MQ=255)
                                                                cTCCGGAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGt                                                                        <  2:857514/143‑1 (MQ=255)
                                                                   cGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGTCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACCTCAAAAGCCAGTCTGCCGAGCAggaaag                                                                     >  2:769565/1‑139 (MQ=255)
                                                                           ccGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGaaaa                                                             <  1:154003/143‑1 (MQ=255)
                                                                               gAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTCGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGc                                                         <  1:576967/143‑1 (MQ=255)
                                                                                   agagGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTc                                                     <  1:1038118/143‑1 (MQ=255)
                                                                                   agagGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGGGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCCGCGCGAAAAGTGCCTTc                                                     >  1:576668/1‑143 (MQ=255)
                                                                                   agagGCAGCGTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTc                                                     <  1:400686/143‑1 (MQ=255)
                                                                                      ggCAGCTTTATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATTAATATTATGTCGCTTTCGCGCTCGCTCAGTGGTGGCAGTTGTACAACAGCAAAAGGCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCAtt                                                  <  1:331503/143‑1 (MQ=255)
                                                                                      ggCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCAtt                                                  <  1:861396/143‑1 (MQ=255)
                                                                                      ggCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCAtt                                                  <  2:884069/143‑1 (MQ=255)
                                                                                       gCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCCGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTg                                                 >  1:274241/1‑143 (MQ=255)
                                                                                               atgtCCCGACACAGGTGCGTGGCAAAGCATTCGGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGc                                         <  2:620887/140‑1 (MQ=255)
                                                                                               aTTTCCCGACATATCTGGATGGCAAAGTATTCTGTCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGc                                         <  2:48484/143‑1 (MQ=255)
                                                                                                    ccGACACATCTGGAGGGCATAGTATTCTGCCATCAACCTTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGTAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGga                                    <  1:671107/143‑1 (MQ=255)
                                                                                                         acaTCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCCGGTCAGCGCGAAAAGTGCCTTCATGGACCATCACCggggagtcg                               >  2:377186/1‑137 (MQ=255)
                                                                                                            tCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCGTCATTGACCATCGCCGGGAGATCGGcg                            >  1:705542/1‑142 (MQ=255)
                                                                                                              tGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCCGGAGATCGGCAtt                          >  1:1015301/1‑143 (MQ=255)
                                                                                                                gATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATtcc                        >  2:389108/1‑142 (MQ=255)
                                                                                                                    gCAAAGTATTCTGCCATCAACGTTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTcgc                    <  1:205263/143‑1 (MQ=255)
                                                                                                                    gCAAAGTATTCTGCCATCAACGTTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTcgc                    <  1:205267/143‑1 (MQ=255)
                                                                                                                    gCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTcgc                    <  1:205262/143‑1 (MQ=255)
                                                                                                                     aaaaGTATTCTGCTTTCAACATTCTGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTcgcg                   <  1:852899/142‑1 (MQ=255)
                                                                                                                      aaaGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGCTCGGCATTCGTCgcgc                  >  1:239614/1‑143 (MQ=255)
                                                                                                                        aGTATGCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGCTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCAt                <  2:1020299/143‑1 (MQ=255)
                                                                                                                                 ccATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGc       <  1:1056232/143‑1 (MQ=255)
                                                                                                                                      aaCATTATGTCGCTTTCGCGCTCGCGCAGAGGTGGCAGTTGTACAACATCAAAAGCAAGTCGATCGAGCAGATCAGCGCGACAAGTGCATTCATCGACCATCGCCGGGCGATCCGCCTTCGTCGCGCATACCAACCGCACCtt  >  1:180023/1‑143 (MQ=255)
                                                                                                                                         |                                                                                                                                           
CCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATT  >  NC_000913/1367216‑1367492

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