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% 219.0 / ‑4.3 64S128L (TCG→TTG) gcvRtranscriptional repressor, regulatory protein accessory to GcvA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (32/31);  total (33/31)
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.94e-01

TGGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTCGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGAGC  >  NC_000913/2600194‑2600383
                                                                                              |                                                                                               
tgGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgca                                                                                            >  1:1092108/1‑100 (MQ=255)
 gggTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcac                                                                                           >  1:751484/1‑100 (MQ=255)
  ggTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcacg                                                                                          >  2:372098/1‑100 (MQ=255)
    ttCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcacgca                                                                                        >  2:1182522/1‑100 (MQ=255)
      cAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAAc                                                                                      >  1:1564863/1‑100 (MQ=255)
       aGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAAcc                                                                                     >  2:183454/1‑100 (MQ=255)
         gTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACctg                                                                                   <  1:1455065/100‑1 (MQ=255)
            gATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGTGCGCAAcaccagc                                                                                <  1:748510/100‑7 (MQ=255)
            gATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACctgctg                                                                                <  2:279838/100‑1 (MQ=255)
            gATGTGGCAGACTCCCCGCATTTAATTAAACGCTTCACAGCACTTTTCAACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACctgct                                                                                 >  1:299897/1‑99 (MQ=255)
               tctgCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGaaa                                                                             <  1:1251237/97‑1 (MQ=255)
               gtgGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGaaa                                                                             <  2:456106/100‑1 (MQ=255)
                tgGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGaaaa                                                                            >  1:810300/1‑100 (MQ=255)
                  gCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGaaaaga                                                                          >  2:1251229/1‑98 (MQ=255)
                   cAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGa                                                                         >  1:1046979/1‑100 (MQ=255)
                      aCTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACg                                                                      <  1:725826/100‑1 (MQ=255)
                          cccGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCt                                                                  >  1:92398/1‑100 (MQ=255)
                            cGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCCACCTGCTGAAAATGAACGGGCTg                                                                 >  2:780990/1‑99 (MQ=255)
                              cATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTgcgc                                                              <  1:793480/100‑1 (MQ=255)
                                tttAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAg                                                            <  2:347192/100‑1 (MQ=255)
                                  tAATTGAACGGGTCACAGCACTTTTGGACGCGCATCATATGAAAATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAACTGAACGGGc                                                                   <  2:1242552/91‑1 (MQ=255)
                                     ttGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCa                                                       >  2:1159765/1‑100 (MQ=255)
                                       gAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCata                                                     >  2:1086467/1‑100 (MQ=255)
                                       gAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCata                                                     >  2:1086732/1‑100 (MQ=255)
                                       gAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCata                                                     >  1:438582/1‑100 (MQ=255)
                                         aCGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATAtt                                                   >  2:1110837/1‑100 (MQ=255)
                                          cGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTc                                                  <  2:1054917/100‑1 (MQ=255)
                                           gCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCGGTTGCATATTCa                                                 <  1:1032242/100‑1 (MQ=255)
                                           gCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCa                                                 <  1:1026577/100‑1 (MQ=255)
                                           gCTTCACAGCACTTTTCGACGCGCATAATATGAAAATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAACATGAACGGGCTGCGCAGTTGCATATTCa                                                 <  1:1254866/100‑1 (MQ=255)
                                            cTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAg                                                >  1:904317/1‑100 (MQ=255)
                                            cTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAg                                                <  1:1182529/100‑1 (MQ=255)
                                             ttCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGa                                               >  1:63820/1‑100 (MQ=255)
                                                 caGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAAc                                           <  2:213557/100‑1 (MQ=255)
                                                 caGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAAc                                           >  2:436117/1‑100 (MQ=255)
                                                 caGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAAc                                           <  2:604865/100‑1 (MQ=255)
                                                   gCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCg                                         >  1:291665/1‑100 (MQ=255)
                                                      cTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGccc                                      >  1:1552967/1‑100 (MQ=255)
                                                        tttCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCcac                                    <  2:1083043/100‑1 (MQ=255)
                                                        tttCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCcac                                    <  1:753573/100‑1 (MQ=255)
                                                         ttCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCcaca                                   <  2:1072135/100‑1 (MQ=255)
                                                          tCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACAGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAg                                  >  1:1101541/1‑100 (MQ=255)
                                                               gcgcATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGccccn                             <  1:887743/100‑2 (MQ=255)
                                                               gcgcATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCg                             <  1:372099/100‑1 (MQ=255)
                                                                   atcatATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATc                         >  1:580064/1‑100 (MQ=255)
                                                                      atatGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGc                      >  2:583710/1‑100 (MQ=255)
                                                                        atGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCgg                    <  2:1137422/100‑1 (MQ=255)
                                                                        atGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCgg                    >  2:1193777/1‑100 (MQ=255)
                                                                        atGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCgg                    <  1:783838/100‑1 (MQ=255)
                                                                         tGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCCTCTGCGGa                   >  1:69754/1‑100 (MQ=255)
                                                                         tGAACATTGCGGAGCTGGTGTCGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGa                   >  2:867354/1‑100 (MQ=255)
                                                                          gAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGAc                  >  2:536153/1‑100 (MQ=255)
                                                                          gAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGAc                  <  2:810295/100‑1 (MQ=255)
                                                                            aCATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACgc                >  2:75428/1‑100 (MQ=255)
                                                                               ttGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACgcagc             >  1:22566/1‑100 (MQ=255)
                                                                                tGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACgcagca            <  1:1086739/100‑1 (MQ=255)
                                                                                tGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACgcagca            <  1:1086474/100‑1 (MQ=255)
                                                                                   ggAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAat         <  1:1509849/100‑1 (MQ=255)
                                                                                      gCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATAtt      >  1:1207840/1‑100 (MQ=255)
                                                                                      gCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATAtt      >  1:1397656/1‑100 (MQ=255)
                                                                                       cTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTg     <  2:1564853/100‑1 (MQ=255)
                                                                                          gtgtTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGAGc  <  2:438580/100‑1 (MQ=255)
                                                                                          gtgtTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGAGc  <  2:480892/100‑1 (MQ=255)
                                                                                          gtgtTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGAGc  <  2:69754/100‑1 (MQ=255)
                                                                                              |                                                                                               
TGGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTCGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGAGC  >  NC_000913/2600194‑2600383

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: