Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,697,952 C→G 8.3% R200P (CGA→CCA)  pgpC ← phosphatidylglycerophosphatase C, membrane bound

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,697,9520CG 8.3% 159.5 / 5.4 48R200P (CGA→CCA) pgpCphosphatidylglycerophosphatase C, membrane bound
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 = 7.85e-01

AGCGCGTGACGCATCCAGTATTCGTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTGCGCTCCAGTTGTGCGACCTTTTCATGTCCCAGACAACGCATCGTCAATACCCAACCAC  >  NC_000913/2697811‑2698099
                                                                                                                                             |                                                                                                                                                   
aGCGCGTGACGCATCCAGTATTCGTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTGGCCAACg                                                                                                                                              >  4:438406/1‑149 (MQ=255)
aGCGCGTGACGCATCCAGTATTCGTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTGGCCAACg                                                                                                                                              >  4:315727/1‑149 (MQ=255)
  cgcgTGACGCATCCAGTATTCGTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGAt                                                                                                                                            <  1:240411/149‑1 (MQ=255)
         cGCATCCAGTATTCGTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACa                                                                                                                                     <  1:68342/149‑1 (MQ=255)
            aTCCAGTATTCGTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGaa                                                                                                                                  >  1:422350/1‑149 (MQ=255)
             tCCAGTATTCGTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGaaa                                                                                                                                 <  1:180641/149‑1 (MQ=255)
                       gTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATaaa                                                                                                                             <  1:326206/143‑1 (MQ=255)
                       gTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATaaa                                                                                                                             >  2:326206/1‑143 (MQ=255)
                            tAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAagcagc                                                                                                                        <  1:416848/143‑1 (MQ=255)
                            tAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAagcagc                                                                                                                        >  2:416848/1‑143 (MQ=255)
                               aTTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTcc                                                                                                               >  4:256548/1‑149 (MQ=255)
                               aTTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTcc                                                                                                               >  3:370491/1‑149 (MQ=255)
                                   gACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGAt                                                                                                           <  1:444937/149‑1 (MQ=255)
                                      ttCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTAc                                                                                                        >  2:76988/1‑149 (MQ=255)
                                      ttCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTAc                                                                                                        >  1:69505/1‑149 (MQ=255)
                                             aaaaaGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCt                                                                                                 >  4:382517/1‑149 (MQ=255)
                                              aaaaGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATCACGCAGCGGATTGTCCTGATTACTGTCGCta                                                                                                >  2:454963/1‑149 (MQ=255)
                                                aaGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGa                                                                                                            <  3:414624/135‑1 (MQ=255)
                                                aaGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTGCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGAAGCGGATTGTCCTGa                                                                                                            >  4:414624/1‑135 (MQ=255)
                                                 aGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAg                                                                                             >  1:329403/1‑149 (MQ=255)
                                                 aGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACCGTCGCTATAg                                                                                             <  4:360584/149‑1 (MQ=255)
                                                         ccGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATaaa                                                                                                                             >  3:454266/1‑109 (MQ=255)
                                                         ccGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATaaa                                                                                                                             <  4:454266/109‑1 (MQ=255)
                                                                                 ataCACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGAtttt                                                             <  1:454963/149‑1 (MQ=255)
                                                                                    cacGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCAc                                                                                         >  3:355193/1‑118 (MQ=255)
                                                                                    cacGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCAc                                                                                         <  4:355193/118‑1 (MQ=255)
                                                                                     acGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTgcgc                                                         >  4:150297/1‑149 (MQ=255)
                                                                                       ggACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTGCGCTc                                                       <  2:69505/149‑1 (MQ=255)
                                                                                       ggACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTGCGCTc                                                       >  1:316454/1‑149 (MQ=255)
                                                                                         aCGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTGCGCTCCa                                                     <  2:422350/149‑1 (MQ=255)
                                                                                           gCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTGCGCTCCAGt                                                   >  4:401213/1‑149 (MQ=255)
                                                                                              aTGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTGCGCTCCAGTtgt                                                <  1:261426/149‑1 (MQ=255)
                                                                                              aTGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTGCGCTCCAGTtgt                                                <  2:329403/149‑1 (MQ=255)
                                                                                                       tctATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTGCGCTCCAGTTGTGCGACCttt                                       >  4:64219/1‑149 (MQ=255)
                                                                                                          aTTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGc                                                                    >  3:324944/1‑117 (MQ=255)
                                                                                                          aTTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGc                                                                    <  4:324944/117‑1 (MQ=255)
                                                                                                           ttCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTGCGCTCCAGTTGTGCGACCTTTTCAt                                   >  1:294402/1‑149 (MQ=255)
                                                                                                            tCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTGCGCTCCAGTTGTGCGACCTTTTCATg                                  >  2:432689/1‑149 (MQ=255)
                                                                                                             ccAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTGCGCTCCAGTTGTGCGACCTTTTCATGt                                 >  2:273725/1‑149 (MQ=255)
                                                                                                                gTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTGCGCTCCAGTTGTGCGACCTTTTCATGTccc                              <  3:4486/149‑1 (MQ=255)
                                                                                                                       gAGTTCACCGCGCGGGGTTACTGGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTGCGCTCCAGTTGTGCGACCTTTTCATGTCCCAGACAAc                       <  3:438406/149‑1 (MQ=255)
                                                                                                                       gAGTTCACCGCGCGGGGTTACTGGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTGCGCTCCAGTTGTGCGACCTTTTCATGTCCCAGACAAc                       <  3:315727/149‑1 (MQ=255)
                                                                                                                          ttCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTGCGCTCCAGTTGTGCGACCTTTTCATGTCCCAGACAACGCa                    <  1:273725/149‑1 (MQ=255)
                                                                                                                             aCCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTGCGCTCCAGTTGTGCGACCTTTTCATGTCCCAGACAACGCATCg                 >  4:464730/1‑149 (MQ=255)
                                                                                                                               cgcgcgGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTGCGCTCCAGTTGTGCGACCTTTTCATGTCCCAGACAACGCATCGTc               <  3:246387/149‑1 (MQ=255)
                                                                                                                                 cgcgGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTGCGCTCCAGTTGTGCGACCTTTTCATGTCCCAGACAACGCATCGTCaa             <  2:453189/149‑1 (MQ=255)
                                                                                                                                    ggggTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTGCGCTCCAGTTGTGCGACCTTTTCATGTCCCAGACAACGCATCGTCAATAc          <  1:76988/149‑1 (MQ=255)
                                                                                                                                            tCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTGCGCTCCAGTTGTGCGACCTTTTCATGTCCCAGACAACGCATCGTCAATACCCAaccac  >  4:113634/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
AGCGCGTGACGCATCCAGTATTCGTGGCTAAATTCGACTTCAGACAAAAAGAAAACTCCGGTTATAAAAGCGCGGCGCATTATACACGGACGCTATGCTTTACTCTATTCCAGTTGCTGGAGTTCACCGCGCGGGGTTACTCGCCAACGATGCTGACAGAAATAAAGCAGCGGATTGTCCTGATTACTGTCGCTATAGCCACTGTACAGCCGCAGCGGAGTGCCGATTTTGCGCTCCAGTTGTGCGACCTTTTCATGTCCCAGACAACGCATCGTCAATACCCAACCAC  >  NC_000913/2697811‑2698099

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