Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,197,502 C→T G11R (GGA→AGA)  ymfD ← e14 prophage; putative SAM‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,197,5020CT92.9% 71.2 / NA 28G11R (GGA→AGA) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (9/17);  minor base A (1/1);  total (10/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCT  >  NC_000913/1197437‑1197563
                                                                 |                                                             
cTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                             >  2:750926/1‑68 (MQ=255)
cTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                             <  1:722063/68‑1 (MQ=255)
cTCGCAAAGATATCTAGAGGGCATCGTACGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                             <  1:540253/68‑1 (MQ=255)
  cGCAAAGATATCTAGAGTGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                           <  2:688014/68‑1 (MQ=255)
  cGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                           <  1:590716/68‑1 (MQ=255)
    cAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCat                                                         <  2:406305/68‑1 (MQ=255)
    cAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCat                                                         <  1:210571/68‑1 (MQ=255)
    cAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCat                                                         <  2:247584/68‑1 (MQ=255)
    cAAAGATATCTAGAGGGCATCGTATGAGGTTTCGCTGCATTCTCTGAGCGAATATTAACTCTGTGCat                                                         <  2:467154/68‑1 (MQ=255)
       aGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtatt                                                      <  1:289239/68‑1 (MQ=255)
            tCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGt                                                 <  2:546029/68‑1 (MQ=255)
               agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCa                                              >  1:289472/1‑68 (MQ=255)
               agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCa                                              >  2:645804/1‑68 (MQ=255)
                  gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                           <  2:535213/68‑1 (MQ=255)
                  gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                           <  2:445070/68‑1 (MQ=255)
                   tgCATCTTAGTAGGTTTTTCTGCATTCTATGATCGAATATTAACTCTGTGCTTTTTATAGTTCAAAGc                                          <  1:168193/67‑1 (MQ=255)
                   ggCAACGAATGAGGTTTTGCTGCATTCTATTCGCGAATAAAAACTCAGACCATATAATAGTTCAAAGc                                          >  2:168193/1‑68 (MQ=255)
                           aTGAGGTTTTGCATCATTCTCTGATCGAATATTAACTCAGTTCATATTATAGTTCAAAGCAAGTACCa                                  <  2:213986/68‑1 (MQ=255)
                               ggTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCa                              <  1:470228/68‑1 (MQ=255)
                               ggTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCa                              >  1:727192/1‑68 (MQ=255)
                                gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtatt                                                      <  2:136392/43‑1 (MQ=255)
                                gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtatt                                                      >  1:136392/1‑43 (MQ=255)
                                gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCaa                             >  1:308616/1‑68 (MQ=255)
                                      ctgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATg                       >  2:279183/1‑68 (MQ=255)
                                            tctcTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTAttt                 >  1:235407/1‑68 (MQ=255)
                                                   gCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGa          <  1:280930/68‑1 (MQ=255)
                                                    cGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGaa         >  2:306759/1‑68 (MQ=255)
                                                           tAACTCTGTTCATATTATAGTTCAAAGCAAGTACCATTCTATTAATGTTATTTTTAGTGAAAAATTCt  <  1:102671/68‑1 (MQ=255)
                                                                 |                                                             
CTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCT  >  NC_000913/1197437‑1197563

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