Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,903,076 T→G 14.1% intergenic (+41/+74) PP_3445 → / ← ilvA‑I hypothetical protein/threonine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,903,0760TG14.1% 102.4 / 9.6 42intergenic (+41/+74)PP_3445/ilvA‑Ihypothetical protein/threonine deaminase
Reads supporting (aligned to +/- strand):  ref base T (14/22);  new base G (4/2);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCGAGCCGGATACCGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCGGTACCCCATCTCGTC  >  NC_002947/3902936‑3903215
                                                                                                                                             |                                                                                                                                           
ccccGAGCCGGATACCGAGAACTTCAAGCAGCTCCGCCGACTTCAAACGCCTGTAAACGTTCCGCAAGTGCGTTGCCGCCTTTTTCATCCGCGAGCGTTGAAATTGGCTCCCCCAGTTCCCTCCAGGAGATCGCCATCCCCTGcgggcg                                                                                                                                      >  1:216055/1‑149 (MQ=255)
 cccGAGCCGGATACCGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTcgggcggg                                                                                                                                    <  1:80829/149‑1 (MQ=255)
     aGCCGGATACCGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGggcggggcgg                                                                                                                                <  2:133632/149‑1 (MQ=255)
          gATACCGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACg                                                                                                                           >  5:310097/1‑149 (MQ=255)
           aTACCGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCCGCCCGAAgggcagggcga                                                                                                                                <  5:293427/143‑2 (MQ=255)
           aTACCGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCCGCCCGAAgggcagggcga                                                                                                                                >  6:293427/1‑142 (MQ=255)
            tACCGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACg                                                                                                                           >  7:182632/1‑147 (MQ=255)
             aCCGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAg                                                                                                                        <  7:312639/149‑1 (MQ=255)
              cctagcaCTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGACACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGc                                                                                                                       <  4:212232/143‑1 (MQ=255)
              ccGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGc                                                                                                                       >  2:224967/1‑149 (MQ=255)
              ccGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGc                                                                                                                       >  5:288284/1‑149 (MQ=255)
               cGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCg                                                                                                                      <  1:224967/149‑1 (MQ=255)
               cGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCg                                                                                                                      <  6:288284/149‑1 (MQ=255)
                  gaACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAac                                                                                                                   <  8:182632/149‑1 (MQ=255)
                        cAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcaggc                                                                                                             >  8:191844/1‑149 (MQ=255)
                        cAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcagg                                                                                                              <  1:81457/148‑1 (MQ=255)
                        cAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcagg                                                                                                              <  8:337084/148‑1 (MQ=255)
                                  cGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcaggcggcaggcgg                                                                                                    >  5:63674/1‑149 (MQ=255)
                                  cGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcaggcggcaggcgg                                                                                                    <  5:348467/149‑1 (MQ=255)
                                      gACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAacac                                                                                                                 >  2:233179/1‑132 (MQ=255)
                                      gACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAacac                                                                                                                 <  1:233179/132‑1 (MQ=255)
                                       aCTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGggcggggcgg                                                                                                                                >  2:75441/1‑116 (MQ=255)
                                       aCTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGggcggggcgg                                                                                                                                <  1:75441/116‑1 (MQ=255)
                                        cTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGAc                                                                                              <  6:94295/149‑1 (MQ=255)
                                          tCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTa                                                                                                 >  5:182918/1‑144 (MQ=255)
                                          tCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTa                                                                                                 <  6:182918/144‑1 (MQ=255)
                                            aaaaGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGc                                                                                          >  8:181365/1‑149 (MQ=255)
                                            aaaaGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGc                                                                                          >  4:14914/1‑149 (MQ=255)
                                              aaGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGt                                                                                        <  6:187330/149‑1 (MQ=255)
                                              aaGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGt                                                                                        >  6:55082/1‑149 (MQ=255)
                                                           ttCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGcac                                                                           <  6:467/149‑1 (MQ=255)
                                                                  aGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGcgcg                                                                    <  3:14914/149‑1 (MQ=255)
                                                                       gTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGc                                                               <  7:181365/149‑1 (MQ=255)
                                                                           ccGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAg                                                           >  1:134000/1‑149 (MQ=255)
                                                                                  tttCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAg                                                     <  6:63674/148‑1 (MQ=255)
                                                                                           cgAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGcc                                           <  2:216055/149‑1 (MQ=255)
                                                                                                            tCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCGTCGGGCGGGGCGGTCACGCAGCGAACACCGGCGGCAGGCGGCTAGCCAGGCGTGCATTGGCCGCACAGGCGCGTCCGCCCAGGAACCGCGTGTACGCCGGGTTGGCAGTCGCGTc                          >  5:300293/1‑149 (MQ=255)
                                                                                                              gcgcAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTcc                         <  6:310097/148‑1 (MQ=255)
                                                                                                                             ggAGATCGCCCTGCCCTTCGGGCGTGGCGGTCACGCAGCGAACACAGTCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCGGTACCCCa         <  4:24143/149‑1 (MQ=255)
                                                                                                                                gaccgccccgcccGAAGGGCAGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCGGTACCCCAtct      >  8:273687/4‑149 (MQ=255)
                                                                                                                                 accgccccgcccGAAGGGCAGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTAc                   <  6:119353/133‑1 (MQ=255)
                                                                                                                                 accgccccgcccGAAGGGCAGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTAc                   >  5:119353/3‑135 (MQ=255)
                                                                                                                                     ccctgccctTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCg                                          >  2:334730/1‑108 (MQ=255)
                                                                                                                                     ccctgccctTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCg                                          <  1:334730/108‑1 (MQ=255)
                                                                                                                                     ccctgcccgTAGGGCCGGGCGGTAACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCCTGCATTTGCCGCCCAGGCCCGTCCGCCCAGGAACAGCTTGTACCCCGGGTTGTGAGTCTCGTACCAGGACCCGGACCCCaacaagta  >  4:2289/1‑141 (MQ=255)
                                                                                                                                             |                                                                                                                                           
CCCCGAGCCGGATACCGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCGGTACCCCATCTCGTC  >  NC_002947/3902936‑3903215

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: ‑