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

CGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑T‑CCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGTCATCACCTGTAA  >  NC_000913/2682018‑2682280
                                                                                                                                      |                                                                                                                                  
cGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑T‑cccccc                                                                                                                              <  1:1050034/139‑1 (MQ=255)
cGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑T‑cccccc                                                                                                                              >  1:3745/1‑139 (MQ=255)
  agaAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑T‑cccccccc                                                                                                                            >  2:358492/1‑138 (MQ=255)
                     gAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGcc                                                                                                          <  1:358492/139‑1 (MQ=255)
                        tagtCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCccac                                                                                                       >  2:791158/1‑139 (MQ=255)
                        tagtCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCccac                                                                                                       >  2:65148/1‑139 (MQ=255)
                        tagtCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCccac                                                                                                       >  2:138851/1‑139 (MQ=255)
                          gtCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCccacca                                                                                                     >  1:647669/1‑139 (MQ=255)
                          gtCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCcca                                                                                                        <  2:501835/136‑1 (MQ=255)
                          gtCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCcca                                                                                                        >  1:501835/1‑136 (MQ=255)
                                tttGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCccac                                                                                                       >  1:1094821/1‑131 (MQ=255)
                                tttGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGaa                                                                                               <  1:161223/139‑1 (MQ=255)
                                tttGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTCATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCccac                                                                                                       <  2:1094821/131‑1 (MQ=255)
                                  tGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAAt                                                                                             >  1:212194/1‑139 (MQ=255)
                                     agtaAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑T‑cccccc                                                                                                                              <  2:428713/102‑1 (MQ=255)
                                     agtaAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑T‑cccccc                                                                                                                              >  1:428713/1‑102 (MQ=255)
                                     agtaAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑T‑ccc                                                                                                                                 <  2:716919/99‑1 (MQ=255)
                                     agtaAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑T‑ccc                                                                                                                                 >  1:716919/1‑99 (MQ=255)
                                                gTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑T‑cccccc                                                                                                                              >  1:624516/1‑91 (MQ=255)
                                                gTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑T‑cccccc                                                                                                                              <  2:624516/91‑1 (MQ=255)
                                                 tAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCggcgtggcgt                                                                              <  2:1090777/139‑1 (MQ=255)
                                                      aTGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCg                                                                         >  2:1204453/1‑139 (MQ=255)
                                                           gAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTtata                                                                    <  1:115347/139‑1 (MQ=255)
                                                                    gaCGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑T‑CCCCCCCcttt                                                                                                                         <  1:714252/76‑2 (MQ=255)
                                                                    gaCGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑T‑CCCCCCCcttt                                                                                                                         >  2:714252/1‑75 (MQ=255)
                                                                         cccGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGa                                                      >  1:1054502/1‑139 (MQ=255)
                                                                                  cGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACcaca                                                                                                   <  2:1314977/85‑1 (MQ=255)
                                                                                  cGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACcaca                                                                                                   >  1:1314977/1‑85 (MQ=255)
                                                                                  cGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAgcg                                             >  1:228312/1‑139 (MQ=255)
                                                                                             acCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCa                                  <  2:1224990/139‑1 (MQ=255)
                                                                                                gccgTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCggg                                                                                                                <  2:646808/58‑1 (MQ=255)
                                                                                                gccgTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCggg                                                                                                                >  1:646808/1‑58 (MQ=255)
                                                                                                gccgTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCggcgtg                                                                                  <  2:1652565/88‑1 (MQ=255)
                                                                                                gccgTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCggcgtg                                                                                  >  1:1652565/1‑88 (MQ=255)
                                                                                                gccgTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAggg                               >  2:1151085/1‑139 (MQ=255)
                                                                                                gccgTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAggg                               >  2:1235307/1‑139 (MQ=255)
                                                                                                gccgTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAggg                               >  2:102227/1‑139 (MQ=255)
                                                                                                 ccgTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGt                                                                                                               <  2:675616/58‑1 (MQ=255)
                                                                                                 ccgTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGt                                                                                                               >  1:675616/1‑58 (MQ=255)
                                                                                                 ccgTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGcc                                                                                                          <  2:496601/63‑1 (MQ=255)
                                                                                                 ccgTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGcc                                                                                                          >  1:496601/1‑63 (MQ=255)
                                                                                                  cgTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGgcg                             <  1:791158/139‑1 (MQ=255)
                                                                                                   gTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGgcgc                            <  2:1054502/139‑1 (MQ=255)
                                                                                                         tCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCcca                                                                                                        <  2:275708/57‑1 (MQ=255)
                                                                                                         tCCGGCGGGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCcca                                                                                                        >  1:275708/1‑57 (MQ=255)
                                                                                                             gcggGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAAc                  <  2:212194/139‑1 (MQ=255)
                                                                                                              cggGAAGACGCTCTGATGCCCTCC‑CCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAAcc                 >  2:1493417/1‑139 (MQ=255)
                                                                                                              cggGAAGACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAAcc                 >  2:1010292/1‑139 (MQ=255)
                                                                                                                    gACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGtcat           <  1:65148/139‑1 (MQ=255)
                                                                                                                    gACGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGtcat           <  1:102227/139‑1 (MQ=255)
                                                                                                                      cGCTCTGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGtcatca         <  2:3745/139‑1 (MQ=255)
                                                                                                                         tctGATGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGTCATCACCt      <  2:647669/139‑1 (MQ=255)
                                                                                                                             aTGCCCT‑TCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGTCATCACCTGTaa  <  1:138851/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                  
CGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCT‑T‑CCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGTCATCACCTGTAA  >  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: ‑