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,5020CT98.0% 157.9 / NA 52G11R (GGA→AGA) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (27/23);  minor base G (1/0);  total (29/23)
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

ACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGA  >  NC_000913_3_pae_tpiA/1197356‑1197646
                                                                                                                                                  |                                                                                                                                                
aCTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                                                                                                                <  1:671051/149‑1 (MQ=255)
aCTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAAGATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                                                                                                                <  2:314002/149‑1 (MQ=255)
  tGATCAATTCATCAGAATATCTAGGTTTTCCGCAGCCAAAATCAACGGCGTGGCCATTTTTCGCAATGCTTCTAATATGCTCGCAAAGATAGCTGGAGAGCATCGTATGACGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                                                                                                              <  1:831940/149‑1 (MQ=255)
    aTTAATTCGTCAGAATATCCAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATGCGCGCAAAGATATCTTGTGGGCATCGGATGAGGTTTTGCTGCATACTATCAGCGAATATTAACTCTGTGCaa                                                                                                                                            >  2:831940/1‑148 (MQ=255)
    aTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCGATTTTTCTCAATGCTTCTAATATACTAGCAAAGATATCTAGAGGGCATAGTATGCGGTTTTGCTGAATTTTCTGTGCGAGTGTTACATCTGTGCat                                                                                                                                            >  2:590373/1‑149 (MQ=255)
         ttCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtattat                                                                                                                                       >  2:603460/1‑149 (MQ=255)
             tcaGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGtt                                                                                                                                   <  2:943441/149‑1 (MQ=255)
                gAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCACTGCTTATAATATACTCACAAAGAAATCTAGAGCGCATCGTATGATGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTTTGCATATTATAGTTCaa                                                                                                                                >  2:581212/1‑149 (MQ=255)
                   tatCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCAACGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCATGCATATTATAGTTCAAAGc                                                                                                                             >  2:875299/1‑149 (MQ=255)
                     tCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCaa                                                                                                                           >  1:314946/1‑149 (MQ=255)
                       tAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCCAAGATATCTGGATGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGCATATTAACTCGCTGCATATTATAGTTCACAGCAAGt                                                                                                                         >  2:694726/1‑149 (MQ=255)
                           ttttCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCa                                                                                                                     >  1:694487/1‑149 (MQ=255)
                                cGCAGCCAAAATCAAGGGCGTGGGCATTTTTCGCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTACACTTGAAGGCAAGTACCATTCaa                                                                                                                <  2:611654/149‑1 (MQ=255)
                                 gCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAAt                                                                                                               >  1:907253/1‑149 (MQ=255)
                                      cAAAATCAAGGGCGTGGCCATTTTTCTCAATGTTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATg                                                                                                          <  2:96738/149‑1 (MQ=255)
                                      cAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATg                                                                                                          >  1:951252/1‑149 (MQ=255)
                                                 gCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGt                                                                                               <  1:861675/149‑1 (MQ=255)
                                                       ccATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGaaaaa                                                                                         <  1:512699/149‑1 (MQ=255)
                                                         aTTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTTTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTTTTTTTAGTGAAAAAtt                                                                                       >  2:922047/1‑149 (MQ=255)
                                                                tcAATGCTTCTAATATAGTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTATTGTGAAAGTCttttat                                                                                <  2:612539/149‑1 (MQ=255)
                                                                  aaTGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttattt                                                                              >  2:797503/1‑149 (MQ=255)
                                                                    tGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGTTGCATTGTGTGAGCGAATATTTGCTCTGTGCATATTATAGTTCAAAGCAGGTAGCTTTAAATTGTTGTTATTTTTAGTGAAGAATTGttttattttt                                                                            >  2:678833/1‑148 (MQ=255)
                                                                    tGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttatttta                                                                            <  2:76618/149‑1 (MQ=255)
                                                                      cTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATc                                                                          <  1:709999/149‑1 (MQ=255)
                                                                          tAATTTATTCGCTATGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTcc                                                                      <  1:155053/149‑1 (MQ=255)
                                                                          tAATATACTCGCAAAGACACCGAGAGGGCATCGTATGAGGTTCTGCTGCATACTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCACTCAATAAATGTTACTTTTAGTGAAAAATCCATTTATTTTATctccg                                                                      >  2:854788/1‑145 (MQ=255)
                                                                           aacataCTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCg                                                                     <  1:663847/146‑1 (MQ=255)
                                                                           aaTATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTACg                                                                     >  1:213664/1‑149 (MQ=255)
                                                                                  tCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTc                                                                                                                                  >  2:843365/1‑81 (MQ=255)
                                                                                  tCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTc                                                                                                                                  <  1:843365/81‑1 (MQ=255)
                                                                                  tCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGAtttt                                                              >  1:704652/1‑149 (MQ=255)
                                                                                    gCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTc                                                            >  1:789727/1‑149 (MQ=255)
                                                                                     cAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGAGATTTTTcc                                                           <  2:803959/149‑1 (MQ=255)
                                                                                     cAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCAGTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTcc                                                           <  1:197180/149‑1 (MQ=255)
                                                                                          atatCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCAGGGATTTTTCCtactt                                                      >  2:955660/1‑149 (MQ=255)
                                                                                                           tATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGc                                     >  2:77289/1‑149 (MQ=255)
                                                                                                                  ttttGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTg                              >  2:857235/1‑149 (MQ=255)
                                                                                                                        tgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACGATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGTATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaag                        >  1:943199/1‑149 (MQ=255)
                                                                                                                        tgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTGAACTAATCTTATGTTTAGTGAAAAGTTGGTGTATTTTATCGTCCTTGATTTTTCCGACTTGACTTTCAGTTCTGGGATTAACTGTcaaag                        >  1:445426/1‑149 (MQ=255)
                                                                                                                        tgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATGTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaag                        >  1:575573/1‑149 (MQ=255)
                                                                                                                           aTTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagcaa                     >  1:880698/1‑149 (MQ=255)
                                                                                                                            ttCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcacaataaa                    <  2:907253/149‑8 (MQ=255)
                                                                                                                             tctcTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagcaaag                   <  1:77289/149‑1 (MQ=255)
                                                                                                                               tctGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTaa                                 >  2:931541/1‑133 (MQ=255)
                                                                                                                               tctGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTaa                                 <  1:931541/133‑1 (MQ=255)
                                                                                                                               tctGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGat                 >  2:280422/1‑149 (MQ=255)
                                                                                                                                 tGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGataa               >  1:269076/1‑149 (MQ=255)
                                                                                                                                  gatcGAATATTAACTCTGTGCATAATATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGataat              <  1:95761/146‑1 (MQ=255)
                                                                                                                                        atatTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcact        <  1:199551/149‑1 (MQ=255)
                                                                                                                                         tatTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACGGTCAAAGCAAAGATAATACcactc       >  1:979748/1‑149 (MQ=255)
                                                                                                                                          atTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTATTGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGTTTTTTCCGACTTGACTTTCAGTTCTGTGCTTAACTGTCAAAGCAAAGATCATACcactcg      <  2:390452/149‑2 (MQ=255)
                                                                                                                                              aCTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTga  <  1:857235/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGA  >  NC_000913_3_pae_tpiA/1197356‑1197646

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