Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 294,093 G→A P236L (CCT→CTT)  yagM ← protein YagM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913294,0930GA94.4% 47.1 / ‑6.1 19P236L (CCT→CTT) yagMprotein YagM
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (9/8);  total (11/8)
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

GAGAGTGATTGCCAGCGATACTTTAACTGGTAAAGTAGTATTTTACCTCTCATAATTGTTTGTCCTTCTTTGTTACCTCTTAGGTTCTGCTCTCTTGTGAGTATGGGTTTATAATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTAT  >  NC_000913/293866‑294337
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
gagagTGATTGCCAGCGATACTTTAACTGGTAAAGTAGTATTTTACCTCTCATAATTGTTTGTCCTTCTTTGTTACCTCTTAGGTTCTGCTCTCTTGTGAGTATGGGTTTATAATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAACGCATAACAATCCCCTTtcttct                                                                                                                                                                                                                                >  1:330488/1‑250 (MQ=255)
      gATTGTCAGCGATACTTTAACTGGTAAAGTAGTATTTGACCTCTCATAATTGTTTGTCCTTCTTTTTTACCTCTTAGGTTCTGCTCTCTTGTGAGTATGGGTTTATAATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAAGATCAATCCAACAAACAACCGGATAAAGCATAACAATCCCCTTTCTTCTAAGAAc                                                                                                                                                                                                                          <  1:356428/250‑1 (MQ=255)
          gCCAGCGATACTTTACCTGGTAAAGTAGCATTTTACCTCTCATAATTGTTTGTCCTTTTTTGTTACCTCTTAGGTTCTGCTCTCTTGTGAGTATGGGTTTATAATGATGTCATACGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAAGCATAACAATCCCCTTTCTTCTAAGAACTTCAt                                                                                                                                                                                                                     >  1:156144/1‑251 (MQ=255)
                                       aTTTTACCTCTCATAAGTGTTTGTCCTTCTTTGTTACCTCTTAGGTTCTGCTCTCTTGTGAGTATGGGTTTATAAGGATGTCATCCGGTACATTAACACGGTTGGAAACATCGGTTTTGCGGTGCACTTGCTGCATAAGCACCAGTTCATTGACATGGTAAATAGCAAGCCAACAAACATCCGGGTGAAGTTGAACAATCCCCTTTCTGCTAAGAACTTCATTCATcccaattgattatattcgggataa                                                                                                                                                                                         <  1:450362/250‑22 (MQ=255)
                                          ttACCTCTCATAATTGTTTTTCCTTCTTTGTTACCTCTTAGGTTCTGCTCTCTTGTGAGTATGGGTTTATAATGATGTCATCCGGTTCATTAACACGGTTGAAAACATCTATTTTGCGTAGCACTGGCAGCATAACCACCATGTCATTGACATGGTAAATATCATTCCAACAAACAACCGGATAAAGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGa                                                                                                                                                                                      <  1:503017/250‑1 (MQ=255)
                                                             gtCCTTCTTTGTTACCTCTTAGGTTCTGCTCTCTTGTGAGTATGGGTTTATAATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAAGCATAACAATCCCCTTTCTTCTAAGATCTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTagagn                                                                                                                                                                          <  1:98748/243‑2 (MQ=255)
                                                                        ttACCTCTTAGGTTCTGCTCTCTTGTGAGTATGGGTTTATAATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTGTTGCGGAGCACTGGCATCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAAGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTagagac                                                                                                                                                                         <  1:456771/233‑3 (MQ=255)
                                                                           ccTCTTAGGTTCTGCTCTCTTGTGAGTATGGGTTTATAATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAAGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCAAAGtaata                                                                                                                                                     >  1:321402/1‑250 (MQ=255)
                                                                              ctTAGGTTNTNNNCNNTTGTGAGTATGGGTTTATAATGATGTCATCCGGTACATTAACACNGTTGGAAACATCTGTNTTNCGGANCACNGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAAGCATAACAATCCCCTTTCTTCTAAGAACTTCa                                                                                                                                                                                                                      >  2:358173/1‑182 (MQ=255)
                                                                                                                   gatgTCATCCGGTACATTAACACGGTTGGAACCATCTGTTTTGCGGAGCACTGCCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAAGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATc                                                                                                                                                                                                                >  1:201478/1‑151 (MQ=255)
                                                                                                                             ggTACATTAACCCTGTTGGAAACATCTTTTGTGCGGAGCACTGGCAGCATAACCACCATGTCATTTACATGGTAAATATCATGCCAACAAACAACCGCATAAATCATAACAATCCCCTTTCTGCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATg                                                                                                                                                                                          <  1:341766/163‑1 (MQ=255)
                                                                                                                                           gTTGGAAACATCTGTTTGGCGGAGCACTGGCAGCATAACCACCATGTCATTGACATGGTAAATATTAATCCAACAAACAACCGGATAAAGCATAGCAATCCCCTTTCTTCTAAGAACGTCATTTATCCCGAATATAATCAATGGTGATGTGGACAAGAGTAGAGCCTGCCCGTCATCAAAGt                                                                                                                                                         <  1:321208/182‑1 (MQ=255)
                                                                                                                                                                                    ccaTTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAAGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGAAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTcc                                           >  1:382098/1‑251 (MQ=255)
                                                                                                                                                                                          cATTGACATGGTAAATATCAATCCAACAAACAACCGGATACAGCATAACCCTCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGcc                                                                                                                                                                         >  1:246164/1‑119 (MQ=255)
                                                                                                                                                                                                        atatCAGTNCNNNANNCAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCANTTATCCCGAATATAANCANTTGGNATGNTGACAAGAGTAGAGc                                                                                                                                                                          >  2:346057/1‑104 (MQ=255)
                                                                                                                                                                                                        anttNAATCCNACAANCAACCGGATAAAGCANAACAATCCCCTTTCTTCTANNANCTTCANTTATc                                                                                                                                                                                                                <  2:201478/63‑1 (MQ=255)
                                                                                                                                                                                                                     aaacaaCCGGATAAAGCATAACAATCCCCTTTCTTCTAAAAACTTCATTAGTCCCGAatat                                                                                                                                                                                                        >  1:450645/1‑61 (MQ=255)
                                                                                                                                                                                                                           ncGGATAAAGCATAACAATCCCCTTTCTTCTAAGAACTTGATTTATCCCGAATATAATCAATTGGGATGTTGACAag                                                                                                                                                                                  >  1:50852/2‑77 (MQ=255)
                                                                                                                                                                                                                              gATAAAGCATAACAATCCCCTTTCTTCTAAGAACTTCATTGATCCCGAATATAATCATTTGGGATGTTGACAAGAGAAAAGCGTGCCCGGCATAAAAGCAATAACCGTCAAAGGATGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCACAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGtat  >  1:229785/1‑250 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
GAGAGTGATTGCCAGCGATACTTTAACTGGTAAAGTAGTATTTTACCTCTCATAATTGTTTGTCCTTCTTTGTTACCTCTTAGGTTCTGCTCTCTTGTGAGTATGGGTTTATAATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTAT  >  NC_000913/293866‑294337

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: