Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 196,002 C→T S109F (TCC→TTC)  cdsA → CDP‑diglyceride synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913196,0020CT96.9% 107.5 / ‑3.9 32S109F (TCC→TTC) cdsACDP‑diglyceride synthetase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (13/18);  total (14/18)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

CGCTGGTTGAAATCTCACTTTGGGCTTCGCTGGGTTGGTGGATTGTCGCGCTATTGCTGGTGCTGTTTTACCCAGGTTCCGCAGCAATCTGGCGTAACTCTAAAACATTGCGCCTTATTTTTGGCGTGCTAACCATTGTTCCCTTCTTCTGGGGC  >  NC_000913/195924‑196078
                                                                              |                                                                            
cGCTGGTTGAAATCTCACTTTGGGCTTCGCTGGGTTGGTGGATTGTCGCGCTATTGCTGGTGCTGTTTTACCCAGGTTTCg                                                                            <  2:45159/81‑1 (MQ=255)
 gCTGGTTGAAATCTCACTTTGGGCTTCGCTGGGTTGGTGGATTGTCGCGCTATTGCTGGTGCTGTTTTACCCAGGTTTCg                                                                            >  2:45202/1‑80 (MQ=255)
  cTGGTTGAAATCTCACTTTGGGCTTCGCTGGGTTGGTGGATTGTCGCGCTATTGCTGGTGCTGTTTTACCCAGGTTCCg                                                                            >  1:45195/1‑79 (MQ=255)
     gTTGAAATCTCACTTTGGGCTTGGCTGGGTTGGTGGATTGTCGCGCTATTGCTGGTGCTGTTTTACCCAGGTTTCgcagca                                                                       <  1:45160/81‑1 (MQ=255)
       tGAAATCTCACTTTGGGCTTCGCTGGGTTGGTGGATTGTCGCGCTATTGCTGGTGCTGTTTTACCCAGGTTTCGCAGCAAt                                                                     >  2:45217/1‑81 (MQ=255)
          aaTCTCACTTTGGGCTTCGCTGGGTTGGTGGATTGTCGCGCTATTGCTGGTGCTGTTTTACCCAGGTTTCGCAGCAATCTg                                                                  >  1:45204/1‑81 (MQ=255)
                aCTTTGGGCTTCGCTGGGTTGGTGGATTGTCGCGCTATTGCTGGTGCTGTTTTACCCAGGTTTCGCAGCAATCTGGc                                                                <  2:45161/77‑1 (MQ=255)
                aCTTTGGGCTTCGCTGGGTTGGTGGATTGTCGCGCTATTGCTGGTGCTGTTTTACCCAGGTTTCGCAGCAATCTGGCGTaa                                                            >  1:45212/1‑81 (MQ=255)
                   ttGGGCTTCGCTGGGTTGGTGGATTGTCGCGCTATTGCTGGTGCTGTTTTACCCAGGTTTCGCAGCAATCTGGCGTAActc                                                         >  2:45193/1‑81 (MQ=255)
                      ggCTTCGCTGGGTTGGTGGATTGTCGCGCTATTGCTGGTGCTGTTTTACCCAGGTTTCGCAGCAATCTGGCGTAACTCTa                                                       <  2:45162/80‑1 (MQ=255)
                      ggCTTCGCTGGGTTGGTGGATTGTCGCGCTATTGCTGGTGCTGTTTTACCCAGGTTTCGCAGCAATCTGGCGTAACTCTa                                                       <  2:45163/80‑1 (MQ=255)
                      ggCTTCGCTGGGTTGGTGGATTGTCGCGCTATTGCTGGTGCTGTTTTACCCAGGTTTCGCAGCAATCTGGCGTAACTCTa                                                       <  1:45164/80‑1 (MQ=255)
                            gCTGGGTTGGTGGATTGTCGCGCTATTGCTGGTGCTGTTTTACCCAGGTCTCGCAGCAATCTGGCGTAACTCTAAAACAtt                                                >  2:45213/1‑81 (MQ=255)
                                  ttGGTGGATTGTCGCGCTATTGCTGGTGCTGTTTTACCCAGGTTTCGCAGCAATCTGGCGTAACTCTAAAACATTGCGCCt                                          <  1:45165/81‑1 (MQ=255)
                                       ggATTGTCGCGCTATTGCTGGTGCTGTTTTACCCAGGTTTCGCAGCAATCTGGCGTAACTCTAAAACATTGCGCCTTAttt                                     >  2:45207/1‑81 (MQ=255)
                                       ggATTGTCGCGCTATTGCTGGTGCTGTTTTACCCAGGTTTCGCAGCAATCTGGCGTAACTCTAAAACATTGCGCCTTAttt                                     >  2:45211/1‑81 (MQ=255)
                                        gATTGTCGCGCTATTGCTGGTGCTGTTTTACCCAGGTTTCGCAGCAATCTGGCGTAACTCTAAAACATTGCGCCTTAtttt                                    <  1:45166/81‑1 (MQ=255)
                                        gATTGTCGCGCTATTGCTGGTGCTGTTTTACCCAGGTTTCGCAGCAATCTGGCGTAACTCTAAAACATTGCGCCTTAttt                                     <  2:45167/80‑1 (MQ=255)
                                             tCGCGCTATTGCTGGTGCTGTTTTACCCAGGTTTCGCAGCAATCTGGCGTAACTCTAAAACATTGCGCCTTATTTTTGGc                                <  1:45168/80‑1 (MQ=255)
                                               gcgcTATTGCTGGTGCTGTTTTACCCAGGTTTCGCAGCAATCTGGCGTAACTCTAAAACATTGCGCCTTATTTTTGGCGTg                             >  2:45229/1‑81 (MQ=255)
                                                 gcTATTGCTGGTGCTGTTTTACCCAGGTTTCGCAGCAATCTGGCGTAACTCTAAAACATTGCGCCTTATTTTTGGCGTGc                            >  2:45220/1‑80 (MQ=255)
                                                     ttGCTGGTGCTGTTTTACCCAGGTTTCGCAGCAATCTGGCGTAACTCTAAAACATTGCGCCTTATTTTTGGCGTGCTAAcc                       <  2:45170/81‑1 (MQ=255)
                                                     ttGCTGGTGCTGTTTTACCCAGGTTTCGCAGCAATCTGGCGTAACTCTAAAACATTGCGCCTTATTTTTGGCGTGCTAAcc                       <  1:45169/81‑1 (MQ=255)
                                                     ttGCTGGTGCTGTTTTACCCAGGTTTCGCAGCAATCTGGCGTAACTCTAAAACATTGCGCCTTATTTTTGGCGTGCTAAc                        >  1:45216/1‑80 (MQ=255)
                                                           gTGCTGTTTTACCCAGGTTTCGCAGCAATCTGGCGTAACTCTAAAACATTGCGCCTTATTTTTGGCGTGCTAACCATTGtt                 >  2:45219/1‑81 (MQ=255)
                                                            tGCTGTTTTACCCAGGTTTCGCAGCAATCTGGCGTAACTCTAAAACATTGCGCCTTATTTTTGGCGTGCTAACCATTGTTc                <  2:45171/81‑1 (MQ=255)
                                                              cTGTTTTACCCAGGTTTCGCAGCAATCTGGCGTAACTCTAAAACATTGCGCCTTATTTTTGGCGTGCTAACCATTGTTcc               <  1:45172/80‑1 (MQ=255)
                                                                 ttttACCCAGGTTTCGCAGCAATCTGGCGTAACTCTAAAACATTGCGCCTTATTTTTGGCGTGCTAACCATTGTTCCct             <  2:45173/79‑1 (MQ=255)
                                                                        cAGGTTTCGCAGCAATCTGGCGTAACTCTAAAACATTGCGCCTTATTTTTGGCGTGCTAACCATTGTTCCCTTCTTCTggg    <  2:45174/81‑1 (MQ=255)
                                                                         aGGTTTCGCAGCAATCTGGCGTAACTCTAAAACATTGCGCCTTATTTTTGGCGTGCTAACCATTGTTCCCTTCTTCTgggg   <  2:45175/81‑1 (MQ=255)
                                                                          ggTTTCGCAGCAATCTGGCGTAACTCTAAAACATTGCGCCTTATTTTTGGCGTGCTAACCATTGTTCCCTTCTTCTgggg   >  2:45209/1‑80 (MQ=255)
                                                                          ggTTTCGCAGCAATCTGGCGTAACTCTAAAACATTGCGCCTTATTTTTGGCGTGCTAACCATTGTTCCCTTCTTCTGGGGc  <  1:45176/81‑1 (MQ=255)
                                                                              |                                                                            
CGCTGGTTGAAATCTCACTTTGGGCTTCGCTGGGTTGGTGGATTGTCGCGCTATTGCTGGTGCTGTTTTACCCAGGTTCCGCAGCAATCTGGCGTAACTCTAAAACATTGCGCCTTATTTTTGGCGTGCTAACCATTGTTCCCTTCTTCTGGGGC  >  NC_000913/195924‑196078

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: