Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,578,609 T→G 11.2% I40L (ATT→CTT)  gntR ← DNA‑binding transcriptional repressor GntR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,578,6090TG11.2% 98.6 / 5.7 44I40L (ATT→CTT) gntRDNA‑binding transcriptional repressor GntR
Reads supporting (aligned to +/- strand):  ref base T (21/18);  new base G (2/3);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

ATTCCGCGTAATACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACAAAAGAGATA  >  NC_000913/3578473‑3578750
                                                                                                                                        |                                                                                                                                             
aTTCCGCGTAATACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTAGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGcc                                                                                                                                         >  2:254241/1‑143 (MQ=255)
   gcgcgTAATACCTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGCTATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAAGATTGCcgcg                                                                                                                                      <  1:29791/142‑1 (MQ=255)
    cgcgTCCTACTTCCGCGAACCCCTGGCTGGTCAGAGAAGGGAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGt                                                                                                                                     <  1:69433/143‑1 (MQ=255)
    cgcgTAATCCTTCCGCGAAAACCTGGTTGGCGAGAGAAGGTAACCGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGt                                                                                                                                     <  1:51745/143‑1 (MQ=255)
    cgcgGAATACTTCCGCGAAAACCTGGGTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTGAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAAGATTGCcgcgg                                                                                                                                     <  1:499075/143‑2 (MQ=255)
     gcgTAATACTTCCGCGAACACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTa                                                                                                                                    <  1:736164/143‑1 (MQ=255)
      cgTAATACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCAGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTag                                                                                                                                   >  1:603005/1‑143 (MQ=255)
        tAATACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTagag                                                                                                                                 <  1:208983/143‑1 (MQ=255)
          aTACTTCCGCGACCACCTGGTTTGCGCGGCAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATAGCGGGCGCACGATTGGGAATAGAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCg                                                                                                                               <  2:709379/143‑1 (MQ=255)
              ttCCGCGAAAACCTGGTTGGTGAGAGAAGGTAAGAGTTGGTCAATCGCCCGGCTGGTGGCGGTAGAGAGGATATCGGGCGCACGATGGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCCATCTTGCGGCGTAGAGCGACgg                                                                                                                           >  2:497850/1‑143 (MQ=255)
              ttCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACgg                                                                                                                           >  1:88525/1‑143 (MQ=255)
                   cgAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAAcc                                                                                                                      >  1:435033/1‑143 (MQ=255)
                    gAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAACCCt                                                                                                                     >  1:560067/1‑143 (MQ=255)
                      acagcgggttggtGATGGAGCGTAACAGGACGCCAAGCGCTCGGCTGGTCGCGTTAGGGAGCACAGCGGGCGCCCGATTGGGAATCTAGCCCAGTTCATCCATTGCCGCGGCTATCTTGACGCGTAGAGCGACGGAAACCTGc                                                                                                                   <  2:330542/137‑1 (MQ=255)
                          cTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGGTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCg                                                                                                               >  1:355871/1‑143 (MQ=255)
                                 gtGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCg                                                                                                        <  1:325511/143‑1 (MQ=255)
                                 gtGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGAATAGGTCATCAAGAGCGGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCg                                                                                                        >  1:513535/1‑143 (MQ=255)
                                       gaAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGCTGCGTaaaaa                                                                                                  >  2:656151/1‑143 (MQ=255)
                                                   aCGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGt                                                                                      >  1:647910/1‑143 (MQ=255)
                                                         ttCGCCCGGCTGGGGGCGTTAGAGAGGATAGCGGGCGCACGATTGGGAATATAGCCCAGTTCAGCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCGGGGTTGCGTAAAAAACGGCTGACCGTCAtttt                                                                                <  2:428523/142‑1 (MQ=255)
                                                          tCGCCCGGCTGGTGGCGTTAGAGAGGATAGCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTg                                                                               <  1:141957/143‑1 (MQ=255)
                                                                              gagagGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGc                                                           >  2:731776/1‑143 (MQ=255)
                                                                                   gagatCGGGCGCACGATTGGGGATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACAt                                                      <  2:314710/140‑1 (MQ=255)
                                                                                        cGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGa                                                 <  2:521398/143‑1 (MQ=255)
                                                                                         gggCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAACCCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGaa                                                <  1:451242/143‑1 (MQ=255)
                                                                                          ggCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAg                                               >  1:825394/1‑143 (MQ=255)
                                                                                              cACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACg                                           >  1:393354/1‑143 (MQ=255)
                                                                                               aCGATTGGGAATATCGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACgg                                          <  1:806530/143‑1 (MQ=255)
                                                                                               aCGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTGACGCATACACGGTCAGCCACATCCTGAAGTACgg                                          >  1:257591/1‑143 (MQ=255)
                                                                                                cGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACCTCCTGAAGTCCggg                                         >  2:638452/1‑143 (MQ=255)
                                                                                                     gggAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCAACATCCTGGAGTACCGGtcttc                                    >  2:424405/1‑142 (MQ=255)
                                                                                                         atatAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGtcttttctt                                >  2:819950/1‑143 (MQ=255)
                                                                                                          tataGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGtcttttcttt                               >  1:484222/1‑143 (MQ=255)
                                                                                                          tataGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGtcttttcttt                               <  2:88525/143‑1 (MQ=255)
                                                                                                           gtggCCCAGTTCGGCAAGCGGGGCGGCAAGTTTGCCGCGTGGCGCGCCGTAATCCTGCTCCGGGTTGCGTAAAAAACGGCTGACTGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTcttttctttt                              <  2:324971/140‑1 (MQ=255)
                                                                                                               cccAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGCAGTACGGGTCTTTTCTTTTTCAt                          >  1:813215/1‑143 (MQ=255)
                                                                                                                 cAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTCGTTCATCg                        >  2:391920/1‑143 (MQ=255)
                                                                                                                   gTTCATCAAGCGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTc                      <  1:181808/143‑1 (MQ=255)
                                                                                                                    ttCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTACAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATGCTGAAGTACGGGTGTTTTCTTTTTCACCGTcc                     >  2:436081/1‑143 (MQ=255)
                                                                                                                        tcaAGAGCCGCGGGACTCCTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGaa                 <  2:532218/143‑1 (MQ=255)
                                                                                                                        tcaAGAGCCGCGGAAATCGTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGaa                 <  2:325814/143‑1 (MQ=255)
                                                                                                                              tacgcgGCAAGATTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACa           >  2:207471/3‑143 (MQ=255)
                                                                                                                               acgcgGCAAGATTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACaa          >  1:489492/2‑143 (MQ=255)
                                                                                                                                cgcgGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGCTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACaaa         <  1:254916/143‑1 (MQ=255)
                                                                                                                                       gtCTTGCCGCGTAGTTGGACTGAAACCTTCTCCGGGTTGCGTAAGCAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATACTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTGCAAAAGagata  <  2:332493/142‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                             
ATTCCGCGTAATACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACAAAAGAGATA  >  NC_000913/3578473‑3578750

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