Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,162,815 A→T 10.1% S223T (TCG→ACG)  plsC ← 1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,162,8150AT10.1% 182.8 / 5.2 57S223T (TCG→ACG) plsC1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (22/29);  new base T (3/3);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATTTGCTCCATGTAAAACTGGCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACGGGAATAAT  >  NC_000913/3162668‑3162962
                                                                                                                                                   |                                                                                                                                                   
atatTTGCTCCATGTAAAACTGGCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGaa                                                                                                                                                    <  1:11000/149‑1 (MQ=255)
atatTTGCTCCATGTAAAACTGGCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGaa                                                                                                                                                    >  2:554222/1‑149 (MQ=255)
   tttGCTCCATGTAAAACTGGCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGa                                                                                                                                                 >  2:565343/1‑149 (MQ=255)
       cTCCATGTAAAACTGGCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAAt                                                                                                                                             <  2:55266/149‑1 (MQ=255)
              tAAAACTGGCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGcagc                                                                                                                                      <  2:273833/149‑1 (MQ=255)
              tAAAACTGGCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGcagc                                                                                                                                      <  2:567331/149‑1 (MQ=255)
                     ggCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCa                                                                                                                               <  3:79997/149‑1 (MQ=255)
                                tACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATc                                                                                                                    <  1:76979/149‑1 (MQ=255)
                                   cGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGCTTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCttt                                                                                                                 <  1:510387/149‑1 (MQ=255)
                                   cGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCttt                                                                                                                 <  4:264403/149‑1 (MQ=255)
                                   cGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCttt                                                                                                                 <  1:308347/149‑1 (MQ=255)
                                   cGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCttt                                                                                                                 <  1:554222/149‑1 (MQ=255)
                                   cGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTCCCATTATCAAACGACAATGGGCAGCCAGCTCACGAACCTGATCgtt                                                                                                                 <  4:340391/149‑1 (MQ=255)
                                       cAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCa                                                                                                             <  2:253711/149‑1 (MQ=255)
                                        aGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCAt                                                                                                            >  2:50725/1‑149 (MQ=255)
                                              gtggTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTAATCGAGCTCGGCGATTGTTTGTTCCATTATCGAACGACAATGGGCAGCGCGCTCACGAACCTGATCTTTGCCATactgac                                                                                                      >  4:231331/1‑149 (MQ=255)
                                                 ttCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATActgactga                                                                                                   <  3:170234/149‑1 (MQ=255)
                                                   ccccGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTGTCGTTCGATAATggaac                                                                                                                                        <  3:285702/110‑4 (MQ=255)
                                                   ccccGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTGTCGTTCGATAATggaac                                                                                                                                        <  3:561599/110‑4 (MQ=255)
                                                   ccccGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTGTCGTTCGATAATggaac                                                                                                                                        >  4:285702/1‑107 (MQ=255)
                                                   ccccGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTGTCGTTCGATAATggaac                                                                                                                                        >  4:561599/1‑107 (MQ=255)
                                                   ccccGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTGTCGTTCGATAATggaac                                                                                                                                        <  1:241063/110‑4 (MQ=255)
                                                   acccGCAAAAAGACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTGTCGTTCGATAATggaac                                                                                                                                        >  2:241063/2‑107 (MQ=255)
                                                    cccGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGt                                                                                                >  4:53956/1‑149 (MQ=255)
                                                       gCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCaa                                                                                             >  1:241745/1‑149 (MQ=255)
                                                       gCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCaa                                                                                             >  1:578336/1‑149 (MQ=255)
                                                         aaaaaTACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAAtt                                                                                           <  3:404413/149‑1 (MQ=255)
                                                            aaTACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTggc                                                                                        <  2:346648/149‑1 (MQ=255)
                                                             aTACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTggcg                                                                                       >  3:38670/1‑149 (MQ=255)
                                                             aTACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTggcg                                                                                       >  1:395015/1‑149 (MQ=255)
                                                             aTACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTggcg                                                                                       <  3:272351/149‑1 (MQ=255)
                                                                cAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGgca                                                                                    <  2:179317/149‑1 (MQ=255)
                                                                    ttGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGACCGACAATGGGCAGCCAGCTCACGAGCCTGATCTTTGCCATACTGACTGCCGCCAATTGACGtcagtca                                                                                <  4:132561/149‑4 (MQ=255)
                                                                        cGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCg                                                                            >  2:574690/1‑149 (MQ=255)
                                                                        cGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCg                                                                            >  3:138690/1‑149 (MQ=255)
                                                                           cTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACa                                                                         <  4:169222/149‑1 (MQ=255)
                                                                                 cTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCAcc                                                                   >  1:75573/1‑149 (MQ=255)
                                                                                     tCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGAc                                                               <  2:226338/149‑1 (MQ=255)
                                                                                      ccGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGAcc                                                              >  2:417920/1‑149 (MQ=255)
                                                                                      ccGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGAcc                                                              >  2:409456/1‑149 (MQ=255)
                                                                                           ggcTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTtgt                                                         >  4:231421/1‑149 (MQ=255)
                                                                                                 gcgTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCTATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTc                                                   <  1:477603/149‑1 (MQ=255)
                                                                                                 gcgTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTc                                                   >  1:105510/1‑149 (MQ=255)
                                                                                                   gTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCga                                                 <  4:209867/149‑1 (MQ=255)
                                                                                                    ttCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCgat                                                <  3:53956/149‑1 (MQ=255)
                                                                                                    ttCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCgat                                                <  1:565343/149‑1 (MQ=255)
                                                                                                            cttcttTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCgattaagatta                                        >  2:334806/1‑149 (MQ=255)
                                                                                                              tcttTATCGAGCTCGGCGATTTTTTGATCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAAt                                      <  3:299062/149‑1 (MQ=255)
                                                                                                                tttATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCgattaaga                                           >  3:325528/1‑142 (MQ=255)
                                                                                                                tttATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCgattaaga                                           <  4:325528/142‑1 (MQ=255)
                                                                                                                  tATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGgca                                                                                    >  2:157406/1‑99 (MQ=255)
                                                                                                                  tATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGgca                                                                                    <  1:157406/99‑1 (MQ=255)
                                                                                                                   aTCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTAt                                 <  3:388262/149‑1 (MQ=255)
                                                                                                                    tCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTAtt                                <  4:38670/149‑1 (MQ=255)
                                                                                                                          tCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGt                          >  1:372492/1‑149 (MQ=255)
                                                                                                                                           cATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACgg         >  3:478083/1‑149 (MQ=255)
                                                                                                                                           cATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCAGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACgg         <  1:105186/149‑1 (MQ=255)
                                                                                                                                            aTTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACggg        >  4:36112/1‑149 (MQ=255)
                                                                                                                                                  gAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACGGGaataat  >  2:284803/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATATTTGCTCCATGTAAAACTGGCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACGGGAATAAT  >  NC_000913/3162668‑3162962

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