Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,495,757 +C coding (470/1344 nt) ydcJ → DUF1338 domain‑containing protein YdcJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,495,7531.C93.3% 44.7 / ‑2.8 15coding (466/1344 nt)ydcJDUF1338 domain‑containing protein YdcJ
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (3/11);  total (4/11)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GGCCCATTGATGATGCTTCTCTGGCGCGTAATCCCTTCCGCGTTTTTACCTCCTTACTCCGCCTTGAGCTTATCGAGAACGAAATTTTGCGCCAGAAAGCGGCGGAGATTCTGCGTCAGCGCGATATCTTCACCCCACGTTGTCGACAACTGTTAGAGGAATATGAGCAGCAGGGCGGTTTTAACGAAACACAGGCACAGGAGTTTGTGCAGGAAGCCCTGGAAACGTTTCGCTGGCA  >  NC_000913/1495622‑1495859
                                                                                                                                    |                                                                                                           
ggCCCATTGATGATGCTTCTCTGGCGCGTAATCCCTTCCGCGTTTTTACCTCCTTACTCCGCCTTGAGCTTATCGAGAACGAAATTTTGCGCCAGAAAGCGGCGGAGATTCTGCGTCAGCGCGATATCTTCAc                                                                                                            <  2:365947/133‑1 (MQ=255)
       tgatgaTGCTTCTCTGGCGCGTAATCCCTTCCGCGTTTTTACCTCCTTACTCCGCCTTGAGCTTATCGAGAACGAAATTTTGCGCCAGAAAGCGGCGGAGATTCTGCGTCAGCGCGATATCTTCAcccccacg                                                                                                    <  1:341056/133‑5 (MQ=255)
                     tGGCGCGTAATCCCTTCCGCGTTTTTACCTCCTTACTCCGCCTTGAGCTTATCGAGAACGAAATTTTGCGCCAGAAAGCGGCGGAGATTCTGCGTCAGCGCGATATCTTCACCCCCACGTTGTCGACAACTGt                                                                                       <  2:405060/133‑1 (MQ=255)
                      ggCGCGTAATCCCTTCCGCGTTTTTACCTCCTTACTCCGCCTTGAGCTTATCGAGAACGAAATTTTGCGCCAGAAAGCGGCGGAGATTCTGCGTCAGCGCGATATCTTCACCCCCACGTTGTCGACAACTGtt                                                                                      <  1:39968/133‑1 (MQ=255)
                        cgcgTAATCCCTTCCGCGTTTTTACCTCCTTACTCCGCCTTGAGCTTATCGATAACGAAATTTTGCGCCAGAAAGCGGCGGAGATTCTGCGTCAGCGCGATATCTTCACCCCCACGTTGTCGACAACTGTTag                                                                                    <  1:118712/133‑1 (MQ=255)
                                           ttttACCTCCTTACTCCGCCTTGAGCTTATCGAGAACGAAATTTTGCGCCAGAAAGCGGCGGAGATTCTGCGTCAGCGCGATATCTTCACCCCCACGTTGTCGACAACTGTTAGAGGAATATGAGCAGCAggg                                                                 <  2:256794/133‑1 (MQ=255)
                                           ttttACCTCCTTACTCCGCCTTGAGCTTATCGAGAACGAAATTTTGCGCCAGAAAGCGGCGGAGATTCTGCGTCAGCGCGATATCTTCACCTCCACGTTGCCGATAACTGTTAGAGGTATATTAGCAGCAggg                                                                 >  2:203882/1‑133 (MQ=255)
                                            tttACCTCCTTACTCCGCCTTGAGCTTATCGAGAACGAAATTTTGCGCCAGAAAGCGGCGGAGATTCTGCGTCAGCGCGATATCTTCACCCCCACGTTGTCGACAACTGTTAGAGGAATATGAGCAGCAGGGc                                                                <  2:273871/133‑1 (MQ=255)
                                                       aCTCCGCCTTGAGCTTATCGAGAACGAAATTTTGCGCCAGAAAGCGGCGGAGATTCTGCGTCAGCGCGATATCTTCACCCCCACGTTGTCGACAACTGTTAGAGGAATATGAGCAGCAGGGCGGTTTTAACGa                                                     >  1:411428/1‑133 (MQ=255)
                                                       aCTCCGCCTTGAGCTTATCGAGAACGAAATTTTGCGCCAGAAAGCGGCGGAGATTCTGCGTCAGCGCGATATCTTCACCCCCACGTTGTCGACAACTGTTAGAGGAATATGAGCAGCAGGGCGGTTTTAACGa                                                     <  2:432378/133‑1 (MQ=255)
                                                          ccGCCTTGAGCTTATCGAGAACGAAATTTTGCGCCAGAAAGCGGCGGAGATTCTGCGTCAGCGCGATATCTTCACCCCCACGTTGTCGACAACTGTTAGAGGAATATGAGCAGCAGGGCGGTTTTAACGAAac                                                  >  1:319067/1‑133 (MQ=255)
                                                            gCCTTGAGCTTATCGTGAACGAAATTTTACGCCAGAAAGCGTCGGAGATTCTGCGTCAGCGCGATATCTTCACCCCCACGTTGTCGACAACTGTTAGAGGAATATGAGCAGCAGGGCGGTTTTAACGAAacac                                                <  2:173968/133‑1 (MQ=255)
                                                                tGAGCTTATCGAGAACGAAATTTTGCGCCAGAAAGCGGCGGAGATTCTGCGTCAGCGCGATATCTTCACCCCCACGTTGTCGACAACTGTTAGAGGAATATGAGCAGCAGGGCGGTTTTAACGAAAcacaggc                                            <  2:59760/133‑1 (MQ=255)
                                                                   gCTTATCGAGAACGAAATTTTGCGCCAGAAAGCGGCGGAGATTCTGCGTCAGCGCGATATCTTCACCCCCACGTTGTCGACAACTGTTAGAGGAATATGAGCAGCAGGGCGGTTTTAACGAAAcacaggcaca                                         <  1:431547/133‑1 (MQ=255)
                                                                            gaACGAAATTTTGCGCCAGAAAGCGGCGGAGATTCTGCGTCAGCGCGATATCTTCACCCCCACGTTGTCGACAACTGTTAGAGGAATATGAGCAGCAGGGCGGTTTTAACGAAACACAGGCACAGGAGTTtgt                                >  2:3568/1‑133 (MQ=255)
                                                                                               aaaGCGGCGGTGATTCTGCGTCAGGGCGGTATGTTCACCCCCACGTTGTCGACAACTGTTAGAGGAATATGAGCAGCAGGGCGGTTTTAACGAAACACAGGCACAGGAGTTTGTGCAGGAAGCCCTGGAAACg             <  2:200332/133‑1 (MQ=255)
                                                                                                          gaTTCTGCGTCAGCGCGATATCTTCACCCCCACGTTGTCGACAACTGTTAGAGGAATATGAGCAGCAGGGCGGTTTTAACGAAACACAGGCACAGGAGTTTGTGCAGGAAGCCCTGGAAACGTTTCGCTGGCa  <  2:317090/133‑1 (MQ=255)
                                                                                                                                    |                                                                                                           
GGCCCATTGATGATGCTTCTCTGGCGCGTAATCCCTTCCGCGTTTTTACCTCCTTACTCCGCCTTGAGCTTATCGAGAACGAAATTTTGCGCCAGAAAGCGGCGGAGATTCTGCGTCAGCGCGATATCTTCACCCCACGTTGTCGACAACTGTTAGAGGAATATGAGCAGCAGGGCGGTTTTAACGAAACACAGGCACAGGAGTTTGTGCAGGAAGCCCTGGAAACGTTTCGCTGGCA  >  NC_000913/1495622‑1495859

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