Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,416,279 A→G S370P (TCC→CCC)  recE ← Rac prophage, exonuclease VIII, 5' to 3' specific dsDNA exonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,416,2790AG91.3% 70.1 / 1.1 23coding (1108/2601 nt)recERac prophage, exonuclease VIII, 5' to 3' specific dsDNA exonuclease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (11/10);  minor base . (1/1);  total (12/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.03e-01

GTTGCGTCCGGTTCCATTGTTTCAGCCTCTCCCTGTTCAACTGCCGTTGTTCCAGATGGTTGCGGTTTTTCTTCATCATCCTGTTTTCCTTCTTCTGTTACTCGCTGCGGCATCGGGGCAGAGGAGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGACCTCGATCCCAATTGGTG  >  NC_000913/1416155‑1416395
                                                                                                                            |                                                                                                                    
gTTGCGTCCGGTTCCATTGTTTCAGCCTCTCCCTGTTCAACTGCCGTTGTTCCAGATGGTTGCGGTTTTTCTTCATCATCCTGTTTTCCTTCTTCTGTTACTCGCTGCGGCATCGGGGCAGAGGGGCGACCGCAGGCaa                                                                                                        >  1:238934/1‑139 (MQ=255)
    cGTCCGGTTCCATTGTTTCAGCCTCTCCCTGTTCAACTGCCGTTGTTCCAGATGGTTGCGGTTTTTCTTCATCATCCTGTTTTCCTTCTTCTGTTACTCGCTGCGGCATCGGGGCAGAGGGGCGACCGCAGGCAATATc                                                                                                    >  1:80449/1‑139 (MQ=255)
      tCCGGTTCCATTGTTTCAGCCTCTCCCTGTTCAACTGCCGTTGTTCCAGATGGTTGCGGTTTTTCTTCATCATCCTGTTTTCCTTCTTCTGTTACTCGCTGCGGCATCGGGGCAGAGGGGCGACCGCAGGCAATATCCa                                                                                                  >  2:457855/1‑139 (MQ=255)
                     tCAGCCTCTCCCTGTTCAACTGCCGTTGTTCCAGATGGTTGCGGTTTTTCTTCATCATCCTGTTTTCCTTCTTCTGTTACTCGCTGCGGCATCGGGGCAGAGGGGCGACCGCAGGCAATATCCACGATTTCCGGAt                                                                                      >  2:50560/1‑136 (MQ=255)
                     tCAGCCTCTCCCTGTTCAACTGCCGTTGTTCCAGATGGTTGCGGTTTTTCTTCATCATCCTGTTTTCCTTCTTCTGTTACTCGCTGCGGCATCGGGGCAGAGGGGCGACCGCAGGCAATATCCACGATTTCCGGAt                                                                                      <  1:50560/136‑1 (MQ=255)
                           tctcCCTGTTCAACTGCCGTTGTTCCAGATGGTTGCGGTTTTTCTTCATCATCCTGTTTTCCTTCTTCTGTTACTCGCTGCGGCATCGGGGCAGAGGGGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTgg                                                                             <  1:44656/139‑1 (MQ=255)
                                 tGTTCAACTGCCGTTGTTCCAGATGGTTGCGGTTTTTCTTCATCATCCTGTTTTCCTTCTTCTGTTACTCGCTGCGGCATCGGGGCAGAGGGGCGACCGCAGGCa                                                                                                         >  1:211574/1‑105 (MQ=255)
                                 tGTTCAACTGCCGTTGTTCCAGATGGTTGCGGTTTTTCTTCATCATCCTGTTTTCCTTCTTCTGTTACTCGCTGCGGCATCGGGGCAGAGGGGCGACCGCAGGCa                                                                                                         <  2:211574/105‑1 (MQ=255)
                                   ttCAACTGCCGTTGTTCCAGATGGTTGCGGTTTTTCTTCATCATCCTGTTTTCCTTCTTCTGTTACTCGCTGCGGCATCGGGGCAGAGGGGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCg                                                                     >  2:95665/1‑139 (MQ=255)
                                             gttgttCCAGATGGTTGCGGTTTTTCTTCATCATCCTGTTTTCCTTCTTCTGTTACTCGCTGCGGCATCGGGGCAGAGGGGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTtcagtca                                                           <  1:457855/139‑1 (MQ=255)
                                                gttCCAGATGGTTGCGGTTTTTCTTCATCATCCTGTTTTCCTTCTTCTGTTACTCGCTGCGGCATCGGGGCAGAGGGGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTa                                                        <  1:472936/139‑1 (MQ=255)
                                                  tCCAGATGGTTGCGGTTTTTCTTCATCATCCTGTTTTCCTTCTTCTGTTACTCGCTGCGGCATCGGGGCAGAGGGGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACt                                                      >  1:648291/1‑139 (MQ=255)
                                                                        tcatcaTCCTGTTTTCCTTCTTCTGTTACTCGCTGCGGCATCGGGGCAGAGGGGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGt                                >  1:645955/1‑139 (MQ=255)
                                                                           tcatcCTGTTTTCCTTCTTCTGTTACTCGCTGCGGCATCGGGGCAGAGGGGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTgcg                             >  1:14936/1‑139 (MQ=255)
                                                                                   ttttCCTTCTTCTGTTACTCGCTGCGGCATCGGGGCAGAGGGGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATg                     <  2:645955/139‑1 (MQ=255)
                                                                                         ttcttcTGTTACTCGCTGCGGCATCGGGGCAGAGGGGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGACCTCg               <  2:238934/139‑1 (MQ=255)
                                                                                           cttctGTTACTCGCTGCGGCATCGGGGCAGAGGGGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGACCTCGAt             <  1:95665/139‑1 (MQ=255)
                                                                                                 ttACTCGCTGCGGCATCGGGGCAGGGGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTa                                                        <  1:117239/88‑1 (MQ=255)
                                                                                                 ttACTCGCTGCGGCATCGGGGCAGGGGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTa                                                        >  2:117239/1‑88 (MQ=255)
                                                                                                     tCGCTGCGGCATCGGGGCAGAGGGGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGACCTCGATCCCAATTGGt   >  1:535801/1‑139 (MQ=255)
                                                                                                      cGCTGCGGCATCGGGGCAGAGGGGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGACCTCGATCCCAATTGGTg  <  2:535801/139‑1 (MQ=255)
                                                                                                         tGCGGCATCGGGGCAGAGGGGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAg                                              >  2:282809/1‑92 (MQ=255)
                                                                                                         tGCGGCATCGGGGCAGAGGGGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAg                                              <  1:282809/92‑1 (MQ=255)
                                                                                                                            |                                                                                                                    
GTTGCGTCCGGTTCCATTGTTTCAGCCTCTCCCTGTTCAACTGCCGTTGTTCCAGATGGTTGCGGTTTTTCTTCATCATCCTGTTTTCCTTCTTCTGTTACTCGCTGCGGCATCGGGGCAGAGGAGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGACCTCGATCCCAATTGGTG  >  NC_000913/1416155‑1416395

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