Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 68,603 C→A 15.5% G224W (GGG→TGG)  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,6030CA15.5% 115.2 / 12.1 45G224W (GGG→TGG) PP_00581‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (21/17);  new base A (5/2);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

TCCAAAAACTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGCGGTACCGCCG  >  NC_002947/68454‑68751
                                                                                                                                                     |                                                                                                                                                    
tCCAAAAACTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGccc                                                                                                                                                      <  2:885957/150‑1 (MQ=255)
   aaaaaCTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGt                                                                                                                                                   <  1:1408042/150‑1 (MQ=255)
          gCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTg                                                                                                                                            <  2:385973/150‑1 (MQ=255)
            acacaAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGCCAACGGCACGGGGCCACGTGCCATTgcc                                                                                                                                          <  1:1305403/150‑3 (MQ=255)
            acacaAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGCCAACGGCACGGGGCCACGTGCCATTgcc                                                                                                                                          <  1:758234/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)
                                                                                                                                     ccaacggCACGGGGCCACGTGCCGTTGGCGTACATCGGCTCACCAATCACCACTTCGACTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGcc                 >  1:898446/6‑150 (MQ=255)
                                                                                                                                      cAACGGCACGGGGCCACGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCa                >  1:962225/1‑150 (MQ=255)
                                                                                                                                      cAACGGCACGGGGCCACGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCa                >  2:848107/1‑150 (MQ=255)
                                                                                                                                       aacggCACGGGGCCACGTGCCATTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGACTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAg               >  1:169391/4‑150 (MQ=255)
                                                                                                                                        acggCACGGGGCCACGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAgg              >  2:891117/3‑150 (MQ=255)
                                                                                                                                        aTGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCg                                                              >  1:346134/1‑102 (MQ=255)
                                                                                                                                        aTGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCg                                                              <  2:346134/102‑1 (MQ=255)
                                                                                                                                          ggCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGc            <  1:685684/150‑1 (MQ=255)
                                                                                                                                                 ggCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGCGGTAccg     >  1:1227206/1‑150 (MQ=255)
                                                                                                                                                   ccccGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGCGGTAccgcc   >  1:463724/1‑150 (MQ=255)
                                                                                                                                                    cccGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGCGGTAccgccg  >  2:818026/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCCAAAAACTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGCGGTACCGCCG  >  NC_002947/68454‑68751

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