Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,068,194 A→T 17.5% intergenic (‑96/+25) 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,1940AT17.5% 77.9 / 5.0 33intergenic (‑96/+25)plsC/PP_18451‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase/enoyl‑CoA hydratase/isomerase family protein
Reads supporting (aligned to +/- strand):  ref base A (13/14);  new base T (3/3);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAACGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTC  >  NC_002947/2068045‑2068332
                                                                                                                                                      |                                                                                                                                          
gcCAACGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCa                                                                                                                                            >  2:175683/1‑150 (MQ=255)
       gCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCcgc                                                                                                                                     >  1:44694/1‑150 (MQ=255)
         aaGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGa                                                                                                                                   >  2:1554/1‑150 (MQ=255)
            aaaaTGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGaaaa                                                                                                                                <  2:44694/150‑1 (MQ=255)
            aaaaTGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGaaaa                                                                                                                                <  2:113090/150‑1 (MQ=255)
             aaaTGCAGCCACAGCAGGAACATACGCAATGAATAAGGCATGGGACGCCTCACAGCAGACGGTCGGGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTTCGCGGGCTTGCCCGCGaaag                                                                                                                                <  2:350174/150‑2 (MQ=255)
                tGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGc                                                                                                                            <  2:133678/150‑1 (MQ=255)
                      cccAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGc                                                                                                                      <  1:1554/150‑1 (MQ=255)
                                cgatctgCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGACGCGCAAGGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCaa                                                                                                            <  2:22327/144‑1 (MQ=255)
                                     cGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtgc                                                                                                       >  1:772136/1‑150 (MQ=255)
                                      gCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCaa                                                                                                            >  2:556370/1‑144 (MQ=255)
                                      gCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCaa                                                                                                            <  1:556370/144‑1 (MQ=255)
                                       cAATGAATAAAGCATGGTACGACTCACATGCGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAatagatcg                                                                                                     >  1:22327/1‑145 (MQ=255)
                                                               cacaCCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTc                                                                             <  2:772136/150‑1 (MQ=255)
                                                                acacCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTcc                                                                            >  1:209425/1‑150 (MQ=255)
                                                                    cgagcGGTGGGGGGGAAGTGTAGGGGTGTGCACCGGCGGGGGTGAAGCTTGGCGGAGGCTTCGTGCGGCTCTTTCGGGGGTTAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAGTATGCTGCTGATCCAGCTCGATCGCCTTGTCCcgcg                                                                        <  2:342078/146‑1 (MQ=255)
                                                                                gggCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGCTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtgctg                                                                                                     <  2:61491/109‑1 (MQ=255)
                                                                                gggCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGCTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtgctg                                                                                                     >  1:61491/1‑109 (MQ=255)
                                                                                       gtACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTTCGCGGGCTTGCCGGCGaaag                                                                                                                                >  1:618903/1‑75 (MQ=25)
                                                                                       gtACGGGTGTGCACTGGCTGGGGCAAACCTTGCTGCAGGCCTGTTACGGCCCTTTTCGCGGGCTTGCCAGCGaaag                                                                                                                                <  2:618903/76‑2 (MQ=25)
                                                                                        tACGGGTGTGCACTGGCTGGGGCAAAGCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTgcgc                                                    >  1:504181/1‑150 (MQ=255)
                                                                                               gtgCACTGGCTGTGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTCGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCa                                             <  1:618883/150‑1 (MQ=255)
                                                                                                 gCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGc                                           >  1:207392/1‑150 (MQ=255)
                                                                                                            ggCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCtt                                <  2:467995/150‑1 (MQ=255)
                                                                                                             gCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATc                                                                                     >  2:714477/1‑96 (MQ=255)
                                                                                                             gCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATc                                                                                     <  1:714477/96‑1 (MQ=255)
                                                                                                                   cTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGcc                         >  2:483201/1‑150 (MQ=255)
                                                                                                                    ttGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCa                        <  2:209425/150‑1 (MQ=255)
                                                                                                                            ggCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTc                <  2:504181/150‑1 (MQ=255)
                                                                                                                                      ggCCCTTTTCGCGGGCTTGCCCGCGAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCa      >  2:1510/1‑150 (MQ=255)
                                                                                                                                          cTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTc  >  1:675849/1‑150 (MQ=255)
                                                                                                                                           ttttCGCGGGCTGGCCCGCGAAAAGGGCCCGTACGGCTAGCCAATATGCTACTGATCCAGCAGGATCGCCTTGTCCAGCGCCTCCAACAGCACTTTTGGAACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCATCTTCTTCAGCTc  >  2:240664/2‑150 (MQ=255)
                                                                                                                                                    gCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGgt                                      <  1:579008/105‑1 (MQ=255)
                                                                                                                                                    gCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGgt                                      >  2:579008/1‑105 (MQ=255)
                                                                                                                                                      |                                                                                                                                          
GCCAACGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTC  >  NC_002947/2068045‑2068332

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: