Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,092,574 G→C 10.9% intergenic (‑270/+67) astC ← / ← PP_4482 acetylornithine aminotransferase/AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,092,5740GC10.9% 101.3 / 2.6 35intergenic (‑270/+67)astC/PP_4482acetylornithine aminotransferase/AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (15/16);  new base C (2/2);  total (17/18)
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.16e-01

TCGGCAATGCGACAAACACTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTG  >  NC_002947/5092436‑5092720
                                                                                                                                          |                                                                                                                                                  
tCGGCAATGCGACAAACACTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGgcggcg                                                                                                                                           <  8:78452/148‑1 (MQ=255)
     aaTGCGACAAACACTGTCGCAATAGCGCAGTGTACGTGAGAGAGGTGGCTTGGTGAACGGGTGCTCACACATGATTGTGAGGGGGCAGGATGTACGTCGATGTTGATGTGTTGACTGGGAGATCGAGCGCCGCGTGGGCGGGGCTCGAt                                                                                                                                     <  8:206352/149‑1 (MQ=255)
         cGACCAACCCTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGGGAACGGGGTCTCACCCAGGATTGTGAGGGGGCCGGATGTACGGCGATGTTGATGTGTTGACGGGCAAATCGAGCGCCGCCCACGCGGCGCTCGAtctcc                                                                                                                                 <  6:218952/149‑5 (MQ=37)
                 acTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAac                                                                                                                         <  8:329723/149‑1 (MQ=255)
                   tGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAaaac                                                                                                                       >  3:38977/1‑149 (MQ=255)
                              cgcAGGCTACGGGAGAGAGGGGGCGGGGGCAACGGGCGCCCACACAGGATCGGCCGGGGCCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAgg                                                                                                            <  3:318809/149‑1 (MQ=255)
                                                 ggggCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCGCGATTTGCGCCCCAACACCCACCCCAAGGCATACCCCCCGCacccccc                                                                                         >  6:244795/1‑149 (MQ=255)
                                                          ctatctggTGCTCACATAGGTTTGTGAGTGGTCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCTAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCACAACAGCAACCCCAAGGCATACACCCATAACCCACGCGATAtct                                                                                 <  4:262062/142‑1 (MQ=255)
                                                           gAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAAtct                                                                                 >  3:191548/1‑147 (MQ=255)
                                                                ggTTCTCACACAGGATTGTGAGGGGGCAGGTTGTACGGCGATGTTGATGTGTTGACGGTGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCa                                                                               <  3:72538/144‑1 (MQ=255)
                                                                ggTGCTCACACAGGATTGTGAGGCGGCAGGAAGTACGGCGATGTTCATCTGATGACGGGGAGATCGAGAGCCGCCTGCGCGGCGCTCGATTTGCCCCCCAACACCAACCCCAAAGCATACACCCAGAACACACCCGAAATCTAAgaga                                                                           >  4:72538/1‑145 (MQ=255)
                                                                 gTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACAACCCGAACCCACCCGAAAtct                                                                                 >  3:262062/1‑141 (MQ=255)
                                                                   gCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCttt                                                                        >  8:176638/1‑148 (MQ=255)
                                                                       acactGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGc                                                                   >  5:73205/1‑149 (MQ=255)
                                                                        cacaGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAAtc                                                                                  <  4:138850/133‑1 (MQ=255)
                                                                        cacaGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAACGCATACACCCCGAACCCACCCGAAAtc                                                                                  >  3:138850/1‑133 (MQ=255)
                                                                               ttGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGcc                                                            >  8:5848/1‑148 (MQ=255)
                                                                                tgtgAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCg                                                           <  4:191548/148‑1 (MQ=255)
                                                                                   gAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCCCACGCGGCGCTc                                                                                                                                        <  2:122844/68‑1 (MQ=37)
                                                                                   gAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCCCACGCGGCGCTc                                                                                                                                        >  1:122844/1‑68 (MQ=37)
                                                                                               tGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACgg                                           <  7:176638/149‑1 (MQ=255)
                                                                                                    ggCGATGTTGATGTGTTGACGGGGAGAGCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGgctgc                                      <  6:17930/149‑1 (MQ=255)
                                                                                                      cGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAAc                                                                                                                   >  6:99622/1‑70 (MQ=255)
                                                                                                      cGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAAc                                                                                                                   <  5:99622/70‑1 (MQ=255)
                                                                                                       gATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGaa                                                                                     <  6:320576/99‑1 (MQ=255)
                                                                                                       gATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGaa                                                                                     >  5:320576/1‑99 (MQ=255)
                                                                                                            tGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGgctgctgc                                   <  7:335767/144‑1 (MQ=255)
                                                                                                            tGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGgctgctgc                                   >  8:335767/1‑144 (MQ=255)
                                                                                                                  gtTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACCCCCCGCACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGa                        >  6:166721/1‑149 (MQ=255)
                                                                                                                     gACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGAtt                      <  6:73205/148‑1 (MQ=255)
                                                                                                                      aCGGGGAGATCGAGCGCCGCGGGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTg                     <  5:244795/148‑1 (MQ=255)
                                                                                                                                cGAGCGCCGCGTGGGCGGCGTTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGACCCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCcg          >  5:198387/1‑149 (MQ=255)
                                                                                                                                cGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCcg          >  4:194296/1‑149 (MQ=255)
                                                                                                                                     gccgcGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGgcg                                                              <  4:7588/92‑1 (MQ=255)
                                                                                                                                     gccgcGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGgcg                                                              >  3:7588/1‑92 (MQ=255)
                                                                                                                                        gcgTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTg  >  2:6838/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
TCGGCAATGCGACAAACACTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTG  >  NC_002947/5092436‑5092720

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: