Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,068,195 A→T 26.9% intergenic (‑97/+24) plsC ← / ← PP_1845 1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase/enoyl‑CoA hydratase/isomerase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,068,1950AT26.9% 62.8 / 27.6 41intergenic (‑97/+24)plsC/PP_18451‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase/enoyl‑CoA hydratase/isomerase family protein
Reads supporting (aligned to +/- strand):  ref base A (17/13);  new base T (7/4);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

CCAACGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTC‑GATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCGCTGCG  >  NC_002947/2068046‑2068342
                                                                                                                                                      |                                                                                                                                                    
ccAACGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCAC‑TTTCGCGGGCaa                                                                                                                                                      <  2:19361/150‑1 (MQ=255)
           aaaaTGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGaaaa                                                                                                                                           <  2:906017/150‑1 (MQ=255)
            aaaTGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTTCGCGGGCTTGCCCGCGaaag                                                                                                                                           <  1:1088610/150‑2 (MQ=255)
                   gCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGt                                                                                                                                   <  2:1117114/150‑1 (MQ=255)
                     cccAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCTTGCCCGCGa                                                                                                                                              <  1:1174557/137‑1 (MQ=255)
                     cccAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCTTGCCCGCGa                                                                                                                                              >  2:1174557/1‑137 (MQ=255)
                          cagGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTa                                                                                                                            >  2:1241014/1‑150 (MQ=255)
                            ggAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCCTGCGGCTAGc                                                                                                                          >  1:422108/1‑150 (MQ=255)
                             aaaCAAACGCAATGAATAAAGCATGGTACGACTCACACCAGGCGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGcc                                                                                                                         <  2:222046/149‑1 (MQ=255)
                               aCATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCaa                                                                                                                       <  1:597438/150‑1 (MQ=255)
                                      cAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGaaaa                                                                                                                                           <  1:964744/123‑1 (MQ=255)
                                      cAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGaaaa                                                                                                                                           >  2:964744/1‑123 (MQ=255)
                                                      gTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTC‑GATc                                                                                               <  1:936363/150‑1 (MQ=255)
                                                       tACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCTTGCCCGCGaaag                                                                                                                                           <  1:600288/106‑2 (MQ=255)
                                                       tACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCTTGCCCGCGaaag                                                                                                                                           >  2:600288/1‑105 (MQ=255)
                                                             tCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATgctgct                                                                                                              <  1:1088618/129‑1 (MQ=255)
                                                             tCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATgctgct                                                                                                              >  2:1088618/1‑129 (MQ=255)
                                                               acacCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTC‑GATCGCCTTGTcc                                                                                      >  2:1023257/1‑150 (MQ=255)
                                                               acacCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGAAGGGGGACACCGATGTATAGGGCGGTTACGGCCC‑TTTCGGGGTCAAGCCCGCGACAACGGCCCGCGCCGCTAGCCAATATGATGCTGATCCAGCTA‑GATCGCCTTGTcc                                                                                      >  1:388493/1‑150 (MQ=255)
                                                                   cAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTC‑GACCGCCTTGTCCAgcg                                                                                  >  2:382434/1‑150 (MQ=255)
                                                                                              gtgCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCACGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTC‑GATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCa                                                       <  1:1241014/150‑1 (MQ=255)
                                                                                               tgCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTC‑GATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCAc                                                            >  2:889519/1‑144 (MQ=255)
                                                                                               tgCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTC‑GATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCAc                                                            <  1:889519/144‑1 (MQ=255)
                                                                                                 cACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCCGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTC‑GATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCt                                                    >  1:535900/1‑150 (MQ=255)
                                                                                                   ctggctggGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTC‑GATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTg                                                  >  2:219429/1‑150 (MQ=255)
                                                                                                     ggctggGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTC‑GATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGgt                                                <  1:1023257/150‑1 (MQ=255)
                                                                                                        tggGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATc                                                                                                          >  1:1412519/1‑90 (MQ=255)
                                                                                                        tggGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATc                                                                                                          <  2:1412519/90‑1 (MQ=255)
                                                                                                            gCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCGAATATCCT‑CTAGCCCGGCGCTGCGCTCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGTTGTTCTtg                                         >  1:1169639/1‑150 (MQ=255)
                                                                                                                aCCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTC‑GATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTgcg                                     >  1:394619/1‑150 (MQ=255)
                                                                                                                   ttGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTC‑GATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCa                                  >  2:906034/1‑150 (MQ=255)
                                                                                                                         cAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTC‑GATc                                                                                               <  2:567972/83‑1 (MQ=255)
                                                                                                                         cAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTC‑GATc                                                                                               >  1:567972/1‑83 (MQ=255)
                                                                                                                              cTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTC‑GATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGc                       >  1:193948/1‑150 (MQ=255)
                                                                                                                                     ggCCCTTTTCGCGGGCTTGCCCGCG‑AAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTC‑GATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCa                <  2:1036422/150‑1 (MQ=255)
                                                                                                                                     ggCCCTTTTCGCGGGCTTGCCCGCG‑AAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTC‑GATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCa                >  1:1036422/1‑150 (MQ=255)
                                                                                                                                         cttttCGCGGGCTTGCCCGCG‑AAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTC‑GATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTc            >  1:1211935/3‑150 (MQ=255)
                                                                                                                                         cttttCGCGGGCTTGCCCGCG‑AAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTC‑GATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTc            >  2:361870/3‑150 (MQ=255)
                                                                                                                                         cttttCGCGGGCTTGCCCGCG‑AAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGATC‑GATCGCCTTCTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCtt            >  1:64789/3‑149 (MQ=255)
                                                                                                                                         c‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTC‑GATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTc            >  1:591778/1‑150 (MQ=255)
                                                                                                                                          ttttCGCGGGCTTGCCCGCG‑AAAGAGCCTGTGCAGCTTGACCATGTGCTCCGGATTCAGGTC‑GAGCGCCTGGTCCAGCGCCGCCAGCAGCCCTTTGCGCAATTTCAGCTTGGTGTTCTTGTGCGCCAGCGTGTTTAGCTGCTTCAGCTcc           >  2:1085683/2‑150 (MQ=255)
                                                                                                                                                    cAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTC‑GATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCGCTGCg  <  2:388493/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CCAACGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTC‑GATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCGCTGCG  >  NC_002947/2068046‑2068342

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