Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,427,156 T→G 9.9% R314R (CGA→CGC)  USA300HOU_2291 ← formate dehydrogenase alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,427,1560TG 9.9% 95.1 / 4.1 38R314R (CGA→CGC) USA300HOU_2291formate dehydrogenase alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (15/19);  new base G (2/2);  total (17/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 = 9.86e-01

AGAAATGAATGATAGTGCATCTGGGCCATGCTTTTCTTTAATAGCTGTAAAATTATCTGCAATGACGTTTAAAGCTTCATCCCATTCTACTTCATGGAACTCACCATTTTTCCTTACTAGTGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATGGTTGTACTTTTAAAATTTCTCTATCTTTAGTCCAAACTTCAAATGAACAACCCACACC  >  CP000730/2427024‑2427296
                                                                                                                                    |                                                                                                                                            
agaaATGAATGATAGTGCATCTGGGCCATGCTTTTCTTTAATAGCTGTAAAATTATCTGCAATGACGTTTAAAGCTTCATCCCATTCTACTTCATGGAACTCACCATTTTTCCTTACTAGTGGTTTAGTTAATCGTTGATc                                                                                                                                      <  1:44179/141‑1 (MQ=255)
agaaATGAATGATAGTGCATCTGGGCCATGCTTTTCTTTAATAGCTGTAAAATTATCTGCAATGACGTTTAAAGCTTCATCCCATTCTACTTCATGGAACTCACCATTTTTCCTTACTAGTGGTTTAGTTAATCGTTGATc                                                                                                                                      <  1:35817/141‑1 (MQ=255)
  aaaTTAATGATATTGCATCTGGGCCTTACTTTTCTTTAGTATATGTAAAGTTATCTGCAATGACATTTAAATCTTCATCCCATTCTACTTCATGGAACTCACCATTTTTCCTTACTAGTGGTTTAGTTAAGCGTTGATCTg                                                                                                                                    <  1:19871/141‑1 (MQ=255)
         tgaTAGTGCATCTGGGCCATGCTTTTCTTTAATAGCTGTAAAATTATCTGCAATGACGTTTAAAGCTTCATCCCATTCTACTTCATGGAACTCACCATTTTTCCTTACTAGTGGTTTAGTTAATCGTTGATCTGAATTAat                                                                                                                             <  2:379080/141‑1 (MQ=255)
               tGCATCTGGGCCATGCTTTTCTTTAATAGCTGTAAAATTATCTGCAATGACGTTTAAAGCTTCATCCCATTCTACTTCATGGAACTCACCATTTTTCCTTACTAGTGGTTTAGTTAATCGTTGATCTGAATTAATATGTcc                                                                                                                       >  2:403188/1‑141 (MQ=255)
                             gCTTTTCTTTAATAGCTGTAAAATTATCTGCAATGACGTTTAAAGCTTCATCCCATTCTACTTCATGGAACTCACCATTTTTCCTTACTAGTGGTTTAGTTAATCGTTGATCCGAATTAATATGTCCCCATGAAAACTTAc                                                                                                         <  2:178958/141‑1 (MQ=255)
                               ttttCTTTAATAGCTGTAAAATTATCTGCAATGACGTTTAAAGCTTCATCCCATTCTACTTCATGGAACTCACCATTTTTCCTTACTAGTGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCt                                                                                                       <  2:192065/141‑1 (MQ=255)
                                 ttCTTTAATAGCTGGAAAATTATCTGCTATGACGTTTAAAGCTTCATCCCATTCTTCTGCATGGAACTCACCATTTGGCCTTTCTATTGTGTTAGTTAATCGTTGATCTGAATTAGTATATCCCCATTAAAACTTACCttt                                                                                                     <  2:318936/141‑1 (MQ=255)
                                   cTTTAATAGCTGTAAAATTATCTGCAATGACGTTTAAAGCTTCATCCCATTCTACTTCATGGAACTCACCATTTTTCCTTACTAGTGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTaa                                                                                                   >  1:201475/1‑141 (MQ=255)
                                         tAGCTGTAAAATTATCTGCAATGACGTTTAAAGCTTCATCCCATTCTACTTCATGGAACTCACCATTTTTCCTTACTAGTGGTTTAGTTAATCGTTGATCTGAATTAAAATGTCCCCATGAAAACTTACCTTTAACACAAg                                                                                             <  1:141124/141‑1 (MQ=255)
                                                   aTTATCTGCAATGACGTTTAAAGCTTCATCCCATTCTACTTCATGGAACTCACCATTTTTCCTTACTAGTGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAAtttt                                                                                   <  1:334873/141‑1 (MQ=255)
                                                      aTCTGCAATGACGTTTAAAGCTTCATCCCATTCTACTTAATGGAACTCACCATTTTTCCTCACTAAGGGGTTAACCACTGGTGGCTGGGAATTAATAAGTGCCCATGAAAACCTACCTTTAACACAAGTCGCAATTTTAtt                                                                                >  1:154532/1‑141 (MQ=255)
                                                                     tAAAGCTTCATCCCATTCTACTTCATGGAACTCACCATTTTTCCTTACTAGTGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGAAATTTTATTTGCTGGAGAATCAtt                                                                 >  1:227443/1‑140 (MQ=255)
                                                                      aaaGCTTCATCCCATTCTACTTCATGGAACTCACCATTTTTCCTTACTAGTGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCAtgt                                                                >  1:50275/1‑141 (MQ=255)
                                                                         gCTTCATCCCATTCTACTTCATGGAACTCACCATTTTTCCTTACTAGTGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGAt                                                             >  1:375724/1‑141 (MQ=255)
                                                                          ttttagcccATTCTACGTGCTGGAACTCAGCGTTTTCCCTTACTTGAGGTTTAGTTAATCGTTGTTCTGAATTAATATGTCCCCTTGAAAACTTACCTTTGGCACCGGTCGCAATTTTATTTGCTGAAGAATCATGTGATg                                                            <  2:393288/135‑1 (MQ=255)
                                                                           ttCATCCCATTCTACTTCATGGAACTCACCATTTTTTTTTTTTTTTTTTTTAGTTAAGCGTTTATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATgg                                                           >  2:192944/1‑141 (MQ=255)
                                                                                cccATTCTACTTCATGGAACTCACCATTTTTCCTTACTAGTGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATGGTTGTa                                                      <  1:190908/141‑1 (MQ=255)
                                                                                cccATTCTACTTCATGGAACTCACCATTTTTCCTTACTAGTGGTTTAGTTAAGCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATGGTTGTa                                                      <  1:136714/141‑1 (MQ=255)
                                                                                    ttCTACTTCATGGAACTCACCATTTTTCCTTACTAGTGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATGGTTGTACttt                                                  <  1:79829/141‑1 (MQ=255)
                                                                                    ttCTACTTCATGGAACTCACCATTTTTCCTTACTAGTGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATGGTTGTACttt                                                  >  1:404746/1‑141 (MQ=255)
                                                                                               ggAACTCACCATTTTTCCTTACTAGTGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATGGTTGTACTTTTAAAATTtctc                                       <  2:290938/141‑1 (MQ=255)
                                                                                                   cTCACCATTTTTCCTTACTAGTGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATGGTTGTACTTTTAAAATTTCTCTATc                                   >  1:288277/1‑141 (MQ=255)
                                                                                                   cTCACCATTTTTCCTTACTAGTGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATGGTTGTACTTTTAAAATTTCTCTATc                                   >  1:402490/1‑141 (MQ=255)
                                                                                                   cTCACCATTTTTCCTTACTAGTGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATGGTTGTACTTTTAAAATTTCTCTATc                                   >  1:402488/1‑141 (MQ=255)
                                                                                                      ccgaTTTTTCCTTGCTAGTTGTTTAGTTGATCGTTGATATCAAGTAATATGTGCGCATGAAATCTTAACATTAACACAAGTCGCTATTTTATTTGCTGGAGAATCATGTGATGGAGGTACTATTAAAACATCTCTATCttg                                <  2:397770/138‑2 (MQ=255)
                                                                                                         aTTTTTCCTTACTAGTGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATGGTTGTACTTTTAAAATTTCTCTATCTTTAGt                             <  2:124110/141‑1 (MQ=255)
                                                                                                            tttCCTTACTAGTGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATGGTTGTACTTTTAAAATTTCTCTATCTTTAGTCCa                          <  1:33457/141‑1 (MQ=255)
                                                                                                             ttCCTTACTAGTGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATGGTTGTACTTTTAAAATTTCTCTATCTTTAGTCCaa                         <  2:346814/141‑1 (MQ=255)
                                                                                                               tcTTATTATTGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATGGTTGTACTTTTAAAATTTCTCTATCTTTAGTCCAAAc                       >  2:121420/2‑141 (MQ=255)
                                                                                                               ccTTACTAGTGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATGGTTGTACTTTTAAAATTTCTCTATCTTTAGTCCAAAc                       >  1:121390/1‑141 (MQ=255)
                                                                                                                  tACTAGTGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAGTTTTATTTGCTGGAGAATCATGTGATGGTTGTACTTTTAAAATTTCTCTATCTTTAGTCCAAACTTc                    >  1:445313/1‑141 (MQ=255)
                                                                                                                     tAGTGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATGGTTGTACTTTTAAAATTTCTCTATCTTTAGTCCAAACTTCaaa                 <  1:420787/141‑1 (MQ=255)
                                                                                                                        tGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATGGTTGTACTTTTAAAATTTCTCTATCTTTAGTCCAAACTTCAAATGa              >  1:279851/1‑141 (MQ=255)
                                                                                                                          gTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATGGTTGTACTTTTAAAATTTCTCTATCTTTAGTCCAAACTTCAAATGaac            <  2:246779/141‑1 (MQ=255)
                                                                                                                              aGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATGGTTGTACTTTTAAAATTTCTCTATCTTTAGTCCAAACTTCAAATGAACAAcc        <  1:440979/141‑1 (MQ=255)
                                                                                                                              aGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATGGTTGTACTTTTAAAATTTCTCTATCTTTAGTCCAAACTTCAAATGAACAAcc        >  1:293924/1‑141 (MQ=255)
                                                                                                                                ttAAGCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATGGTTGTACTTTTAAAATTTCTCTATCTTTAGTCCAAACTTCAAATGAACAACcca      >  1:127550/1‑141 (MQ=255)
                                                                                                                                    tCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATGGTTGTACTTTTAAAATTTCTCTATCTTTAGTCCAAACTTCAAATGAACAACccacacc  >  2:298264/1‑141 (MQ=255)
                                                                                                                                    |                                                                                                                                            
AGAAATGAATGATAGTGCATCTGGGCCATGCTTTTCTTTAATAGCTGTAAAATTATCTGCAATGACGTTTAAAGCTTCATCCCATTCTACTTCATGGAACTCACCATTTTTCCTTACTAGTGGTTTAGTTAATCGTTGATCTGAATTAATATGTCCCCATGAAAACTTACCTTTAACACAAGTCGCAATTTTATTTGCTGGAGAATCATGTGATGGTTGTACTTTTAAAATTTCTCTATCTTTAGTCCAAACTTCAAATGAACAACCCACACC  >  CP000730/2427024‑2427296

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