Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 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_000913_3_pae_tpiA1,197,5020CT97.0% 92.5 / NA 33G11R (GGA→AGA) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (14/18);  minor base G (0/1);  total (14/19)
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

ATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGA  >  NC_000913_3_pae_tpiA/1197364‑1197646
                                                                                                                                          |                                                                                                                                                
aTTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGCGCAtatta                                                                                                                                        >  2:219051/1‑149 (MQ=255)
aTTCATCAGAATATCTAAGTTTTCCCCAGCCAAAATCAAGGGCGTGTCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGTTTTTGCTGCATTCTCTGGGCGAATATTAACTCTGTGCAtatta                                                                                                                                        <  2:360429/149‑1 (MQ=255)
  tcatcaGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAACGTTTATACTTTACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATtata                                                                                                                                      <  2:417571/149‑1 (MQ=255)
     tcaGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTCATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGtt                                                                                                                                   >  1:268947/1‑149 (MQ=255)
      caGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTc                                                                                                                                  <  1:573944/149‑1 (MQ=255)
      caGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTc                                                                                                                                  <  1:573857/149‑1 (MQ=255)
         aaTATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGAATCTAATATACTCGAAAAGATATCTAGATGACATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCaaa                                                                                                                               <  1:647423/149‑1 (MQ=255)
               tAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAAGTTTTGCTACATTCTCTGAGCGAATATTAACTCTGTGGATATTATAGTTTAAATCAAGt                                                                                                                         >  2:371955/1‑149 (MQ=255)
                   ttttCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCa                                                                                                                     >  2:448527/1‑149 (MQ=255)
                                     aaGGGCGTGGCCATTTTTCTCAATGCTTCAAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTAtttt                                                                                                   <  1:637283/149‑1 (MQ=255)
                                       gggCGGTGCTATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCGCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTa                                                                                                 <  1:266670/149‑1 (MQ=255)
                                                 aTTTTTCTCAAGGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACGCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAAtt                                                                                       >  1:420418/1‑149 (MQ=255)
                                                       ctcATTGCTTCTAATATACTCGCAAAGACATTTAGAGGGCAACGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTGTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCtttta                                                                                 <  2:198622/149‑1 (MQ=255)
                                                        tcAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttat                                                                                >  1:311201/1‑149 (MQ=255)
                                                              gttttAATATACTCGCAAAGATATCTAGAGGTCAACGTATGATGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATc                                                                          <  2:653187/145‑1 (MQ=255)
                                                                      ataCTCGAAAAAAAATATAAAGGGTATCGTATGAGGTTTTGATGCATTCTATGAGCGAATATTAACTCTCTGCATATTATAGTACAAATCAAGTACCATTTTTTTAATGTAAAAAAAAAAAAAAAATTCTTTTATTTTATCGTCCGGGa                                                                  <  2:580444/149‑1 (MQ=255)
                                                                            gCAAAGATATCTAGAGGGCATCGTATGAGGTTTTTCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTc                                                            >  1:165632/1‑149 (MQ=255)
                                                                                aGATATTTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTCAAAACTTCTTTTATTTTATCGTCCGGGATTTTTCCgac                                                        <  2:116463/149‑1 (MQ=255)
                                                                                aGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgac                                                        <  1:448527/149‑1 (MQ=255)
                                                                                aGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGAATTTTTCCgac                                                        >  2:545893/1‑149 (MQ=255)
                                                                                     tCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgacttgac                                                   >  2:325823/1‑149 (MQ=255)
                                                                                        agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACttt                                                <  1:209360/149‑1 (MQ=255)
                                                                                         gagGGCATCTTATGAGGTGTGGCTGCATTCTCTGGGCGAATATTAACTCTGTGCAGATTATAGTTCGAAGCAAGTACCATTCAGTTAAGGTTATTTTTAGTGAAGAGTACTTTCATTTTATCGTCCGGGATTTTTCGTACTTGACTTTc                                               >  2:540704/1‑149 (MQ=255)
                                                                                                     tGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCtt                                   >  1:359527/1‑149 (MQ=255)
                                                                                                     tGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGAATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTACGGGATTTTTCCGACTTGACTTTCAGTTCTGGGGtt                                   >  1:263903/1‑149 (MQ=255)
                                                                                                       aGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTaa                                 <  2:372935/149‑1 (MQ=255)
                                                                                                          ttttGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttat                                                                                <  2:25289/99‑1 (MQ=255)
                                                                                                          ttttGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttat                                                                                >  1:25289/1‑99 (MQ=255)
                                                                                                                     tcgctGAGCGAATATTAACTCTGTGTATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTTTTTTATTTTATCGTCCGGGATTTTTTCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagcaaag                   <  2:280720/146‑1 (MQ=255)
                                                                                                                      gtctGAGCGAAGATGAACTTGGCGTATATTATAGTTCAAAGCATGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTATTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGa                  <  2:70471/148‑1 (MQ=255)
                                                                                                                       tctGAGCGTATATTAACGCTGTGCATATTATAGTGCAAAGCAAGTACCATTCAAATAATGTTATTTTTAGTGAAAAATCCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGat                 <  1:418631/149‑1 (MQ=255)
                                                                                                                               aaTATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTGTTGAAACGGCCGGGATCCATCGTTCTTTTCTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcac         >  1:130305/1‑149 (MQ=255)
                                                                                                                                      aCTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTga  <  1:219051/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                
ATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGA  >  NC_000913_3_pae_tpiA/1197364‑1197646

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