Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,092,577 G→C 10.3% intergenic (‑273/+64) 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,5770GC10.3% 107.0 / 2.5 37intergenic (‑273/+64)astC/PP_4482acetylornithine aminotransferase/AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (16/17);  new base C (2/2);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

TCGGCAATGCGACAAACACTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCAC  >  NC_002947/5092436‑5092723
                                                                                                                                             |                                                                                                                                                  
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)
                                                                                                                                           tGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCAc  <  5:166721/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
TCGGCAATGCGACAAACACTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCAC  >  NC_002947/5092436‑5092723

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: