Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,162,820 C→T 10.3% C221Y (TGT→TAT)  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,8200CT10.3% 172.3 / 6.5 58C221Y (TGT→TAT) plsC1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (23/29);  new base T (3/3);  total (26/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 = 9.84e-01

CTCCATGTAAAACTGGCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACGGGAATAATCGGGA  >  NC_000913/3162675‑3162967
                                                                                                                                                 |                                                                                                                                                   
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                                                                                            >  1:395015/1‑149 (MQ=255)
                                                      aTACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTggcg                                                                                            >  3:38670/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                                                     <  1:565343/149‑1 (MQ=255)
                                                                                             ttCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCgat                                                     <  3:53956/149‑1 (MQ=255)
                                                                                                     cttcttTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCgattaagatta                                             >  2:334806/1‑149 (MQ=255)
                                                                                                       tcttTATCGAGCTCGGCGATTTTTTGATCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAAt                                           <  3:299062/149‑1 (MQ=255)
                                                                                                         tttATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCgattaaga                                                <  4:325528/142‑1 (MQ=255)
                                                                                                         tttATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCgattaaga                                                >  3:325528/1‑142 (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)
                                                                                                                                                aCAATGGGCAGTCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACGGGAATAATCGGGa  >  1:156198/1‑149 (MQ=255)
                                                                                                                                                aCAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACGGGAATAATCGGGa  >  4:467369/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTCCATGTAAAACTGGCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACGGGAATAATCGGGA  >  NC_000913/3162675‑3162967

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: ‑