Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 849,322 C→G 8.3% G366R (GGC→CGC)  pcs ← phosphatidylcholine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947849,3220CG 8.3% 142.4 / 4.6 48G366R (GGC→CGC) pcsphosphatidylcholine synthase
Reads supporting (aligned to +/- strand):  ref base C (23/21);  new base G (2/2);  total (25/23)
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.10e-01

AAACACACGCGCCAATGTGGGCATATTGTTCTCCTTGAAGAAATTATTTAGACAGGTTGTTTTCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGCATACTGGACAGATGCCATCAACGGCATGGCGCTGGTAGGTGTGCTCCTGACAAGATCGGCTTGCAATTTGAAAAATGC  >  NC_002947/849176‑849463
                                                                                                                                                  |                                                                                                                                             
aaaCACACGCGCCAATGTGGGCATATTGTTCTCCTTGAAGAAATTATTTAGACAGGTTGTTTTCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCaa                                                                                                                                             >  1:224734/1‑149 (MQ=255)
      acGCGCCAATGTGGGCATATTGTTCTCCTTGAAGAAATTATTTAGACAGGTTGTTTTCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAAc                                                                                                                                            >  7:48919/1‑144 (MQ=255)
      acGCGCCAATGTGGGCATATTGTTCTCCTTGAAGAAATTATTTAGACAGGTTGTTTTCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAAc                                                                                                                                            <  8:48919/144‑1 (MQ=255)
          gcCAATGTGGGCATATTGTTCTCCTTGAAGAACTTATCTAGACAGGTGCTTTTCCCTTCATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAAc                                                                                                                                   <  4:94538/149‑1 (MQ=255)
               tgtgGGCATATTGTTCTCCTTGAAGAAATTATTTAGACAGGTTGTTTTCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTa                                                                                                                              <  2:178962/149‑1 (MQ=255)
                      atatTGTTCTCCTTGAAGAAATTATTTAGACAGGTTGTTTTCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGcc                                                                                                                       >  8:22936/1‑149 (MQ=255)
                        atTGTTCTCCTTGAAGAAATTATTTAGACAGGTTGTTTTCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCGACATACTTAACCAGTATTGAGCCgg                                                                                                                     >  5:175865/1‑149 (MQ=255)
                                 cTTGAAGAAATTATTTAGACAGGTTGTTTTCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATATAAAGCAACGAAGATAATGCTCTTGCTGGACAGGCGAACATACTTAACCAGTATTGAGCCGGCGAAGATgg                                                                                                            >  7:86936/1‑149 (MQ=255)
                                 cTTGAAGAAATTATTTAGACAGGTTGTTTTCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCGAACATACTTAACCAGTATTGAGCCGGCGAAGATgg                                                                                                            >  7:91487/1‑149 (MQ=255)
                                 cTTGAAGAAATTATTTAGACAGGTTGTTTTCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGc                                                                                                                    >  4:60251/1‑141 (MQ=255)
                                 cTTGAAGAAATTATTTAGACAGGTTGTTTTCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGc                                                                                                                    <  3:60251/141‑1 (MQ=255)
                                 cTTGAAGAAATTATTTAGACAGGTTGTTTTCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATgg                                                                                                            >  3:278384/1‑149 (MQ=255)
                                 cTTGAAGAAATTATTTAGACAGGTTGTTTTCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATgg                                                                                                            >  5:58156/1‑149 (MQ=255)
                                              tttAGACAGGTTGTTTTCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGcc                                                                                               >  7:283096/1‑149 (MQ=255)
                                                tAGACAGGTTGTTTTCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCt                                                                                              >  1:302929/1‑148 (MQ=255)
                                                  gACAGGTTGTTTTCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAg                                                                                           >  8:20933/1‑149 (MQ=255)
                                                   aCAGGTTGTTTTCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGt                                                                                          >  4:280452/1‑149 (MQ=255)
                                                    cAGGTTGTTTTCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTa                                                                                         <  8:30777/149‑1 (MQ=255)
                                                    cAGGTTGTTTTCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTa                                                                                         <  6:58156/149‑1 (MQ=255)
                                                       gttgttTTCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGAc                                                                                      >  5:93154/1‑149 (MQ=255)
                                                           ttttCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCGACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGAccccc                                                                                  <  6:175865/149‑1 (MQ=255)
                                                           ttttCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGAccccc                                                                                  <  2:224734/149‑1 (MQ=255)
                                                              tCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGc                                                                                <  3:280452/148‑1 (MQ=255)
                                                                ccTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGCATa                                                                             <  7:20933/149‑1 (MQ=255)
                                                                 cTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCa                                                                                                   <  5:260242/126‑1 (MQ=255)
                                                                 cTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCa                                                                                                   >  6:260242/1‑126 (MQ=255)
                                                                  ttGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGCATACt                                                                           >  5:274433/1‑149 (MQ=255)
                                                                   tGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAGGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGCATACTg                                                                          >  6:194721/1‑149 (MQ=255)
                                                                                tCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGCATACTGGACAGATGCCAt                                                              >  1:212516/1‑148 (MQ=255)
                                                                                  gTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGCATACTGGACAGATGCCATCaa                                                           >  7:228756/1‑149 (MQ=255)
                                                                                   ttCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGCATACTGGACAGATGCCATCAAc                                                          >  1:289044/1‑149 (MQ=255)
                                                                                          ggCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGCATACTGGACAGATGCCATCAACGGCATgg                                                   <  2:302929/149‑1 (MQ=255)
                                                                                                 gATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGTCGAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGCATACTGGACAGATGCCATCAACGGCATGGCGCTGGt                                            <  8:86936/149‑1 (MQ=255)
                                                                                                 gATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCGAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGCATACTGGACAGATGCCATCAACGGCATGGCGCTGGt                                            <  8:91487/149‑1 (MQ=255)
                                                                                                  aTGCGACGAAATAGAAAGCAACGAATATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGCATACTGGACAGATGCCATCAACGGCATGGCGCTGGTa                                           <  6:138273/149‑1 (MQ=255)
                                                                                                      gacgaAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGCATACTGGACAGATGCCATCAACGGCATGGCGCTGGTAGgtg                                       >  1:143198/1‑149 (MQ=255)
                                                                                                        cgaAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGCATACTGGACAGATGCCATCAACGGCATGGCGCTGGTAGGtgtg                                     <  4:278384/149‑1 (MQ=255)
                                                                                                           aaTAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGCATACTGGACAGATGCCATCAACGGCATGGCGCTGGTAGGTGTGCTc                                  >  4:227145/1‑149 (MQ=255)
                                                                                                                  aGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGCATACTGGACAGATGCCATCAACGGCATGGCGCTGGTAGGTGTGCTCCTGAc                             >  6:63918/1‑147 (MQ=255)
                                                                                                                  aGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGCATACTGGACAGATGCCATCAACGGCATGGCGCTGGTAGGTGTGCTCCTGACaa                           <  7:334/149‑1 (MQ=255)
                                                                                                                     aaCGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGAccccc                                                                                  <  7:51903/91‑1 (MQ=255)
                                                                                                                     aaCGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGAccccc                                                                                  >  8:51903/1‑91 (MQ=255)
                                                                                                                     aaCGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGCATACTGGACAGATGCCATCAACGGCATGGCGCTGGTAGGTGTGCTCCTGACAAgtt                        <  2:5800/149‑3 (MQ=255)
                                                                                                                     aaCGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGAAGGTGATTATCATGCCTAAGTAGACCCCCGCATACTGGACAGATGCCATCAACGGCATGGCGCTGGTAGGTGTGCTCCTGACAAGAt                        <  1:15310/149‑1 (MQ=255)
                                                                                                                          aGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGCATACTGGACAGATGCCATCAACGGCATGGCGCTGGTAGGTGTGCTCCTGACAAGATCGGCt                   <  3:227145/149‑1 (MQ=255)
                                                                                                                                           aaCAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGCATACTGGACAGATGCCATCAACGGCATGGCGCTGGTAGGTGTGCTCCTGACAAGATCGGCTTGCAATTTGAAAAAt    <  5:63918/147‑1 (MQ=255)
                                                                                                                                           aaCAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGCATACTGGACAGATGCCATCAACGGCATGGCGCTGGTAGGTGTGCTCCTGACAAGATCGGCTTGCAATTTGAAAAATGc  >  4:287125/1‑149 (MQ=255)
                                                                                                                                            aCAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGCATACTGGACAGATGCCATCAACGGCATGGCGCTGGTAGGTGTGCTCCTGACAAGATCGGCTTGCAATTTGAAAAATGc  <  2:289044/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
AAACACACGCGCCAATGTGGGCATATTGTTCTCCTTGAAGAAATTATTTAGACAGGTTGTTTTCCCTTGATCGCGAACTGTCGTTCGCAGGGCGTTGGATGCGACGAAATAGAAAGCAACGAAGATAATGCTCGTGCTGAACAGGCCAACATACTTAACCAGTATTGAGCCGGCGAAGATGGTGATTATCATGCCTAAGTAGACCCCCGCATACTGGACAGATGCCATCAACGGCATGGCGCTGGTAGGTGTGCTCCTGACAAGATCGGCTTGCAATTTGAAAAATGC  >  NC_002947/849176‑849463

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