Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,578,609 T→G 8.3% 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,6090TG 8.3% 127.2 / 4.7 55I40L (ATT→CTT) gntRDNA‑binding transcriptional repressor GntR
Reads supporting (aligned to +/- strand):  ref base T (32/18);  new base G (2/3);  total (34/21)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-02

GATTCCGCGTAATACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACAAAAGAGATA  >  NC_000913/3578472‑3578750
                                                                                                                                         |                                                                                                                                             
gactCCGGGTAACACTTCCGCGAAAACCTGGTTGGTCAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGc                                                                                                                                          <  1:112009/140‑1 (MQ=255)
     cgcgTGATACTGCCGCGAAACCCGGGTGGGGGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAAGATTGCcgcgg                                                                                                                                     <  2:86084/143‑2 (MQ=255)
      gcgTAATACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGCATATAGCCGCGTTCATCAAGAGCCGCGGCAATCTTGCCGCGtc                                                                                                                                    >  2:820219/1‑142 (MQ=255)
      gcgTAATACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTa                                                                                                                                    <  1:529661/143‑1 (MQ=255)
          caTACTTCCGCGAACACCGGCTTGGCGAGAGAACGTAACAGGACGACAATCGCCCGGCTGGTGGCGTTAGAGAGGACATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAAGATTGCcgcggctctt                                                                                                                                <  1:331390/142‑7 (MQ=255)
          aaTACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGc                                                                                                                                <  1:56859/143‑1 (MQ=255)
               ttCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCCCGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCACTATTGCCGCGTAGAGCGACgg                                                                                                                           >  2:632248/1‑143 (MQ=255)
               ttCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACgg                                                                                                                           >  1:472268/1‑143 (MQ=255)
                        aaCCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCt                                                                                                                  >  1:518604/1‑143 (MQ=255)
                           cTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGCGCGACGGAACCCTGCTCCg                                                                                                               >  2:316059/1‑143 (MQ=255)
                                    gagagaAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTa                                                                                                      >  2:288731/1‑143 (MQ=255)
                                     agagaAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTaa                                                                                                     >  2:873636/1‑143 (MQ=255)
                                     agagaAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTaa                                                                                                     >  2:64923/1‑143 (MQ=255)
                                      gagaACGGAACAGGACGCCAATCGCCCGGCTGCTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTaaa                                                                                                    <  1:126979/143‑1 (MQ=255)
                                         aaGGTAACAGGACGCGACTCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTaaaaaa                                                                                                 <  2:334758/143‑1 (MQ=255)
                                         aaGCTACCAGGCAGCCAATCGCCCGGCTGGTCGCGTTAGAGAGAATATCGGGCGCACGATTGTGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTaaaaaa                                                                                                 <  2:823525/143‑1 (MQ=255)
                                              aaCAGGGCGCCAATGGCCCGGCTGGGGGCGTGAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTGGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCt                                                                                            <  1:617419/143‑1 (MQ=255)
                                              aaCAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGGTGCGTAAAAAACGGCt                                                                                            >  1:562124/1‑143 (MQ=255)
                                                cAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGa                                                                                          <  1:513812/143‑1 (MQ=255)
                                                cAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTGCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGGAACCCGCCCCCGGTCGCGTAACAAActgacgg                                                                                          >  2:521643/1‑137 (MQ=255)
                                                 aGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGAc                                                                                         >  2:831517/1‑143 (MQ=255)
                                                         gaTCGCCCGGGGTGGGGGGTTACGGCGGCTAGCGGGCCCGCGATTGGGAATAGAGGCGAGTGGGTCCGGAGCCGCGGCAAGCTTGCCGCGTAGCGCGACGGAAACCTGCTCCGGGTCGCGTAAAAAACGGCTGACCGTCAttt                                                                                 <  2:223159/142‑1 (MQ=255)
                                                         aaTCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCAttt                                                                                 >  2:158149/1‑143 (MQ=255)
                                                         aaTCG‑CCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGCTTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGGTGCGTAAAGAACGGCTGCCCGTCAtttt                                                                                >  1:145058/1‑143 (MQ=255)
                                                           tCGCCGGGGTGGGGCCGTTAGAGAGGTGAGCGGGCGCGCGATTGGGAAGATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTg                                                                               <  2:834744/143‑1 (MQ=255)
                                                           tCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGCCTGACCGTCATTTTg                                                                               >  1:728197/1‑143 (MQ=255)
                                                              cccGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTc                                                                            >  1:513047/1‑143 (MQ=255)
                                                               ccGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCa                                                                           >  1:416518/1‑143 (MQ=255)
                                                                      gtggCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTac                                                                    <  2:610789/143‑1 (MQ=255)
                                                                         gCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACg                                                                 >  1:429777/1‑143 (MQ=255)
                                                                         gCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCACAGTTCATCCAGAGCCGCGGCAATCTAGCCGCGTAGAGCGAACGAAACATGCCACGGGTCGCGTAACACACGGCTGACCGTCTTTTTGGCCACGCCTACACg                                                                 >  2:219690/1‑143 (MQ=255)
                                                                               gagagGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGGCACGCCTACACGGTCacc                                                           >  2:555646/1‑141 (MQ=255)
                                                                                  agGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCCATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGCGTTGCGTAAAAAACAGCTGACCGTCAGCTTGGTCACGCCCCCACCGCCAGCcac                                                        >  1:607747/1‑143 (MQ=255)
                                                                                  agGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCCCACCGGCATTTTGGTCACGCCTACACGGTCAGCcac                                                        >  1:407433/1‑143 (MQ=255)
                                                                                       gtCGGGCGCACGATGGGGAATAGAGCCCCGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCCTCATTTTGGTCACGCCTACACGGTCAGCCACATCCt                                                   <  1:569842/142‑1 (MQ=255)
                                                                                        tCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTg                                                  <  1:150517/143‑1 (MQ=255)
                                                                                         cGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGATGCGTAAAAAACGGCTGACCGTCCTTCTGGTCACGCCTACCCCCTCAGCCCCATCCTGa                                                 >  1:295232/1‑143 (MQ=255)
                                                                                          gggCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGaa                                                >  1:484131/1‑143 (MQ=255)
                                                                                              gcACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCGTGCCGCGTAGAGCGACGGAAACCTGCTCCCGGGTGCGTAAAAAACCACTGACCGTCACTGTGGTCACGCCTACACGGTCACCCACATCCTGCAGTAc                                            >  2:300480/1‑143 (MQ=255)
                                                                                                aCGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGCTGCGTAAAAAACGGGAGACGGACCTTTTGGTCACGCCTACACGGTCAGCCCCATACCGAAATACgg                                          >  2:348396/1‑143 (MQ=255)
                                                                                                  gATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGt                                        >  1:754294/1‑143 (MQ=255)
                                                                                                  gATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGt                                        <  1:438394/143‑1 (MQ=255)
                                                                                                         aaTATACCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAATCGGCTCACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGtcttttct                                 <  1:288731/143‑1 (MQ=255)
                                                                                                                cccAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGACTTTTCTTTTTCAt                          >  2:443363/1‑143 (MQ=255)
                                                                                                                 ccAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTCAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGtcttttctttgtagtc                         >  1:58149/1‑137 (MQ=255)
                                                                                                                         tcaAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTCCGGGTCTTTTCTTTTTCATCGTCCTGaa                 >  1:507323/1‑143 (MQ=255)
                                                                                                                           aaGAGCCGCGGCAATCTTGCCGCGTAGAGCGTCAGAAACCGGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAgg               <  1:831143/143‑1 (MQ=255)
                                                                                                                            ctctacgcgGCAAGATTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATTGTCCTGGAGGt              >  1:318188/6‑143 (MQ=255)
                                                                                                                                ccGCGGCTATTTTTCCTCGTATAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACaa          <  2:49381/143‑1 (MQ=255)
                                                                                                                                acgcgGCAAGATTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACaa          >  2:475493/2‑143 (MQ=255)
                                                                                                                                   cgGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCATGGAGGTACAAAag       >  1:537546/1‑143 (MQ=255)
                                                                                                                                     gCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTCTCTTTTTCATCGTCCTGACGCTACAAAagag     >  2:112975/1‑143 (MQ=255)
                                                                                                                                      cAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACAAAagagg    >  1:647490/1‑142 (MQ=255)
                                                                                                                                      cAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACAAAAgaga    <  1:34613/143‑1 (MQ=255)
                                                                                                                                        aTCTGGCCGCGTATAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACAAAAGagata  <  2:813948/143‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                             
GATTCCGCGTAATACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACAAAAGAGATA  >  NC_000913/3578472‑3578750

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