Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_vch_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_vch_tpiA1,197,5020CT97.1% 113.6 / ‑5.9 35G11R (GGA→AGA) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (22/12);  minor base G (1/0);  total (23/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACT  >  NC_000913_3_vch_tpiA/1197356‑1197644
                                                                                                                                                  |                                                                                                                                               
aCTGAATAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                                                                                                               <  2:695524/149‑1 (MQ=255)
 cTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTATAGGGCATCGTCTGATCTGGTGCTGCATTCTCTGAGCGAATATTAACTCtgtg                                                                                                                                              <  1:385242/149‑1 (MQ=255)
 cTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCtgtg                                                                                                                                              <  1:351932/149‑1 (MQ=255)
  tGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAAAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                                                                                                             >  1:442995/1‑149 (MQ=255)
            atcaGAATATCTAAGTTTTCCGCTGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGt                                                                                                                                   <  2:157724/149‑1 (MQ=255)
                     tCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCaa                                                                                                                          <  1:467775/149‑1 (MQ=255)
                         aGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTAc                                                                                                                      <  1:294834/149‑1 (MQ=255)
                             ttCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCAtt                                                                                                                  <  1:42908/149‑1 (MQ=255)
                                           tCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAAGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTAAATTAATGTTAtt                                                                                                    >  1:320901/1‑149 (MQ=255)
                                                   gTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCACTTAATGTTATTTTTATTGa                                                                                            >  1:14865/1‑149 (MQ=255)
                                                      gCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATTACGGTATTTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGaaaa                                                                                         >  1:741192/1‑149 (MQ=255)
                                                               ctcAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTATTGAAAAATTCtttta                                                                                >  2:542761/1‑149 (MQ=255)
                                                                     gCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTtttattttat                                                                          >  2:463197/1‑149 (MQ=255)
                                                                          tAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTcc                                                                     >  1:597165/1‑149 (MQ=255)
                                                                              ataCTCGCAAATATATCTAGAGTGCATCGTATCAGGTTTTGCTGAATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGa                                                                 <  1:454157/149‑1 (MQ=255)
                                                                                  tCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑tttt                                                            >  2:246464/1‑149 (MQ=255)
                                                                                    gCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTc                                                          >  1:617014/1‑149 (MQ=255)
                                                                                       aaGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATAGTCCGGGA‑TTTTTCCgc                                                       >  2:443140/1‑148 (MQ=255)
                                                                                                    ggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGt                                          >  1:93572/1‑149 (MQ=255)
                                                                                                      cATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCGGTGCATATTATAGTTCAAAGCAAGTACCAGTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCAGACTTGACTTTCAGTTc                                        >  2:589786/1‑149 (MQ=255)
                                                                                                        tCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCTGGA‑TTTGTCCGACTTGACTTTCAGTTCTg                                      <  2:268371/149‑1 (MQ=255)
                                                                                                                  ttttGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGTAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTg                            <  1:265296/149‑1 (MQ=255)
                                                                                                                   tttGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCAGACTTGACTTTCAGTTCTGGGTTTAACtat                           >  2:183069/1‑147 (MQ=255)
                                                                                                                       ctgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaa                       >  2:654038/1‑149 (MQ=255)
                                                                                                                       ctgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaa                       <  2:85926/149‑1 (MQ=255)
                                                                                                                           aTTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCATTTCTATGCTTAACTTTaaaagcaa                   >  2:73653/1‑149 (MQ=255)
                                                                                                                               tctGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCg                                                                        <  1:29288/93‑1 (MQ=255)
                                                                                                                               tctGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCg                                                                        >  2:29288/1‑93 (MQ=255)
                                                                                                                               tctGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGa                >  1:695879/1‑149 (MQ=255)
                                                                                                                               tctGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGat               <  2:617014/149‑1 (MQ=255)
                                                                                                                                 tGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGataa             >  1:245051/1‑149 (MQ=255)
                                                                                                                                   aGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCACTTAATGTTATTTTTAGTGACAACTTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGTCTTGACTTTCAGTTCTGGGCTTAACTGTCACAGCACAGatcata           >  1:554157/1‑149 (MQ=255)
                                                                                                                                    gCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATAc          >  2:525931/1‑149 (MQ=255)
                                                                                                                                    gCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATAc          <  1:46030/149‑1 (MQ=255)
                                                                                                                                    gCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATAc          >  1:154863/1‑149 (MQ=255)
                                                                                                                                         tatTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactc     >  1:310897/1‑149 (MQ=255)
                                                                                                                                            tAACTATGTACAAATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactcact  >  2:291968/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
ACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGA‑TTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACT  >  NC_000913_3_vch_tpiA/1197356‑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: ‑