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,5020CT96.5% 179.6 / NA 57G11R (GGA→AGA) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (32/23);  minor base G (1/1);  total (33/24)
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

TGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACT  >  NC_000913_3_pae_tpiA/1197358‑1197644
                                                                                                                                                |                                                                                                                                               
tGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTTATATACTCGCAAAGATATCTAGAGGGCATAGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTATGTGc                                                                                                                                             >  1:221048/1‑149 (MQ=255)
  aTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCat                                                                                                                                           <  1:80726/149‑1 (MQ=255)
         catcaGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAg                                                                                                                                    <  2:103698/149‑1 (MQ=255)
           tcaGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGtt                                                                                                                                  >  1:1356090/1‑149 (MQ=255)
               aaTATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTACCTGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCaaa                                                                                                                              <  1:1359361/149‑1 (MQ=255)
                atatCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCGATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCACCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGAATATAATAGTTCAAAg                                                                                                                             <  2:1307507/149‑1 (MQ=255)
                   tcgtaGTTATCTGCAACCAAAGTCACGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCGCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCaa                                                                                                                          <  2:877129/145‑1 (MQ=255)
                   tCTAAGTTTTCCGGAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATTCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCaa                                                                                                                          <  2:1227939/149‑1 (MQ=255)
                    cTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAg                                                                                                                         >  1:25455/1‑149 (MQ=255)
                      aaGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGTTTACAATGCGCTCACAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTa                                                                                                                       <  1:1132033/149‑1 (MQ=255)
                         ttttCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCa                                                                                                                    >  2:97251/1‑149 (MQ=255)
                         ttttCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCa                                                                                                                    >  1:872347/1‑149 (MQ=255)
                           ttCCGCAGCCAAAATCAAGGGCGGGGCCATTTTTCTCAATGCTTCTAATATGCTCGCAAAGAGATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCAtt                                                                                                                  <  2:582117/149‑1 (MQ=255)
                              cGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAg                                                                                                                                    >  1:541573/1‑128 (MQ=255)
                              cGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAg                                                                                                                                    <  2:541573/128‑1 (MQ=255)
                               gCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAATTACCATTCAAt                                                                                                              >  1:463023/1‑149 (MQ=255)
                                   ccAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGCGAGTGCCAGGACGCGGTGCATCTTTTGGTTGAAAGCAAGTACCATTCAATTAAt                                                                                                          >  1:1413867/1‑149 (MQ=255)
                                      aaaTCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCAAGCGCTAAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGtt                                                                                                       >  1:701401/1‑149 (MQ=255)
                                           aaGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAACAAAGTACCATTCAATTAATGTTAtttt                                                                                                  >  2:20967/1‑149 (MQ=255)
                                                    gCCATTTTTCTCAATTCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGaaaa                                                                                         >  2:919530/1‑149 (MQ=255)
                                                      cATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGaccagt                                                                                       >  1:1051110/1‑144 (MQ=255)
                                                       aTTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGATGCATTCTCTGAACGAATATTAACTCTATGCATATTATAGTTCAAAACAAATACCAATCAATTAATATTATTTTTAGTAAAACAtt                                                                                      >  2:629067/1‑149 (MQ=255)
                                                        tttttCTCAATGCTTCTAATATACTCGCAAAGACATCTAGAGGGCATCGTATGAGTTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGCGAAAAATTc                                                                                     <  1:938552/149‑1 (MQ=255)
                                                         ttttCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCt                                                                                    >  2:477511/1‑149 (MQ=255)
                                                          tttCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCtt                                                                                   >  1:718539/1‑149 (MQ=255)
                                                             ctcCATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTCTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCtttta                                                                                <  1:1259306/149‑1 (MQ=255)
                                                              tcAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttat                                                                               >  2:477872/1‑149 (MQ=255)
                                                                   gCTTTTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTtttattttat                                                                          <  2:51850/149‑1 (MQ=255)
                                                                     ttCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATCATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTATTTTATTTTATCg                                                                        >  2:1161316/1‑149 (MQ=255)
                                                                      tCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGt                                                                       >  1:128870/1‑149 (MQ=255)
                                                                       cTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCAGTCGATTAATGTTATTTTTAGTCAAAAATTCTTTTATTTTATCGTc                                                                      <  1:1146566/149‑1 (MQ=255)
                                                                        tAATATACTCGTAAATATATCTAGATGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAAATCTGTGCATATCATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTATTAGTGAAAAATTCTTTTATTTTATCGTcc                                                                     <  2:806242/149‑1 (MQ=255)
                                                                                tCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTGATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑tttt                                                            >  1:919824/1‑149 (MQ=255)
                                                                                  gCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTc                                                          <  2:257638/149‑1 (MQ=255)
                                                                                  gCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTc                                                          <  2:980127/149‑1 (MQ=255)
                                                                                  gCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTc                                                          <  1:1097407/149‑1 (MQ=255)
                                                                                  gCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTc                                                          >  1:116728/1‑149 (MQ=255)
                                                                                   cAAAGATATCTAGAGGGGATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTcc                                                         <  2:459074/149‑1 (MQ=255)
                                                                                      aGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCgac                                                      <  1:636809/149‑1 (MQ=255)
                                                                                      aGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGAGT‑ATTTTcctcc                                                      >  2:1168502/1‑146 (MQ=255)
                                                                                           tCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGTGA‑TTTTTCCgacttgac                                                 >  2:249368/1‑149 (MQ=255)
                                                                                               gagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTc                                             >  1:568998/1‑149 (MQ=255)
                                                                                                agGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAt                                                                                                          <  2:1125317/88‑1 (MQ=255)
                                                                                                agGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAt                                                                                                          >  1:1125317/1‑88 (MQ=255)
                                                                                                agGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCa                                            <  2:25455/149‑1 (MQ=255)
                                                                                                  ggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGt                                          >  1:967300/1‑149 (MQ=255)
                                                                                                    cATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTTTGCTGAGTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTc                                        >  1:706913/1‑149 (MQ=255)
                                                                                                    cATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTc                                        >  1:376959/1‑149 (MQ=255)
                                                                                                       cGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTCTgg                                     <  2:128870/149‑1 (MQ=255)
                                                                                                         tATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTCTGGGc                                   >  1:453883/1‑149 (MQ=255)
                                                                                                               gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTCTTGGCTTAACt                             >  2:386007/1‑149 (MQ=255)
                                                                                                               gTTTTGCTGCATTCTCGTAGCGCATATTAACTCGGTGCATATTCTAGTTAAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACt                             <  1:193773/149‑1 (MQ=255)
                                                                                                                     ctgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTGTAGTAAAAAATTATTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTCTGGGCTTAAcgtttaaa                       >  2:1245516/1‑142 (MQ=255)
                                                                                                                     ctgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACAATACCATTAATGTTATTTTTGGTTAAACATTTTTTGCCGTGAT‑GTCCCTGATTTTTTCCGTCTTGACTTTTATTTTTGTGCTTAACCGTcaaa                       >  2:305962/1‑149 (MQ=255)
                                                                                                                                        atTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactca    <  1:835519/149‑1 (MQ=255)
                                                                                                                                         ttAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACccctcac   >  1:628622/1‑149 (MQ=255)
                                                                                                                                          tAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactcact  >  1:725875/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACT  >  NC_000913_3_pae_tpiA/1197358‑1197644

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