Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,600,288 C→T S128L (TCG→TTG)  gcvR → transcriptional repressor, regulatory protein accessory to GcvA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,600,2880CT98.4% 212.9 / ‑4.4 62S128L (TCG→TTG) gcvRtranscriptional repressor, regulatory protein accessory to GcvA
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (31/29);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GTCTGGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTCGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGA  >  NC_000913/2600191‑2600381
                                                                                                 |                                                                                             
gtctgGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTtgc                                                                                             >  2:345457/1‑100 (MQ=255)
gtctgGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTtgc                                                                                             <  1:1445264/100‑1 (MQ=255)
  ctgGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTgcgc                                                                                           <  1:722319/100‑1 (MQ=255)
  ctgGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTgcgc                                                                                           >  1:85339/1‑100 (MQ=255)
    gggTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcac                                                                                         >  1:1394488/1‑100 (MQ=255)
     ggTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcacg                                                                                        <  1:729994/100‑1 (MQ=255)
      gTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcacgc                                                                                       >  2:1179661/1‑100 (MQ=255)
      gTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcacgc                                                                                       <  2:361263/100‑1 (MQ=255)
      gTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcacgc                                                                                       >  2:1075675/1‑100 (MQ=255)
       ttaaGGTCGATGTGGCAGCCTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGGCGCTCATCATCTGAACATTGCGGAGCTGGTGTGGcgcacg                                                                                        <  2:524333/95‑1 (MQ=255)
        tCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCaa                                                                                     <  2:1532495/100‑1 (MQ=255)
          aGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAAcc                                                                                   <  1:1494657/100‑1 (MQ=255)
          aGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAAcc                                                                                   >  2:1129887/1‑100 (MQ=255)
            gTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACctg                                                                                 <  2:1098813/100‑1 (MQ=255)
              cGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGTGCGCAAcaccag                                                                               <  1:775456/100‑6 (MQ=255)
              cGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACctgct                                                                               >  1:973275/1‑100 (MQ=255)
              cGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACctgct                                                                               >  1:47262/1‑100 (MQ=255)
               gATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGTGCGCAAcaccagc                                                                              <  1:1572846/100‑7 (MQ=255)
               gATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACctgctg                                                                              >  1:270084/1‑100 (MQ=255)
                aTGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGa                                                                             <  1:196649/100‑1 (MQ=255)
                    ggCAGACTCCCCGCATTTAATTGTGCGCTTCACAGAACTTTTCGAAGCGCATCATATGAACATTGCGGAGCTGGTGTCGCGAAAGGATCCTGCTGaaaa                                                                          <  2:1082229/99‑1 (MQ=255)
                    ggCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAAt                                                                         <  1:810580/100‑1 (MQ=255)
                    ggCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAAt                                                                         >  2:194777/1‑100 (MQ=255)
                    ggCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAAt                                                                         <  1:434322/100‑1 (MQ=255)
                      cAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGa                                                                       >  2:1185399/1‑100 (MQ=255)
                          cTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACgg                                                                   <  1:371429/100‑1 (MQ=255)
                            ccccGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGc                                                                 <  1:1045269/100‑1 (MQ=255)
                            ccccGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGc                                                                 <  1:1435985/100‑1 (MQ=255)
                             cccGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCt                                                                >  2:458284/1‑100 (MQ=255)
                             cccGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCt                                                                >  1:513603/1‑100 (MQ=255)
                             cccGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCt                                                                >  1:864714/1‑100 (MQ=255)
                              ccGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTg                                                               >  1:1060883/1‑100 (MQ=255)
                                  aTTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCa                                                           >  2:572299/1‑100 (MQ=255)
                                      aaTTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTg                                                       <  2:216307/100‑1 (MQ=255)
                                      aaTTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTg                                                       <  1:933511/100‑1 (MQ=255)
                                       aTTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGc                                                      >  1:686739/1‑100 (MQ=255)
                                         tGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCat                                                    <  2:439694/100‑1 (MQ=255)
                                          naaCGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCata                                                   >  2:566905/2‑100 (MQ=255)
                                            aCGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATAtt                                                 >  1:121919/1‑100 (MQ=255)
                                             cGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTc                                                <  2:270084/100‑1 (MQ=255)
                                             atcTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTc                                                <  2:44701/98‑1 (MQ=255)
                                                ttCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGa                                             <  1:955084/100‑1 (MQ=255)
                                                ttCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGa                                             <  1:1336226/100‑1 (MQ=255)
                                                 tCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGAt                                            >  1:44701/1‑100 (MQ=255)
                                                  cacaGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGAt                                            >  2:99430/1‑99 (MQ=255)
                                                  cacaGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATa                                           <  1:1433728/100‑1 (MQ=255)
                                                    caGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAAc                                         >  1:593473/1‑100 (MQ=255)
                                                    caGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAAc                                         >  1:1497000/1‑100 (MQ=255)
                                                               gACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGcc                              >  1:314412/1‑100 (MQ=255)
                                                                aCGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGccc                             <  1:1556129/100‑1 (MQ=255)
                                                                  gcgcATCATATGAACATTGCGGAGCTATTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCg                           <  1:345457/100‑1 (MQ=255)
                                                                        catATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTg                     >  1:1028757/1‑100 (MQ=255)
                                                                          tatGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCg                   <  2:589263/100‑1 (MQ=255)
                                                                          tatGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCg                   >  1:1178274/1‑100 (MQ=255)
                                                                          tatGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCg                   >  1:1015532/1‑100 (MQ=255)
                                                                           atGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCgg                  <  2:1155662/100‑1 (MQ=255)
                                                                           atGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCgg                  >  1:1159681/1‑100 (MQ=255)
                                                                             gAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGAc                >  2:1370735/1‑100 (MQ=255)
                                                                               aCATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACgc              <  2:510799/100‑1 (MQ=255)
                                                                                cATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACgca             <  1:1002761/100‑1 (MQ=255)
                                                                                cATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACgca             >  2:933174/1‑100 (MQ=255)
                                                                                cATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACgca             >  1:1541321/1‑100 (MQ=255)
                                                                                   tGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACgcagca          <  1:194777/100‑1 (MQ=255)
                                                                                           tGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGa  >  1:6086/1‑100 (MQ=255)
                                                                                                 |                                                                                             
GTCTGGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTCGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGA  >  NC_000913/2600191‑2600381

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