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.C97.1% 120.4 / ‑2.5 35coding (596/3282 nt)yphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (20/14);  total (21/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

CAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTT‑‑CCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGTCATCACCTGTAAACCATGCATCGGC  >  NC_000913/2682001‑2682293
                                                                                                                                                        |                                                                                                                                                
cAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCtt                                                                                                                                                   >  1:673147/1‑150 (MQ=255)
        tttttATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTT‑‑cccccccct                                                                                                                                        >  1:40286/1‑149 (MQ=255)
         ttttAAAGGGAAAAATGCCCACCCCAGGGGAGAAGCCCCTTTTGTAGTAGGTGCGGGTGGGGATGGGGAAAGGGGGGCGGGCCCTTTGGCGGGGCCAAACCCCCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTT‑‑ccccccc                                                                                                                                          <  2:501541/148‑1 (MQ=255)
           ttATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACACTCTGATGCCCTT‑‑CCCCCCCctttc                                                                                                                                     >  1:237320/1‑149 (MQ=255)
                gCGAGAAATGTCCACCCCAGCGGAGTAGTCCATTTTGTAGTAAGTGCGGGTAGCGATGGGGAAAGGGGAGACGGCCTTTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGc                                                                                                                                 <  1:53052/151‑1 (MQ=255)
                        tGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTgg                                                                                                                         >  1:132930/1‑151 (MQ=255)
                        tGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTgg                                                                                                                         >  1:216803/1‑151 (MQ=255)
                           ccACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGccc                                                                                                                      >  2:209635/1‑151 (MQ=255)
                                ccAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCcaccac                                                                                                                 >  2:411746/1‑151 (MQ=255)
                                   gCGGAGTAGTCCATTTTTTAGTATGTGCGGGAGGGGGGGGGAAAAGGGGAGACGCCCTTTTTGCGGTGATAAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGtt                                                                                                                           <  2:7165/138‑1 (MQ=255)
                                          agtCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGAc                                                                                                       >  2:448299/1‑151 (MQ=255)
                                              cACTTTGTAGTAATGCCGGGTAGCGAGGGGAAAAGCGGAAACGGCCCGTTTGCCGTGATCAAACTCCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGAcg                                                                                                      <  1:404758/148‑1 (MQ=255)
                                                        taAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTc                                                                                         >  1:490734/1‑151 (MQ=255)
                                                                 ggaggggaggggaAAAGGGGGGGCGCCCTGTTTGCGTTGACAAACCACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTtat                                                                                  <  2:178858/141‑1 (MQ=255)
                                                                   tAGCGATGGGGAAAGGGGAGACGGCCCGTTTGCCGTGATAAACCCCCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAg                                                                                <  1:183885/150‑1 (MQ=255)
                                                                      cGAGGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGAcg                                                                             >  2:663405/1‑151 (MQ=255)
                                                                           gggAAGGCGGAGACGGCCCTTTTGCGGTGATAAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTg                                                                       <  1:50339/151‑1 (MQ=255)
                                                                           gggAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTg                                                                       >  1:304744/1‑151 (MQ=255)
                                                                             gAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTgg                                                                      >  2:452305/1‑150 (MQ=255)
                                                                                         gcccttttggctgTGATAAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAgc                                                           <  1:642799/137‑1 (MQ=255)
                                                                                                   ccGTAATCAACCACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTcc                                               <  1:11/151‑2 (MQ=255)
                                                                                                     gTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAgg                                             >  1:642684/1‑151 (MQ=255)
                                                                                                     gTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAgg                                             >  1:152795/1‑151 (MQ=255)
                                                                                                      ggctacaacggccgccgTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCa                                               <  1:123603/139‑1 (MQ=255)
                                                                                                        aacAAACCCCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGt                                                 <  2:296757/142‑1 (MQ=255)
                                                                                                         tCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGgcg                                          >  2:97863/1‑150 (MQ=255)
                                                                                                            aaCCCCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGgc                                           <  2:493818/146‑1 (MQ=255)
                                                                                                               cacCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAggg                                     >  2:620659/1‑150 (MQ=255)
                                                                                                                       ttACATCCGGCGGGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTg                            >  1:55597/1‑151 (MQ=255)
                                                                                                                               ggcggGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGTCATCAcc                    >  1:43541/1‑151 (MQ=255)
                                                                                                                               ggcggGAAGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGTCATCAcc                    <  1:297154/151‑1 (MQ=255)
                                                                                                                                     aaGACGCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGTCATCACCTGTaa               >  2:157805/1‑150 (MQ=255)
                                                                                                                                          gCTCTGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGTCATCACCTGTAAACCATg         >  2:125779/1‑151 (MQ=255)
                                                                                                                                            tctGATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGTCATCACCTGTAAACCATGc        >  2:243190/1‑150 (MQ=255)
                                                                                                                                               gATGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGTCATCACCTGTAAACCATGCa       <  1:197052/148‑1 (MQ=255)
                                                                                                                                                aTGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGTCATCACCTGTAAACCATGCATCgg   >  1:92948/1‑151 (MQ=255)
                                                                                                                                                aTGCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGTCATCACCTGTAAACCATGCATCgg   >  1:455906/1‑151 (MQ=255)
                                                                                                                                                  gCCCTTCCCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGTCATCACCTGTAAACCATGCATCGGc  <  1:349142/150‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                
CAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGTTACATCCGGCGGGAAGACGCTCTGATGCCCTT‑‑CCCCCCCTTTCACTGCCGGGTTGGCCCACCACAAGAAATGACGCGGCGTGGCGTTCCCGTTATAGACGCGGCTGGCGATTTCCAGCGCCGCCCGGTCAGGGCGCAGGGTGAAACCTGTCATCACCTGTAAACCATGCATCGGC  >  NC_000913/2682001‑2682293

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