Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,789,596 C→T 9.2% A105V (GCC→GTC)  cdsA → CTP:phosphatidate cytidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,789,5960CT 9.2% 149.9 / 3.7 54A105V (GCC→GTC) cdsACTP:phosphatidate cytidyltransferase
Reads supporting (aligned to +/- strand):  ref base C (23/26);  new base T (3/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCCTATGCAGCGGTGGTCGCGGGGGCGCTGATGCTGCTGCACATCCTTCCGGAACTGGCGCCCTGGGTGCTGGGCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTGGGCAACTGGCTGATCCTGTCGGTCATGGTGCTGGTGTGGGCCGCCGACATCGGTGCGTACTTC  >  NC_002947/1789451‑1789744
                                                                                                                                                  |                                                                                                                                                    
gCCTATGCAGCGGTGGTCGCGGGGGCGCTGATGCTGCTGCACATCCTTCCGGAACTGGCGCCCTGGGTGCTGGGCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAg                                                                                                                                                   >  3:3013/1‑149 (MQ=255)
     tGCAGCGGTGGTCGCGGGGGCGCTGATGCTGCTGCACATCCTTCCGGAACTGGCGCCCTGGGTGCTGGGCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCgg                                                                                                                                              >  3:34462/1‑149 (MQ=255)
          cGGTGGTCGCGGGGGCGCTGATGCTGCTGCACATCCTTCCGGAACTGGCGCCCTGGGTGCTGGGCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTAGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGt                                                                                                                                         >  8:271673/1‑149 (MQ=255)
               gtCGCGGGGGCGCTGATGCTGCTGCACATCCTTCCGGAACTGGCGCCCTGGGTGCTGGGCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGgtt                                                                                                                                    >  1:139127/1‑149 (MQ=255)
                  gcgGGGGCGCTGATGCTGCTGCACATCCTTCCGGAACTGGCGCCCTGGGTGCTGGGCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGgttgtt                                                                                                                                 >  2:104640/1‑149 (MQ=255)
                  gcgGGGGCGCTGATGCTGCTGCACATCCTTCCGGAACTGGCGCCCTGGGTGCTGGGCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGgttgtt                                                                                                                                 <  8:132376/149‑1 (MQ=255)
                  gcgGGGGCGCTGATGCTGCTGCACATCCTTCCGGAACTGGCGCCCTGGGTGCTGGGCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGgttgt                                                                                                                                  <  2:203273/148‑1 (MQ=255)
                   cgGGGGCGCTGATGCTGCTGCACATCCTTCCGGAACTGGCGCCCTGGGTGCTGGGCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTAGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGttgttg                                                                                                                                <  7:271673/149‑1 (MQ=255)
                       ggCGCTGATGCTGCTGCACATCCTTCCGGAACTGGCGCCCTGGGTGCTGGGCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATc                                                                                                                             <  1:71005/149‑1 (MQ=255)
                        gcgcTGATGCTGCTGCACATCCTTCCGGAACTGGCGCCCTGGGTGCTGGGCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCgg                                                                                                                           >  2:35240/1‑149 (MQ=255)
                         cgcTGATGCTGCTGCACATCCTTCCGGAACTGGCGCCCTGGGTGCTGGGCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGc                                                                                                                          <  4:34462/149‑1 (MQ=255)
                                 ctgctgcACATCCTTCCGGAACTGGCGCCCTGGGTGCTGGGCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGgt                                                                                                                  <  2:191766/149‑1 (MQ=255)
                                  tgctgcACATCCTTCCGGAACTGGCGCCCTGGGTGCTGGGCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGgtg                                                                                                                 >  7:248854/1‑149 (MQ=255)
                                      gcACATCCTTCCGGAACTGGCGCCCTGGGTGCTGGGCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTtgc                                                                                                             >  3:74221/1‑149 (MQ=255)
                                       cacaTCCTTCCGGAACTGGCGCCCTGGGTGCTGGGCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTtgct                                                                                                            <  1:104640/149‑1 (MQ=255)
                                             cTTCCGGAACTGGCGCCCTGGGTGCTGGGCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGCCCTGTGGGTCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGc                                                                                                      >  5:221453/1‑149 (MQ=255)
                                              ttCCGGAACTGGCGCCCTGGCTGCTGGGCGCTGCCGTGAACTGGAGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCACCAGCGACCTGTGCGCCAGTGCGGCCTGACGGGTCTTGATCGGCATGCTGGTGTTGCTGCCGGc                                                                                                      >  6:10143/1‑148 (MQ=255)
                                                ccGGAACTGGCGCCCTGGGTGCTGGGCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGcc                                                                                                     <  6:8443/147‑1 (MQ=255)
                                                ccGGAACTGGCGCCCTGGGTGCTGGGCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTg                                                                                                   >  5:323861/1‑149 (MQ=255)
                                                 cGGAACTGGCGCCCTGGGTGCTGGGCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTg                                                                                                   <  7:285540/148‑1 (MQ=255)
                                                          cgcCCTGGGTGCTGGGCGCTGGCGTGATCTGGTGGGGGCGGGCCTGCTGGCTGGTGCTTACCTAGCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCGTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTTCCGGCCTGGCAAGGGCTg                                                                                         <  5:10143/149‑1 (MQ=255)
                                                             ccTGGGTGCTGGGCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGtg                                                                                      >  4:180449/1‑149 (MQ=255)
                                                                        ggCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCa                                                                           >  4:26429/1‑149 (MQ=255)
                                                                           gcTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGACCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTg                                                                        <  4:3013/149‑1 (MQ=255)
                                                                               ccGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGTCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGccc                                                                    <  4:107801/149‑1 (MQ=255)
                                                                               ccGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGccc                                                                    >  5:177009/1‑149 (MQ=255)
                                                                                 gggATCTGGGGGGGGCGTGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTTATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCtt                                                                  <  8:237378/147‑1 (MQ=255)
                                                                                         gtggGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTggg                                                               <  1:35240/144‑1 (MQ=255)
                                                                                         gtggGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTGGGCaa                                                            <  6:177009/147‑1 (MQ=255)
                                                                                         gtggGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTGGGCAACt                                                          >  8:236923/1‑149 (MQ=255)
                                                                                            ggggCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTGGGCAACTGGc                                                       <  7:236923/149‑1 (MQ=255)
                                                                                            ggggCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTGGGCAACTGGc                                                       <  3:180449/149‑1 (MQ=255)
                                                                                               gCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTGGGCAACTGGCTGa                                                    <  4:74221/149‑1 (MQ=255)
                                                                                                cTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTGGGCAACTGGCTGAt                                                   <  2:76124/149‑1 (MQ=255)
                                                                                                     cACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGTCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTGGGCAACTgg                                                        >  5:19845/1‑139 (MQ=255)
                                                                                                     cACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGTCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTGGGCAACTgg                                                        <  6:19845/139‑1 (MQ=255)
                                                                                                     cACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTGGGCAACTGGCTGATCCTGt                                              >  6:243519/1‑149 (MQ=255)
                                                                                                             tggTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTgg                                                                >  8:2897/1‑123 (MQ=255)
                                                                                                             tggTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTgg                                                                <  7:2897/123‑1 (MQ=255)
                                                                                                             tggTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTGGGCAACTGGCTGATCCTGTCGGTCATg                                      >  3:281132/1‑149 (MQ=255)
                                                                                                             gggggCTTACCTACCCGCGGAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTTGTGCTGCTCAAGCACTGGCCCTTGGGCAACTGGCTGATCCTGTGGGTCATg                                      <  8:248854/145‑1 (MQ=255)
                                                                                                                   gttcctaccCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCt                                                                   <  6:140863/111‑1 (MQ=255)
                                                                                                                     acctaccCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAg                                                                             <  2:119446/102‑1 (MQ=255)
                                                                                                                     acctaccCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAg                                                                             >  1:119446/1‑102 (MQ=255)
                                                                                                                       ctaccCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTGGGCAACTGGCTGATCCTGTCGGTCATGGTGCTGgtgt                            >  4:156171/1‑149 (MQ=255)
                                                                                                                          cccGCGCAGCAGCGACCTGTGGGCCTGTGCGGCGTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTGGGCAACTGGCTGATCCTGTCGGTCATGGTGCTGGtgtg                           <  6:226142/147‑1 (MQ=255)
                                                                                                                           ccGCGCAGCAGCGACCTGTGGGTCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTGGGCAACTGGCTGATCCTGTCGGTCATGGTGCTGGTGTGGgc                        >  4:129067/1‑149 (MQ=255)
                                                                                                                                gcagcaGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTGGGCAACTGGCTGATCCTGTCGGTCATGGTGCTGGTGTGGgccgcc                    <  2:139127/149‑1 (MQ=255)
                                                                                                                                gcagcaGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGGGCTGCTCAAGCACTGGCCCTTGGGCAACTGGCTGATCCTGTCGGTCATGGTGCTGGTGTGGgccgcc                    <  6:323861/149‑1 (MQ=255)
                                                                                                                                 cagcagCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTGGGCAACTGGCTGATCCTGTCGGTCATGGTGCTGGTGTGGgccgcc                    >  2:166487/1‑148 (MQ=255)
                                                                                                                                  agcagcGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTGGGCAACTGGCTGATCCTGTCGGTCATGGTGCTGGTGTGGGCCGCCGa                  >  3:276880/1‑149 (MQ=255)
                                                                                                                                              tgGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACt                                                                         >  2:54552/1‑82 (MQ=255)
                                                                                                                                              tgGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACt                                                                         <  1:54552/82‑1 (MQ=255)
                                                                                                                                              tgGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTGGGCAACTGGCTGATCCTGTCGGTCATGGTGCTGGTGTGGGCCGCCGACATCGGTGCGTa      >  5:322188/1‑149 (MQ=255)
                                                                                                                                                  ccAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTGGGCAACTGGCTGATCCTGTCGGTCATGGTGCTGGTGTGGGCCGCCGACATCGGTGCGTACTtc  >  1:135337/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCCTATGCAGCGGTGGTCGCGGGGGCGCTGATGCTGCTGCACATCCTTCCGGAACTGGCGCCCTGGGTGCTGGGCGCTGCCGTGATCTGGTGGGGGCTGGCCACCTGGCTGGTGCTTACCTACCCGCGCAGCAGCGACCTGTGGGCCAGTGCGGCCTGTCGGTTGTTGATCGGCCTGCTGGTGTTGCTGCCGGCCTGGCAAGGGCTGGTGCTGCTCAAGCACTGGCCCTTGGGCAACTGGCTGATCCTGTCGGTCATGGTGCTGGTGTGGGCCGCCGACATCGGTGCGTACTTC  >  NC_002947/1789451‑1789744

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: