Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,205,303 G→A 15.4% L273F (CTT→TTT)  ddlA ← D‑alanine‑‑D‑alanine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,205,3030GA15.4% 60.7 / 5.1 26L273F (CTT→TTT) ddlAD‑alanine‑‑D‑alanine ligase
Reads supporting (aligned to +/- strand):  ref base G (12/10);  new base A (2/2);  total (14/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 = 7.99e-01

ACATACTGAAAGCCGTAAATCCAGGCATTGCATTTGTTTCATTAATATATATTTGGTTGTCTTCTGTTACAAAGAAATCAGCACGGACTAAACCAGAACAATCTGTCGCTTTGAATGCCTCTAATGCCATATTTCTAAGCGTTAATTGAACATCTTCGTCTAAGTCAGCTGGAATTTGTAATTGAACCTTACCATCTTTATATTTTGATTTGTAATCGTAAAACGCGACATCTTTTACGACTTCACCTGGCCATGTCGCTTCAGGATAGTCATTTCCTA  >  CP000730/2205165‑2205443
                                                                                                                                          |                                                                                                                                            
ctatacTGAAAACCGTAAATCCAGGCATTGCATTTGTTTCATTAATATATATTTGGTTGTCTTCTGTTACAAAGAAATCAGCACGGACTAAACCAGAACAATCTGTCGCTTTGAATGCCTCTAATGCCATATTTCTAAGCg                                                                                                                                            <  2:253315/139‑1 (MQ=255)
                           ttgcattTGTTTCATTAATATATATTTGGTTGTCTTCTGTTACAAAGAAATCAGCACGGACTAAACCAGAACAATCTGTCGCTTTGAATGCCTCTCATGCCATCTTTCTAAGCGTTAATTGAACATCTTCGTCTAAGTCAg                                                                                                                 >  1:11918/1‑141 (MQ=255)
                           ttgcattTGTTTCATTAATATATATTTGGTTGTCTTCTGTTACAAAGAAATCAGCACGGACTAAACCAGAACAATCTGTCGCTTTGAATGCCTCTAATGCCATATTTCTAAGCGTTAATTGAACGTCTTCGTCTAATTCAg                                                                                                                 >  1:13172/1‑141 (MQ=255)
                              cattTGTTTCATTAATATATATTTGGTTGTCTTCTGTTACAAAGAAATCAGCACGGACTAAACCAGAACAATCTGTCGCTTTGAATGCCTCTAATGCCATATTTCTAAGCGTTAATTGAACATCTTCGTCTAAGTCAGCTg                                                                                                              >  1:175577/1‑141 (MQ=255)
                                                        ttGTCTTCTGTTACAAAGAAATCAGCACGGACTAAACCAGAACAATCTGTCGCTTTGAATGCCTCTAATGCCATATTTCTAAGCGTTAATTGAACATCTTCGTCTAAGTCAGCTGGAATTTGTAATTGAACCTTACCATCt                                                                                    <  2:272258/141‑1 (MQ=255)
                                                         tGTCTTCTGTTACAAAGAAATCAGCACGGACTAAACCAGAACAATCTGTCGCTTTGAATGCCTCTAATGCCATATTTCTAAGCGTTAATTGAACATCTTCGTCTAAGTCAGCTGGAATTTGTAATTGAACCTTACCATCtt                                                                                   <  2:218137/141‑1 (MQ=255)
                                                          gTCTTCTGTTACAAAGAAATCAGCACGGACTAAACCAGAACAATCTGTCGCTTTGAATGCCTCTAATGCCATATTTCTAAGCGTTAATTGAACATCTTCGTCTAAGTCAGCTGGAATTTGTAATTGAACCTTACCATCttt                                                                                  >  2:285544/1‑141 (MQ=255)
                                                           tcttctGTTACAAAGAAATCAGCACGGACTAAACCAGAACAATCTGTCGCTTTGAATGCCTCTAATGCCATATTTCTAAGCGTTAATTGAACATCTTCGTCTAAGTCAGCTGGAATTTGTAATTGAACCTTACCATCTTta                                                                                 >  2:53408/1‑141 (MQ=255)
                                                              tctGTTACAAAGAAATCAGCACGGACTAAACCAGAACAATCTGTCGCTTTGAATGCCTCTAATGCCATATTTCTAAGCGTTAATTGAACATCTTCGTCTAAGTCAGCTGGAATTTGTAATTGAACCTTACCATCTTTatat                                                                              >  1:258798/1‑141 (MQ=255)
                                                               ctGTTACAAAGAAATCAGCACGGACTAAACCAGAACAATCTGTCGCTTTGAATGCCTCTAATGCCATATTTCTAAGCGTTAATTGAACATCTTCGTCTAAGTCAGCTGGAATTTGTAATTGAACCTTACCATCTTTATAtt                                                                             <  2:49432/141‑1 (MQ=255)
                                                               ctGTTACAAAGAAATCAGCACGGACTAAACCAGAACAATCTGTCGCTTTGAATGCCTCTAATGCCATATTGCTAAGCGTTAATTCCACCACGTCGGCTAAGTCAGCTGGAAGTTGTAATGGAACCTTACCATCATtataat                                                                             >  1:268449/1‑139 (MQ=255)
                                                                   tACAAAGAAATCAGCACGGACTAAACCAGAACAATCTGTCGCTTTGAATGCCTCTAATGCCATATTTCTAAACGTTAATTGAACATCTTCGTCTAAGTCAGCTGGAATTTGTAATTGAACCTTACCATCTTTATATTTTGa                                                                         >  2:91091/1‑141 (MQ=255)
                                                                   tACAAAGAAATCAGCACGGACTAAACCAGAACAATCTGTCGCTTTGAATGCCTCTAATGCCATATTTCTAAACGTTAATTGAACATCTTCGTCTAAGTCAGCTGGAATTTGTAATTGAACCTTACCATCTTTATATTTTGa                                                                         >  2:91052/1‑141 (MQ=255)
                                                                      aaaGAAATCAGCACGGACTAAACCAGAACAATCTGTCGCTTTGAATGCCTCTAATGCCATATTTCTAAGCGTTAATTGAACATCTTCGTCTAAGTCAGCTGGCATTTGTAATTGAACCTTACCATCTTTATATTTTGAttt                                                                      >  1:166588/1‑141 (MQ=255)
                                                                            aTCAGCACGGACTAAACCAGAACAATCTGGCGCTTTGAATGCCTCTAATGCCATATTTCTAAGCGTTAATTGAACATCTTCGTCTAAGTCAGCTGGAATTTGAAATTGAACCTTACCATCTTTATATTTTGATTTGTACTc                                                                <  2:74706/141‑1 (MQ=255)
                                                                                  aCGGACTAAACCAGAACAATCTGTCGCTTTGAATGCCTCTAATGCCATATTTCTAAGCGTTAATTGAACATCTTCGTCTAAGTCAGCTGGAATTTGTAATTGAACCTTACCATCTTTATATTTTGATTTGTAATCGTaaaa                                                          >  1:61251/1‑141 (MQ=255)
                                                                                  aCGGACTAAACCAGAACAATCTGTCGCTTTGAATGCCTCTAATGCCATATTTCTAAGCGTTAATTGAACATCTTCGTCTAAGTCAGCTGGAATTTGTAATTGAACCTTACCATCTTTATATTTTGATTTGTAATCGTaaaa                                                          <  1:22933/141‑1 (MQ=255)
                                                                                          acccAGAACAATCTGTCGCTTTGAATGCCGCTAATGCCATATTTCTAAGCGTTAATTGAACATCTTCGTCTAAGTCAGCTGGAATTTGTAATTGAACCTTACCATCTTTATATTTTGATTTGTAATCGTAAAACGCGACAt                                                  <  2:128779/139‑1 (MQ=255)
                                                                                                             tttGAATGCCTCTAATGCCATATTTCTAAGCGTTAATTGAACATCTTCGTCTAAGTCAGCTGGAATTTGTAATTGAACCTTACCATCTTTATATTTTGATTTGTAATCGTAAAACGCGACATCTTTTACGACTTCACCTgg                               <  1:104477/141‑1 (MQ=255)
                                                                                                                  aTGCCTCTAATGCCATATTTCTAAGCGTTAATTGAACATCTTCGTCTAAGTCAGCTGGGATTTGTAATTGAACCTTTCCATCTTTATATTTTGATTTGTCATCGTAAAACGCGACATCTTTTACGGCTTCAACTGGcgcag                          >  1:71565/1‑137 (MQ=255)
                                                                                                                             gCCATATTTCTAAGCGTTAATTGAACATCTTCGTCTAAGTCAGCTGGAATTTGTAATTGAACCTTACCATCTTTATATTTTGATTTGTAATCGTAAAACGCGACATCTTTTACGACTTCACCTGGCCATGTCGCTTCGGGa               >  2:21102/1‑141 (MQ=255)
                                                                                                                              ccATATTTCTAAACGTTAATTGAACATCTTCGTCTAAGTCAGCTGGAATTTGTAATTGAACCTTACCATCTTTATATTTTGATTTGTAATCGTAAAACGCGACATCTTTTACGACTTCACCTGGCCATGTCGCTTCAGGAt              <  1:91052/141‑1 (MQ=255)
                                                                                                                              ccATATTTCTAAACGTTAATTGAACATCTTCGTCTAAGTCAGCTGGAATTTGTAATTGAACCTTACCATCTTTATATTTTGATTTGTAATCGTAAAACGCGACATCTTTTACGACTTCACCTGGCCATGTCGCTTCAGGAt              <  1:91091/141‑1 (MQ=255)
                                                                                                                                    ttCTAAGCGTTAATTGAACATCTTCGTCTAAGTCAGCTGGAATTTGTAATTGAACCTTACCATCTTTATATTTTGATTTGTAATCGTAAAACGCGACATCTTTTACGACTTCACCTGGCCATGTCGCTTCAGGATAGTCAt        <  2:175577/141‑1 (MQ=255)
                                                                                                                                       tAAGCGTTAATTGAACATCTTCGTCTAAGTCAGCTGGAATTTGTAATTGAACCTTACCATCTTTATATTTTGATTTGTAATCGTACAACGCGACATCTTTTTCTACTGCACCTGGCCATGTCGCTTCAGGATAGTCATTTc     >  1:194225/1‑141 (MQ=255)
                                                                                                                                         aGCGTTAATTGAACATCTTCGTCTAAGTCAGCTGGAATTTGTAATTGAACCTTACCATCTTTATATTTTGATTTGTAATCGTAAAACGCGACATCTTTTACGACTTCACCTGGCCATGTCGCTTCAGGATAGTCATTTCCt   <  1:285544/141‑1 (MQ=255)
                                                                                                                                          gCGTTAATTGCACATCTTCGTCTAATTCAGTTGGAATTTGTAATTGAACCTTACCATCTTTATATTTTGATTTGTAATCGTAAAACGCGACATCTTTTACGACTTCACCTGGCCATGTCGCTTCAGGATAGTCATTTCCTa  <  1:166199/141‑1 (MQ=255)
                                                                                                                                          gCGTTAATTGAACATCTTCGTCTAAGTCAGTTGGAATTTGTAATTGACCCTTACCATCTTTATATTTTGATTTGTAATCGTAAAACGCGACATCTTTTACGACTTCACCTGGCCATGTCGCTTCAGGATAGTCATTTCCTa  <  1:166214/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                            
ACATACTGAAAGCCGTAAATCCAGGCATTGCATTTGTTTCATTAATATATATTTGGTTGTCTTCTGTTACAAAGAAATCAGCACGGACTAAACCAGAACAATCTGTCGCTTTGAATGCCTCTAATGCCATATTTCTAAGCGTTAATTGAACATCTTCGTCTAAGTCAGCTGGAATTTGTAATTGAACCTTACCATCTTTATATTTTGATTTGTAATCGTAAAACGCGACATCTTTTACGACTTCACCTGGCCATGTCGCTTCAGGATAGTCATTTCCTA  >  CP000730/2205165‑2205443

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