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

TAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCA  >  NC_000913_3_pae_tpiA/1197362‑1197642
                                                                                                                                            |                                                                                                                                             
tAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCat‑at                                                                                                                                      >  1:47691/1‑149 (MQ=255)
tAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCAAGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAAGATTATCTCTGTGAat‑at                                                                                                                                      >  1:274794/1‑149 (MQ=255)
   ttCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCTTTTGAGGTTTTGCTGCATTCTCTGAGCGAATATTAATTTTGTGCTt‑attat                                                                                                                                   >  2:54814/1‑149 (MQ=255)
   ttCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAt‑attat                                                                                                                                   >  1:637430/1‑149 (MQ=255)
         aGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCa                                                                                                                             >  2:153877/1‑149 (MQ=255)
          gAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGCTATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAAGATTAACTCTGGACGT‑ATTATAGTTCaa                                                                                                                            >  2:444315/1‑149 (MQ=255)
              atCTAAGTTTTACGCAGCCAAAATAAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGATGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCa                                                                                                                        <  2:152084/149‑1 (MQ=255)
                 tAAGTTTACCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTACTAATATAATCGCAAAGGTGTCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGt                                                                                                                     >  1:241172/1‑149 (MQ=255)
                               gcAAAATCAAGGGCGGCGCGATTTTTCTCAATGCCTCTAATATACTCGCAAAGATAACTAGAGGGCATCGTATGAGGTTTTGCTGCATTATCTGAGCGAATATTCA‑GCCTTTCATCCTTTCAGGTCAAAGCAAGTACCATTCAATTAAt                                                                                                       >  2:319803/2‑149 (MQ=255)
                                cAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATg                                                                                                      >  1:499353/1‑149 (MQ=255)
                                      cAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTAttt                                                                                                >  1:57744/1‑149 (MQ=255)
                                       aaGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTAtttt                                                                                               <  1:153877/149‑1 (MQ=255)
                                           gtgTGGCCTTTTTTCTCAATTCTTCGAATAGATTCGCAAAGAGATCTAGAGGGTATCGTATGAGGTTTTGCTGCATTCTCTTAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGt                                                                                           <  2:198418/147‑1 (MQ=255)
                                                 caaTTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTCTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGaaaaa                                                                                     <  1:140592/147‑1 (MQ=255)
                                                   aTTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAAtt                                                                                   >  2:7787/1‑149 (MQ=255)
                                                      tttCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCtt                                                                                <  2:499353/149‑1 (MQ=255)
                                                             aTGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttatttt                                                                         >  1:658582/1‑149 (MQ=255)
                                                                   cTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTc                                                                   <  2:179071/149‑1 (MQ=255)
                                                                    tAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTATATCGTcc                                                                  >  2:68568/1‑149 (MQ=255)
                                                                        ataCTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATCTTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGa                                                              >  2:204462/1‑149 (MQ=255)
                                                                            tCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATACTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTCTCGTCCGGGAtttt                                                          >  2:187894/1‑149 (MQ=255)
                                                                                   gATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGTTTTTTCCgcct                                                   >  2:131503/1‑149 (MQ=255)
                                                                                                          ggTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAAc                            <  1:121004/149‑1 (MQ=255)
                                                                                                           gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTaa                             <  1:2487/147‑1 (MQ=255)
                                                                                                           gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTaa                             >  2:2487/1‑147 (MQ=255)
                                                                                                           gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACt                           <  1:575751/149‑1 (MQ=255)
                                                                                                                 gtaaaatctcTGAGCAAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACGATTCAATTAATGTTATTTTGAGTGAAAAATTCTTTTATTTTATCGCCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaa                     <  1:67002/143‑1 (MQ=255)
                                                                                                                  tgcATTCTCTGAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTCTCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaag                    >  1:638659/1‑149 (MQ=255)
                                                                                                                            gAGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGataat          >  1:224478/1‑149 (MQ=255)
                                                                                                                             aGCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGataata         >  2:106497/1‑149 (MQ=255)
                                                                                                                              gCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATc                                                                      >  2:21317/1‑87 (MQ=255)
                                                                                                                              gCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATc                                                                      <  1:21317/87‑1 (MQ=255)
                                                                                                                              gCGAATATTAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATAc        <  1:7787/149‑1 (MQ=255)
                                                                                                                                    actAACTCTGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGAGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactca  <  1:346466/147‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                             
TAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCAT‑ATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCA  >  NC_000913_3_pae_tpiA/1197362‑1197642

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: ‑