Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 68,591 T→C 8.5% I228V (ATT→GTT)  PP_0058 ← 1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294768,5910TC 8.5% 131.1 / 2.7 47I228V (ATT→GTT) PP_00581‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base T (22/21);  new base C (2/2);  total (24/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.39e-01

TCAAGCATTTAATCCAAAAACTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCA  >  NC_002947/68442‑68737
                                                                                                                                                     |                                                                                                                                                  
tCAAGCATTTAATCCAAAAACTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAAt                                                                                                                                                    >  2:437494/1‑150 (MQ=255)
 cAAGCATTTAATCCAAAAACTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATg                                                                                                                                                   <  1:437494/150‑1 (MQ=255)
  aaGCATTTAATCCAAAAACTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATgg                                                                                                                                                  >  1:768347/1‑150 (MQ=255)
  aaGCATTTAATCCAAAAACTCCACACAAGTTATCCACAGCAGGTCAAGCCATCTGCTCTTGAGGCATTTCCGCAGCTCGCGCCAGCGATCCCATGGCCCGCTGGGTCCCTTCGTTCCAAGCTGCTGCGCGGTCGTTGAGTTCGGCAATgg                                                                                                                                                  <  2:192203/150‑1 (MQ=255)
    gCATTTAATCCAAAAACTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCa                                                                                                                                                >  2:1408042/1‑150 (MQ=255)
          aaTCGAAAAACTAGACACGAGTGATCCGCACAAGGTTAAGCCATCTGCTCCTGAGGTTTTTCCGCAGGTGGCGGCACCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGGAATGGCACGTGGc                                                                                                                                          <  2:63916/150‑1 (MQ=255)
           agccAAATACTGCTCACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGcc                                                                                                                                         <  2:627482/148‑1 (MQ=255)
            tCCAAAAACTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGccc                                                                                                                                        <  2:885957/150‑1 (MQ=255)
               aaaaaCTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGt                                                                                                                                     <  1:1408042/150‑1 (MQ=255)
                      gCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTg                                                                                                                              <  2:385973/150‑1 (MQ=255)
                        acacaAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGCCAACGGCACGGGGCCACGTGCCATTgcc                                                                                                                            <  1:758234/150‑3 (MQ=255)
                        acacaAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGCCAACGGCACGGGGCCACGTGCCATTgcc                                                                                                                            <  1:1305403/150‑3 (MQ=255)
                              gTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCCCGTGGCCCCGTGCCGTTGGCGTACAt                                                                                                                      >  2:1397393/1‑150 (MQ=255)
                               ttATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCGTGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATc                                                                                                                     <  2:1229725/150‑1 (MQ=255)
                                  tCCACAGAAGGTCAATCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGc                                                                                                                  <  1:128803/150‑1 (MQ=255)
                                     acaGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTGGGCGTACATCGGCTCa                                                                                                               >  2:926206/1‑150 (MQ=255)
                                           ggTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAAtt                                                                                                         <  2:768347/150‑1 (MQ=255)
                                                  ccATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACtt                                                                                                  >  2:858440/1‑150 (MQ=255)
                                                  ccATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACtt                                                                                                  >  1:19488/1‑150 (MQ=255)
                                                          tCCTGAGGTGTTTCAGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTa                                                                                          <  1:1360659/150‑1 (MQ=255)
                                                           ccTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTCCATCGGCTCACCAATTACCACTTCGATTGTAc                                                                                         >  1:1263414/1‑150 (MQ=255)
                                                              gAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCgg                                                                                      >  2:848388/1‑150 (MQ=255)
                                                                 gTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCAc                                                                                   >  2:685684/1‑150 (MQ=255)
                                                                      tCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTg                                                                              <  1:926206/150‑1 (MQ=255)
                                                                             cTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAAc                                                                       <  2:10382/150‑1 (MQ=255)
                                                                                       cGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTggg                                                             >  2:1003746/1‑150 (MQ=255)
                                                                                            ccATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGa                                                        <  1:1003746/150‑1 (MQ=255)
                                                                                                ggCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACtt                                                    <  1:848388/150‑1 (MQ=255)
                                                                                                 gCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCa                                                                                                               <  2:1193401/90‑1 (MQ=255)
                                                                                                 gCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCa                                                                                                               >  1:1193401/1‑90 (MQ=255)
                                                                                                  cccGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGaa                                                       >  1:52030/1‑145 (MQ=255)
                                                                                                  cccGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGaa                                                       <  2:52030/145‑1 (MQ=255)
                                                                                                  cccGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGc                                                  >  1:321661/1‑150 (MQ=255)
                                                                                                        gggTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGt                                            >  1:1080562/1‑150 (MQ=255)
                                                                                                           tCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTtgt                                         >  2:325109/1‑150 (MQ=255)
                                                                                                                   ttCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATc                                 >  2:682666/1‑150 (MQ=255)
                                                                                                                   ttCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGACAGAACTTGCCGGCGTTGTGGGCAATc                                 <  2:19488/150‑1 (MQ=255)
                                                                                                                      cAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGgc                              <  1:682666/150‑1 (MQ=255)
                                                                                                                       aaGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGgca                             >  1:435582/1‑150 (MQ=255)
                                                                                                                                    gTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAAc                                                                       <  1:862083/95‑1 (MQ=255)
                                                                                                                                    gTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAAc                                                                       >  2:862083/1‑95 (MQ=255)
                                                                                                                                     tCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCg               >  1:1265253/1‑150 (MQ=255)
                                                                                                                                                ggCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGc    >  2:951391/1‑150 (MQ=255)
                                                                                                                                                 gCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGcc   >  1:4436/1‑150 (MQ=255)
                                                                                                                                                  cAACGGCACGGGGCCACGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCa  >  1:962225/1‑150 (MQ=255)
                                                                                                                                                  cAACGGCACGGGGCCACGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCa  >  2:848107/1‑150 (MQ=255)
                                                                                                                                                    aTGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCg                                                <  2:346134/102‑1 (MQ=255)
                                                                                                                                                    aTGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCg                                                >  1:346134/1‑102 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TCAAGCATTTAATCCAAAAACTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCA  >  NC_002947/68442‑68737

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