Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,162,104 T→A 20.0% intergenic (‑48/+29) PP_3663 ← / ← pssA GGDEF domain‑containing protein/CDP‑diacylglycerol‑‑serine O‑phosphatidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,162,1040TA20.0% 45.2 / 3.1 24intergenic (‑48/+29)PP_3663/pssAGGDEF domain‑containing protein/CDP‑diacylglycerol‑‑serine O‑phosphatidyltransferase
Reads supporting (aligned to +/- strand):  ref base T (10/9);  new base A (3/2);  total (13/11)
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.96e-01

TATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAACC  >  NC_002947/4161964‑4162243
                                                                                                                                            |                                                                                                                                           
tATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATTTCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGc                                                                                                                                       <  1:89175/147‑1 (MQ=255)
                 gTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGc                                                                                                                     >  8:99745/1‑148 (MQ=255)
                         gggCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCaga                                                                                                            <  5:178598/149‑1 (MQ=255)
                             cTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGCCGACAGGATCAGGACATGCGATGACagcacac                                                                                                        >  6:224826/1‑144 (MQ=255)
                                  aCTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGAAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGc                                                                                                    <  1:191620/148‑1 (MQ=255)
                                               caTTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCaga                                                                                                            >  4:185347/1‑127 (MQ=255)
                                                          aGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCCGGCAAGCCAGCGACAGGGCCAga                                                                                                                           <  2:55965/101‑2 (MQ=37)
                                                          aGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCCGGCAAGCCAGCGACAGGGCCAga                                                                                                                           >  1:55965/1‑100 (MQ=37)
                                                                    aaaGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGGTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCACCGGCTCCAACCGGCTATAGCGCCACCg                                                                 >  2:118799/1‑149 (MQ=255)
                                                                         gCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCag                                                                                                             <  5:103507/100‑1 (MQ=255)
                                                                         gCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCag                                                                                                             >  6:103507/1‑100 (MQ=255)
                                                                                  agcaATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTgc                                                   <  7:99745/149‑1 (MQ=255)
                                                                                           aTAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCActtgct                                          >  7:65424/1‑149 (MQ=255)
                                                                                                   aTCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGt                                  >  4:269917/1‑149 (MQ=255)
                                                                                                      aaTGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCTGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCg                               <  5:224826/149‑1 (MQ=255)
                                                                                                                     ctgtactGACCCTGTCGCCGGCAAGCCAGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACATGCGCACTTGCGTGGGATGTTCGGCCTTCTCCCCCAgc                >  5:214752/6‑149 (MQ=38)
                                                                                                                     ctgtactGACCCTGTCGCCGGCAAGCCAGCGACAGGGTCAGTACAGGCGATATCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACtt                                                     >  7:138676/6‑112 (MQ=37)
                                                                                                                     ctgtactGACCCTGTCGCCGGCAAGCCAGCGACAGGGTCAGTACAGGCGATATCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACtt                                                     <  8:138676/107‑1 (MQ=37)
                                                                                                                       acTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAgcgc              <  8:65424/149‑1 (MQ=255)
                                                                                                                            gacccTGTCGCGGGAATTCCAGCGAAAGGGTAAGTACAGGCGATTTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCctct          <  6:214752/146‑1 (MQ=255)
                                                                                                                                 tGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGaa    >  3:58539/1‑149 (MQ=255)
                                                                                                                                   tCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGa                                     >  4:97581/1‑114 (MQ=255)
                                                                                                                                   tCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGa                                     <  3:97581/114‑1 (MQ=255)
                                                                                                                                   tCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAAcc  >  5:126838/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                           
TATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAACC  >  NC_002947/4161964‑4162243

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