Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,682,157 (C)7→8 coding (589/3282 nt) yphG ← DUF4380 domain‑containing TPR repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,682,1501.C93.2% 144.8 / 3.1 44coding (596/3282 nt)yphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (20/21);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

CGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGTCATCACCTGTAA  >  NC_000913/2682018‑2682280
                                                                                                                                      |                                                                                                                                  
cGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTcccccc                                                                                                                              <  1:1050034/139‑1 (MQ=255)
cGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTcccccc                                                                                                                              >  1:3745/1‑139 (MQ=255)
  agaAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTcccccccc                                                                                                                            >  2:358492/1‑138 (MQ=255)
                     gAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGcc                                                                                                          <  1:358492/139‑1 (MQ=255)
                        tagtCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCccac                                                                                                       >  2:791158/1‑139 (MQ=255)
                        tagtCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCccac                                                                                                       >  2:65148/1‑139 (MQ=255)
                        tagtCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCccac                                                                                                       >  2:138851/1‑139 (MQ=255)
                          gtCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCccacca                                                                                                     >  1:647669/1‑139 (MQ=255)
                          gtCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCcca                                                                                                        <  2:501835/136‑1 (MQ=255)
                          gtCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCcca                                                                                                        >  1:501835/1‑136 (MQ=255)
                                tttGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCccac                                                                                                       >  1:1094821/1‑131 (MQ=255)
                                tttGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGaa                                                                                               <  1:161223/139‑1 (MQ=255)
                                tttGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTCATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCccac                                                                                                       <  2:1094821/131‑1 (MQ=255)
                                  tGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAAt                                                                                             >  1:212194/1‑139 (MQ=255)
                                     agtaAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTcccccc                                                                                                                              <  2:428713/102‑1 (MQ=255)
                                     agtaAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTcccccc                                                                                                                              >  1:428713/1‑102 (MQ=255)
                                     agtaAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTccc                                                                                                                                 <  2:716919/99‑1 (MQ=255)
                                     agtaAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTccc                                                                                                                                 >  1:716919/1‑99 (MQ=255)
                                                gTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTcccccc                                                                                                                              >  1:624516/1‑91 (MQ=255)
                                                gTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTcccccc                                                                                                                              <  2:624516/91‑1 (MQ=255)
                                                 tAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCggcgtggcgt                                                                              <  2:1090777/139‑1 (MQ=255)
                                                      aTGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCg                                                                         >  2:1204453/1‑139 (MQ=255)
                                                           gAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTtata                                                                    <  1:115347/139‑1 (MQ=255)
                                                                    gaCGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCcttt                                                                                                                         <  1:714252/76‑2 (MQ=255)
                                                                    gaCGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCcttt                                                                                                                         >  2:714252/1‑75 (MQ=255)
                                                                         cccGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGa                                                      >  1:1054502/1‑139 (MQ=255)
                                                                                  cGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACcaca                                                                                                   <  2:1314977/85‑1 (MQ=255)
                                                                                  cGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACcaca                                                                                                   >  1:1314977/1‑85 (MQ=255)
                                                                                  cGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAgcg                                             >  1:228312/1‑139 (MQ=255)
                                                                                             acCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCa                                  <  2:1224990/139‑1 (MQ=255)
                                                                                                gccgTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCggg                                                                                                                <  2:646808/58‑1 (MQ=255)
                                                                                                gccgTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCggg                                                                                                                >  1:646808/1‑58 (MQ=255)
                                                                                                gccgTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCggcgtg                                                                                  <  2:1652565/88‑1 (MQ=255)
                                                                                                gccgTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCggcgtg                                                                                  >  1:1652565/1‑88 (MQ=255)
                                                                                                gccgTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAggg                               >  2:1151085/1‑139 (MQ=255)
                                                                                                gccgTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAggg                               >  2:1235307/1‑139 (MQ=255)
                                                                                                gccgTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAggg                               >  2:102227/1‑139 (MQ=255)
                                                                                                 ccgTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGt                                                                                                               <  2:675616/58‑1 (MQ=255)
                                                                                                 ccgTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGt                                                                                                               >  1:675616/1‑58 (MQ=255)
                                                                                                 ccgTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGcc                                                                                                          <  2:496601/63‑1 (MQ=255)
                                                                                                 ccgTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGcc                                                                                                          >  1:496601/1‑63 (MQ=255)
                                                                                                  cgTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGgcg                             <  1:791158/139‑1 (MQ=255)
                                                                                                   gTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGgcgc                            <  2:1054502/139‑1 (MQ=255)
                                                                                                         tCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCcca                                                                                                        <  2:275708/57‑1 (MQ=255)
                                                                                                         tCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCcca                                                                                                        >  1:275708/1‑57 (MQ=255)
                                                                                                             gcggGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAAc                  <  2:212194/139‑1 (MQ=255)
                                                                                                              cggGAAGACGCTCTGATGCCCTCCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAAcc                 >  2:1493417/1‑139 (MQ=255)
                                                                                                              cggGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAAcc                 >  2:1010292/1‑139 (MQ=255)
                                                                                                                    gACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGtcat           <  1:65148/139‑1 (MQ=255)
                                                                                                                    gACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGtcat           <  1:102227/139‑1 (MQ=255)
                                                                                                                      cGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGtcatca         <  2:3745/139‑1 (MQ=255)
                                                                                                                         tctGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGTCATCACCt      <  2:647669/139‑1 (MQ=255)
                                                                                                                             aTGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGTCATCACCTGTaa  <  1:138851/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                  
CGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGTCATCACCTGTAA  >  NC_000913/2682018‑2682280

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: