Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,067,7830GA20.7% 54.2 / 10.8 29T91T (ACC→ACT) plsC1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (10/13);  new base A (3/3);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCTGCAGAACCCGGGTGGTCTTCTGTAAGGCCCTGCGCGCCTGATAGGCATTCTTGCGGTCGACCAGCACGTTGCCGCCCAGCCAGAACAGTTGGCCGAACAGCGGGATCCAGCCCAGGCTCTTCTTGCCGATGGCGACGGTGCGCTGCGGCACCACTTGGCCCAGCACGAACAGATCGAAGTTGGACTGGTGGTTGGCCACGATCACGCAGCCGGGCGGGTGGTCCCACAATGGGCCGACTTCGGCCTTGACCTCGATGCGCATCAGCCAGGTGGCCGGCAA  >  NC_002947/2067644‑2067926
                                                                                                                                           |                                                                                                                                               
ccTGCAGAACCCGGGTGGTCTTCTGTAAGGCCCTGCGCGCCTGATAGGCATTCTTGCGGTCGACCAGCACGTTGCCGCCCAGCCAGAACAGTTGGCCGAACAGCGGGATCCAGCCCAGGCTCTTCTTGCCGATGGCGACGGTGCGCTg                                                                                                                                         <  5:145362/148‑1 (MQ=255)
  tGCAGAACCCGGGTGGTCTTCTGTAAGGCCCTGCGCGCCTGATAGGCATTCTTGCGGTCGACCAGCACGTTGCCGCCCAGCCAGAACAGTTGGCCGAACAGCGGGATCCAGCCCAGGCTCTTCTTGCCGATGGCGACGGTGCGCTGCg                                                                                                                                       >  2:25342/1‑148 (MQ=255)
     aGAACCCGGGTGGTCTTCTGTAAGGCCCTGCGCGCCTGATAGGCATTCTTGCGGTCGACCAGCACGTTGCCGCCCAGCCAGAACAGTTGGCCGAACAGCGGGATCCAGCCCAGGCTCTTCTTGCCGATGGCGACGGGGCGCTGCGGc                                                                                                                                     <  6:139371/147‑1 (MQ=255)
           cGGGTGGTCTTCTGTAAGGCCCTGCGCGCCTGATAGGCATTCTTGCGGTCGACCAGCACGTTGCCGCCCAGCCAGAACAGTTGGCCGAACAGCGGGATCCAGCCCAGGCTCTTCTTGCCGATGGCGACGGTGCGCTGCGGCACCACtt                                                                                                                              <  3:32679/148‑1 (MQ=255)
            gggTGGTCTTCTGTAAGGCCCTGCGCGCCTGATAGGCATTCTTGCGGTCGACCAGCACGTTGCCGCCCAGCCAGAACAGTTGGCCGAACAGCGGGATCCAGCCCAGGCTCTTCTTGCCGATGGCGACGGTGCGCTGCGGCACCACTTg                                                                                                                             <  8:74991/148‑1 (MQ=255)
                  tcttctGTAAGGCCCTGCGCGCCTGATAGGCATTCTTGCGGTCGACCAGCACGTTGCCGCCCAGCCAGAACAGTTGGCCGAACAGCGGGATCCAGCCCAGGCTCTTCTTGCCGATGGCGACAGTGCGCTGCGGCACCACTTGGCCCa                                                                                                                        >  6:3730/1‑147 (MQ=255)
                  tcttctGTAAGGCCCTGCGCGCCTGATAGGCATTCTTGCGGTCGACCAGCACGTTGCCGCCCAGCCAGAACAGTTGGCCGAACAGCGGGATCCAGCCCAGGCTCTTCTTGCCGATGGCGACAGTGCGCTGCGGCACCACTTGGCCCAg                                                                                                                       >  6:5843/1‑148 (MQ=255)
                  tcttatGTAAGGCCCTGCGCGCCTGATAGGCATTCTTGCGGTCGACCAGCACGTTGCCGCCGAGTCAGAACAGTTGGACGAACAGCGGGATCCAGCACAGGCTCTTCTTGCCGATGGCGACAGAGCGCTGAGGCACCACTTGGGCCa                                                                                                                        >  6:3732/1‑147 (MQ=255)
                    ttctGTAAGGCCCTGCGCGCCTGATAGGCATTCTTGCGGTCGACCAGCACGTTGCCGCCCAGCCAGAACAGTTGGCCGAACAGCGGGATCCAGCCCAGGCTCTTCTTGCCGATGGCGACGGTGCGCTGCGGCACCACTTGGCCCAGCa                                                                                                                     >  5:8678/1‑148 (MQ=255)
                           aGGCCCTGCGCGCCTGATAGGCATTCTTGCGGTCGACCAGCACGTTGCCGCCCAGCCAGAACAGTTGGCCGAACAGCGGGATCCAGCCCAGGCTCTTCTTGCCGATGGCGACGGTGCGCTGCGGCACCACTTGGCCCAGCACGa                                                                                                                  >  5:115063/1‑144 (MQ=255)
                              cccTGCGCGCCTGATAGGCATTCTTGCGGTCGACCAGCACGTTGCCGCCCAGCCAGAACAGTTGGCCGAACAGCGGGATCCAGCCCAGGCTCTTCTTGCCGATGGCGACGGTGCGCTGCGGCACCACTTGGCCCAGCACGAACAGAt                                                                                                            >  6:144916/1‑147 (MQ=255)
                              cccTGCGCGCCTGATAGGCATTCTAGCGGTCGACCAGCCCGTTGCCGCCCAGCCAGAACAGTTGGCCGACCAGCGGGATCCAGCCCACGCTCTTCTTGCCGCTGGCGACGGTGCGCTGCGCCACCCCAAGGCCCCGCACCAACACAt                                                                                                            >  2:157609/1‑147 (MQ=255)
                                    gcgcCTGATAGGCATTCTTGCGGTCGACCAGCACGTTGCCGCCCAGCCAGAACAGTTGGCCGAACAGCGGGATCCAGCCCAGGCTCTTCTTGCCGATGGCGACGGTGCGCTGCGGCACCACTTGGCCCAGCACGAACAGATCGAAg                                                                                                       <  1:25342/146‑1 (MQ=255)
                                    gcgcCTGATAGGCATTCTTGCGGTCGACCAGCACGTTGCCGCCCAGCCAGAACAGTTGGCCGAACAGCGGGATCCAGCCCAGGCTCTTCTTGCCGATGGCGACGGTGCGCTGCGGCACCACTTGGCCCAGCACGAACAGATCGAAGtt                                                                                                     <  7:112398/148‑1 (MQ=255)
                                                               ccAGCACGGGGCCGCCCAGCCAGAACAGTTGGCCGAACAGCGGGATCCAGCCCAGGCTCTTCTGGCCGATGGCGACGGTGCGCTGCGGCACCACTTGGCCCAGCACGAACAGAGCGAAGTTGGACTGGTGGTTGGCCACGATCACGCa                                                                          <  8:142504/148‑1 (MQ=255)
                                                                    aCGTTGCCGCCCAGCCAGAACAGTTGGCCGAACAGCGGGATCCAGCCCAGGCTCTTCTTGCCGATGGCGACGGTGCGCTGCGGCACCACTTGGCCCAGCACGAACAGATCGAAGTTGGACTGGTGGTTGGCCACGATCACGCAGCcgg                                                                     <  7:147560/148‑1 (MQ=255)
                                                                              ccagccagAACAGTTGGCCGAACAGCGGGATCCAGCCCAGGCGCTTCTTGCCGATGGCGACGGTGCGCTGCGGCACCACTTGGCCCAGCACGAACAGATCGAAGTTGGACTGGTGGTTGGCCACGATCACGCAGCCGGGCGGGTGGTc                                                           <  2:50037/148‑1 (MQ=255)
                                                                                     gAACAGTTGGCCGAACAGCGGGATCCAGCCCAGGCTCTTCTTGCCGATGGCGACGGTGCGCTGCGGCACCACTTGGCCCAGCACGAACAGATCGAAGTTGGACTGGTGGTTGGCCACGATCACGCAGCCGGGCGGGTGGTCCCACAAt                                                    >  2:17360/1‑148 (MQ=255)
                                                                                                cGAACAGCGGGATCCAGCCCAGGCTCTTCTTGCCGATGGCGACGGTGCGCTGCGGCACCACTTGGCCCAGCACGAACAGATCGAAGTTGGACTCGTCGTTGCCCACGATCACCCAGCCGGCCGGGTGGTCCCACAATCGGACGCCTTc                                         >  6:93324/1‑148 (MQ=255)
                                                                                                     aGCGGGATCCAGCCCAGGCTCTTCTTGCCGATGGCGACGGTGCGCTGCGGCACCACTTGGCCCAGCACGAACAGATCGAAGTTGGACTGGTGGTTGGCCACGATCACGCAGCCGGGCGGGTGGTCCCACAATGGGCCGACTTCGGCCt                                    >  1:19536/1‑148 (MQ=255)
                                                                                                        gggATCCAGCCCAGGCTCTTCTTGCCGATGGCGACGGTGCGCTGCGGCACCACTTGGCCCAGCACGAACAGATCGAAGGTGGACTGGTGGTTGGCCACGATCACGCAGCCGGGCGGGTGGTCCCACAATGGGCCGACTTCGGCCTTGa                                 <  5:144916/148‑1 (MQ=255)
                                                                                                          gATCCAGCCCAGGCTCTTCTTGCCGATGGCGACGGTGCGCTGCGGCACCACTTGGCCCAGCACGAACAGATCGAAGTTGGACTGGTGGTTGGCCACGATCACGCAGCCGGGCGGGTGGTCCCACAATGGGCCGACTTCGGCCTTGAcc                               <  2:144451/148‑1 (MQ=255)
                                                                                                               aGCCCAGGCTCTTCTTGCCGATGGCGACGGTGCGCTGCGGCTCCACTTGGCCGTGGACGAACAGATCGAAGTTGGACTGGTGGTTGGCCACGATCACGCAGCCGGGCGGGTGGTCCCACAATGGGCCGACTTCGGCCTTGACCTCGAt                          <  5:93324/148‑1 (MQ=255)
                                                                                                                             ttGCCGATGGCGACGGTGCGCTGCGGCACCACTTGGCCCAGCACGAACAGATCGAAGTTGGACTGGTGGTTGGCCACGATCACGCAGCCGGGCGGGTGGTCCCACAATGGGCCGACTTCGGCCTTGACCTCGATGCGCATCAGCCAgg            >  5:7025/1‑148 (MQ=255)
                                                                                                                                  gATGGCGACGGTGCGCTGCGGCACCACTTGGCCCAGCACGAACAGATCGAAGTTGGACTGGTGGTTGGCCACGATCACGCAGCCGGGCGGGTGGTCCCACAATGGGCCGACTTCGGCCTTGACCTCGATGCGCATCAGCCAGGTGGcc       <  2:110416/148‑1 (MQ=255)
                                                                                                                                    tGGCGACGGTGCGCTGCGGCACCACTTGGCCCAGCACGAACAGATCGAAGTTGGACTGGTGGTTGGCCACGATCACGCAGCCGGGCGGGTGGTCCCCCAATGGGCCGACTTCGGCCTTGCCCTCGATGCGCATCAGCCAGGTGGCCg      >  5:76156/1‑147 (MQ=255)
                                                                                                                                       cGACAGTGCGCTGCGGCACCACTTGGCCCAGCACGAACAGATCGAAGTTGGACTGGTGGTTGGCCACGATCACGCAGCCGGGCGGGTGGTCCCACAATGGGCCGACTTCGGCCTTGACCTCGATGCGCATCAGCCAGGTGGCCGGCaa  <  5:5843/148‑1 (MQ=255)
                                                                                                                                       cGACAGTGCGCTGCGGCACCACTTGGCCCAGCACGAACAGATCGAAGTTGGACTGGTGGTTGGCCACGATCACGCAGCCGGGCGGGTGGTCCCACAATGGGCCGACTTCGGCCTTGACCTCGATGCGCATCAGCCAGGTGGCCGGCaa  <  5:3730/148‑1 (MQ=255)
                                                                                                                                       agACAGTGCGCTGCGCCACCACTAGGCCCAGCACCAACAGATCGAAGTTGCACTGCTTGATCGCCACGATCACGAAGCCGGGCGGGTGGTCCCACAATGGGCAGACTTCGGACTTGACCTCGATGCGCATCAGCCAGGTGGCCGGCaa  <  5:3732/147‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                               
CCTGCAGAACCCGGGTGGTCTTCTGTAAGGCCCTGCGCGCCTGATAGGCATTCTTGCGGTCGACCAGCACGTTGCCGCCCAGCCAGAACAGTTGGCCGAACAGCGGGATCCAGCCCAGGCTCTTCTTGCCGATGGCGACGGTGCGCTGCGGCACCACTTGGCCCAGCACGAACAGATCGAAGTTGGACTGGTGGTTGGCCACGATCACGCAGCCGGGCGGGTGGTCCCACAATGGGCCGACTTCGGCCTTGACCTCGATGCGCATCAGCCAGGTGGCCGGCAA  >  NC_002947/2067644‑2067926

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