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,5020CT95.7% 132.8 / NA 48G11R (GGA→AGA) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (26/19);  minor base A (1/1);  total (27/21)
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

CTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAG  >  NC_000913_3_pae_tpiA/1197357‑1197647
                                                                                                                                                 |                                                                                                                                                  
cTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCtg‑tg                                                                                                                                                >  1:316272/1‑149 (MQ=255)
   aTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAAGGCATCGTATAAGGTTTTGCTTCATTCTCTGAGCGAATATTAACTCTG‑TGCat                                                                                                                                             >  2:636129/1‑149 (MQ=255)
     tAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTTAGCGAATATTAACTGTG‑TGCatat                                                                                                                                           >  1:740866/1‑149 (MQ=255)
        ttCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCAtattat                                                                                                                                        >  1:294900/1‑149 (MQ=255)
           atcaGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGt                                                                                                                                     >  2:589017/1‑149 (MQ=255)
            tcaGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTACATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGtt                                                                                                                                    >  1:42030/1‑149 (MQ=255)
              aGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCa                                                                                                                                  >  1:396349/1‑149 (MQ=255)
               gAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCaa                                                                                                                                 >  1:123400/1‑149 (MQ=255)
                    tCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCaa                                                                                                                            <  2:215264/149‑1 (MQ=255)
                      tAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGt                                                                                                                          <  1:306267/149‑1 (MQ=255)
                          ttttCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCa                                                                                                                      >  2:146632/1‑149 (MQ=255)
                                    ccAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAt                                                                                                            <  1:572785/149‑1 (MQ=255)
                                    ccAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTAGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAt                                                                                                            >  1:669959/1‑149 (MQ=255)
                                    ccAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTAGCAAAGATATCTAAAGGGCATCGTATGAGGTTTTGCTGCATTATCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTAAATTAAt                                                                                                            >  2:56405/1‑149 (MQ=255)
                                     cAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCTTATGAGTTTTTGCTGCATTCTCTGAGCGAATANTAACTCTG‑TGAATATTATAGTNTNAAGCANTNACCATTCAATTAAtt                                                                                                           >  1:167658/1‑148 (MQ=255)
                                          tCAAGAGCGGCGCCATTTTTCTCAATGCTTATAATAGACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTAtt                                                                                                      <  2:642931/149‑1 (MQ=255)
                                               ggagtagCCATTTTTCTCAAAGCATCAAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAg                                                                                                 <  1:971883/143‑1 (MQ=255)
                                               ggCGTGGCCATTTTTCTTAGTGCTTCTAATATACTCGACAAGATCTCTGGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGTGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAg                                                                                                 <  1:611214/149‑1 (MQ=255)
                                                   tGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAAAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTAACAACAAATTAATGTTATTTTTAGTGaa                                                                                             >  1:281384/1‑149 (MQ=255)
                                                             tctcAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCtttt                                                                                   <  2:387293/149‑1 (MQ=255)
                                                              ctcAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCtttta                                                                                  <  1:146632/149‑1 (MQ=255)
                                                               tcAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttat                                                                                 <  2:855231/149‑1 (MQ=255)
                                                                  aTGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttatttt                                                                              >  2:421841/1‑149 (MQ=255)
                                                                       tCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTAAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGt                                                                         >  1:179/1‑149 (MQ=255)
                                                                           atatACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCgg                                                                     >  1:871003/1‑149 (MQ=255)
                                                                                  cGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGAttttt                                                              <  2:294900/149‑1 (MQ=255)
                                                                                   gaaaaGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCG‑ACAACAACCCCGTTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTc                                                             <  1:240062/147‑1 (MQ=255)
                                                                                       aGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgac                                                         >  1:106181/1‑149 (MQ=255)
                                                                                        gATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgact                                                        <  1:920561/149‑1 (MQ=255)
                                                                                               agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACttt                                                 <  2:316272/149‑1 (MQ=255)
                                                                                                      aTCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCAGAGATTTTTTCGAGTTGACTTTCACTTCt                                          >  2:277390/1‑149 (MQ=255)
                                                                                                      aTCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCAG‑TGCATATTATAGTTCAAAACAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTATTTTATTTTATCGTCCGTGATTTTTCCGACTTGACttttattgct                                          >  1:93043/1‑142 (MQ=255)
                                                                                                       tCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTg                                         <  1:421841/149‑1 (MQ=255)
                                                                                                        cGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTgg                                        <  1:653238/149‑1 (MQ=255)
                                                                                                         ttagtaGGTCTTGCTGCATTCTCTTAGCGAACATTAACTATG‑TGCATATTATAGTTCCAAGCAAGTACCATTCAATGCATGCAATTTTGAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTggg                                       <  2:787891/144‑1 (MQ=255)
                                                                                                           aTGAGGTTTTGCTGCATTCACTTAGCGTATATTAACTCAG‑TGCATCGGCGAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCt                                     <  1:923645/149‑1 (MQ=255)
                                                                                                             gAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTa                                   >  2:814376/1‑149 (MQ=255)
                                                                                                             gAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTa                                   >  1:675779/1‑149 (MQ=255)
                                                                                                                gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAGGTACCATTCAATTAATGTTATTTTTAGTGACTAATATGATTAGTTTATCGTCCGGGTTTTTTCCGACTTGAATTTCAATTCTGGGCTTAACt                                >  1:764994/1‑149 (MQ=255)
                                                                                                                 ttttGCTGCATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACtt                               >  2:1025494/1‑148 (MQ=255)
                                                                                                                    tgctgcATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTca                            <  2:123400/149‑1 (MQ=255)
                                                                                                                      ctgcATTCTCTGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaa                          >  2:733568/1‑149 (MQ=255)
                                                                                                                               ctGAGCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGata                 <  2:605281/149‑1 (MQ=255)
                                                                                                                                   gCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCa                                               <  2:267977/115‑1 (MQ=255)
                                                                                                                                   gCGAATATTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCa                                               >  1:267977/1‑115 (MQ=255)
                                                                                                                                     gAATATTGACGCTG‑TGCACTTTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGGGAAAAATTCTTTTATTTTCTCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACca           <  2:3668/149‑1 (MQ=255)
                                                                                                                                        tatTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactc        >  1:61210/1‑149 (MQ=255)
                                                                                                                                         atTAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGAGTTTTCGTTCTGGACATAACTGTCAAAGAAAAAATAATACcactct       >  1:810735/1‑148 (MQ=255)
                                                                                                                                           tAACTCTG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactcact     >  1:526885/1‑149 (MQ=255)
                                                                                                                                              ctctG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTGTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTgag  >  1:198446/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCG‑TGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAG  >  NC_000913_3_pae_tpiA/1197357‑1197647

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