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.7% 134.4 / NA 43G11R (GGA→AGA) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (21/21);  total (21/22)
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

ATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACT  >  NC_000913_3_pae_tpiA/1197360‑1197644
                                                                                                                                              |                                                                                                                                              
aTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAATATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCat                                                                                                                                          >  1:983779/1‑149 (MQ=255)
  tAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCatat                                                                                                                                        <  1:732664/149‑1 (MQ=255)
     ttCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtattat                                                                                                                                     >  2:49284/1‑149 (MQ=255)
       cattcgcatatCTAAGTTTTCCGCAGTCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAg                                                                                                                                   <  2:526893/142‑1 (MQ=255)
            gAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCaa                                                                                                                              >  2:920808/1‑149 (MQ=255)
             aaTATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCaaa                                                                                                                             <  1:920808/149‑1 (MQ=255)
                  cTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAg                                                                                                                        >  2:168834/1‑149 (MQ=255)
                   tAAGTTTTCCGCAGCCAAAATCAAGAGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTATCCATATTATAGTTCAAACAAAGt                                                                                                                       >  2:931290/1‑149 (MQ=255)
                       ttttCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTTTTATTGTTTCTAATATAAACGCAAAGATATCAAAAACGCATCGTATGAGGTTTTCCTGCATTCTCTTATCGAATATTAACTCTGTGCATCATATATATCAAACCGAGTCCCa                                                                                                                   >  2:483692/1‑149 (MQ=255)
                                 ccAAAATCAAGGGCTTGGCCTATTTTCTCAATTTTTCTAATATTCTCGCAAAGATATCTATAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCACGTACCATTCAATTAAt                                                                                                         <  2:594400/149‑1 (MQ=255)
                                 ccAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAAATCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAt                                                                                                         <  1:542345/149‑1 (MQ=255)
                                           gggCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTa                                                                                               <  2:1013523/149‑1 (MQ=255)
                                                        tttCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCtt                                                                                  <  1:168834/149‑1 (MQ=255)
                                                              aaTGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTATTGAAAAATTCttttattt                                                                            >  2:277602/1‑149 (MQ=255)
                                                                tGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                                                                                                            <  2:681512/83‑1 (MQ=255)
                                                                tGCTTCTAATACACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                                                                                                            >  1:681512/1‑83 (MQ=255)
                                                                  cTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATc                                                                        >  2:859475/1‑149 (MQ=255)
                                                                    tCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGt                                                                      >  2:63566/1‑149 (MQ=255)
                                                                        atatGCTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCgg                                                                  <  2:734054/149‑1 (MQ=255)
                                                                                gCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTACATCTTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTc                                                          >  2:951821/1‑149 (MQ=255)
                                                                                            agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTAAAAGCAAGTACTATTCAATTAATGTTATTTTTAGTGACAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACttt                                              <  2:383433/149‑1 (MQ=255)
                                                                                             gagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTc                                             >  2:685083/1‑149 (MQ=255)
                                                                                                   aTCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGAATTTTTCCGACTTGACTTTCAGTTCt                                       >  1:717506/1‑149 (MQ=255)
                                                                                                    tCGTATGGGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTg                                      <  1:255509/149‑1 (MQ=255)
                                                                                                    tCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTg                                      <  2:765558/149‑1 (MQ=255)
                                                                                                         tGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCtt                                 <  1:185042/149‑1 (MQ=255)
                                                                                                               tttGCTTCATTCTCTGAGCGAATAGTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGt                           <  1:746024/149‑1 (MQ=255)
                                                                                                                   aggagttCTCTAATCGAATATTAACTCTGTGCATATAATAGTGCAAAGCAAGTACCAGTCAATTAATGTAATTTTTAGTGCAAAATTCGTTTATGTGATCGTCCGGGAGTTTTCCGAGTTGACTTTCAGTTCTGGGCTTAACTGTcaaa                       <  2:1030555/144‑1 (MQ=255)
                                                                                                                    tgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCg                                                                   >  1:184126/1‑104 (MQ=255)
                                                                                                                    tgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCg                                                                   <  2:184126/104‑1 (MQ=255)
                                                                                                                     gcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGGAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTGCAGTTCTGGGCTTAACTGTcaaagc                     >  1:997760/1‑149 (MQ=255)
                                                                                                                     gcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagc                     <  1:277602/149‑1 (MQ=255)
                                                                                                                       aTTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcacagcaa                   >  2:36141/1‑149 (MQ=255)
                                                                                                                        ttCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagcaaa                  <  1:1003/149‑1 (MQ=255)
                                                                                                                           tctGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGat               >  2:845945/1‑149 (MQ=255)
                                                                                                                               aGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGataata           <  1:449544/149‑1 (MQ=255)
                                                                                                                               aGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGataata           <  1:63566/149‑1 (MQ=255)
                                                                                                                                gCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATAc          <  1:224600/149‑1 (MQ=255)
                                                                                                                                      atTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAACAATCTTTTATTTTATCGTTCGGGATTATCACTAGAAGACGTGCAATGCTGGGCTTAACTGTCAAAGCAAAGATAATACcactca    >  2:1061618/1‑149 (MQ=255)
                                                                                                                                      atTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactca    >  1:752864/1‑149 (MQ=255)
                                                                                                                                      atTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactca    >  1:752896/1‑149 (MQ=255)
                                                                                                                                       ttAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATGTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAACGCAAAGATAATACcactcac   >  2:1039691/1‑149 (MQ=255)
                                                                                                                                        tAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactcact  <  1:36141/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                              
ATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACT  >  NC_000913_3_pae_tpiA/1197360‑1197644

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