Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,059,872 C→G 8.8% G321A (GGA→GCA)  PP_4458 ← opine ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,059,8720CG 8.8% 157.0 / 5.5 55G321A (GGA→GCA) PP_4458opine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (29/21);  new base G (2/3);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

CGTTGGAAACAACCGATTTGAGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAAGAGAGTACGGAATTCGGCAGGC  >  NC_002947/5059732‑5060009
                                                                                                                                            |                                                                                                                                            
cgTTGGAACCAACCGATTTGAGTTTTAACCCGCTGGGATAACCCGCGCTGGCCAACAAGGCGTTCGCCTGCGCAGGATCGTACGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGCGGATGC‑‑‑TCCATCa                                                                                                                                   <  4:121681/149‑1 (MQ=255)
 gTTGGAAACAACCGATTTGAGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAg                                                                                                                                  >  2:169517/1‑149 (MQ=255)
    ggAAACAACCGATTTGAGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTa                                                                                                                               >  2:122123/1‑149 (MQ=255)
     gAAACAACCGATTTGAGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTa                                                                                                                               >  1:61299/1‑148 (MQ=255)
     gAAACAACCGATTTGAGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTa                                                                                                                               >  7:4602/1‑148 (MQ=255)
           acCGATTTGAGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCttt                                                                                                                        <  4:148715/149‑1 (MQ=255)
            ccGATTTGAGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGCGGATGC‑‑‑TCCATCatac                                                                                                                                >  4:263210/1‑137 (MQ=255)
            ccGATTTGAGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGCGGATGC‑‑‑TCCATCatac                                                                                                                                <  3:263210/140‑4 (MQ=255)
               aTTTGAGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTAccc                                                                                                                    <  4:99540/149‑1 (MQ=255)
                    aGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTa                                                                                                               <  1:255470/149‑1 (MQ=255)
                           gcccGTTGGGATAAACCGCGCTGGACAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATAGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCACCGAACCGAAt                                                                                                         >  2:306772/2‑148 (MQ=255)
                           aCCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAAt                                                                                                         >  7:18440/1‑148 (MQ=255)
                            cccGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAg                                                                                                              <  8:117119/142‑1 (MQ=255)
                            cccGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAg                                                                                                              >  7:117119/1‑142 (MQ=255)
                              cGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGAtc                                                                                                     <  8:4602/149‑1 (MQ=255)
                                        aCCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATg                                                                                            <  4:226937/148‑1 (MQ=255)
                                         cccGCGCTGGCCAAGAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCATTCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGtt                                                                                          <  5:288624/149‑1 (MQ=255)
                                         cccGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGCGGATGC‑‑‑TCCATCatac                                                                                                                                >  6:266436/1‑108 (MQ=255)
                                         cccGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGCGGATGC‑‑‑TCCATCatac                                                                                                                                <  5:266436/111‑4 (MQ=255)
                                          ccGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGtt                                                                                          >  3:290778/1‑148 (MQ=255)
                                            gcgcTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGAc                                                                                       <  7:10208/149‑1 (MQ=255)
                                                             tttCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCtt                                                                      >  7:229886/1‑149 (MQ=255)
                                                                  cgTGCGCAGGATCGTTAGCGTATTGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCCCTTTCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGc                                                                    <  1:122123/149‑1 (MQ=255)
                                                                       gcAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGAGCTCATTGATGGTGACCGACGCAGCGATAGCGTGCCGAACCtt                                                            >  4:305916/1‑149 (MQ=255)
                                                                          tgATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGAc                                                         <  1:367298/148‑1 (MQ=255)
                                                                          ggATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGAc                                                         <  2:61299/149‑1 (MQ=255)
                                                                          ggATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGAc                                                         <  1:361565/149‑1 (MQ=255)
                                                                           gATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGATACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGAcg                                                        >  5:207816/1‑149 (MQ=255)
                                                                             tCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGAcgtc                                                      >  1:197160/1‑149 (MQ=255)
                                                                                tAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCATAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACgtcgtc                                                   >  4:383046/1‑149 (MQ=255)
                                                                                       aTGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGCCAGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTACCGGc                                            >  2:248243/1‑149 (MQ=255)
                                                                                          gCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCtt                                          <  2:197160/148‑1 (MQ=255)
                                                                                           cccGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACGTCTTTACCCGCGTAGCGAATGATCTCATTGATGTTTACCGACGCAGCGATAGCTTGc                                                                    <  2:200801/121‑1 (MQ=255)
                                                                                           cccGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACGTCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGc                                                                    >  1:200801/1‑121 (MQ=255)
                                                                                            ccGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGa                                        >  7:255338/1‑148 (MQ=255)
                                                                                            ccGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGa                                        >  2:282927/1‑148 (MQ=255)
                                                                                            ccGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGa                                        >  5:77345/1‑148 (MQ=255)
                                                                                                  cTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCg                                 <  1:185828/149‑1 (MQ=255)
                                                                                                    gCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGt                               <  8:229886/149‑1 (MQ=255)
                                                                                                         cAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTaag                           >  6:356054/1‑148 (MQ=255)
                                                                                                             ccccTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTaagaaaga                      <  8:18440/149‑1 (MQ=255)
                                                                                                               ccTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTaagaaaga                      >  2:45507/1‑147 (MQ=255)
                                                                                                               ccTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTaagaaaga                      >  7:16083/1‑147 (MQ=255)
                                                                                                               ccTGGTAGCCATTGGGTATGATGGAGAATCC‑‑‑GCCATCAGCTACATCTTTACCCTCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAagaga                    <  3:333568/149‑1 (MQ=255)
                                                                                                                cTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAagagag                   >  6:11206/1‑149 (MQ=255)
                                                                                                                cTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACAGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTaagaaaga                      >  5:33844/1‑146 (MQ=255)
                                                                                                                    tAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAAGAGAGTACg               >  5:123722/1‑149 (MQ=255)
                                                                                                                     aGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGATTGATATCATTGATGTTGACCGACGCAGAGATAGCTTTCCGAACCTTGACGTCGTATAACGTCTTGATATTCAGGTTAAGAAATAGATTAct               >  5:352134/1‑147 (MQ=255)
                                                                                                                             ggTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAAGAGGGTACGGAATTCgg       >  7:91462/1‑148 (MQ=255)
                                                                                                                             ggTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAAGAGAGTACGGAATTCgg       >  3:2924/1‑148 (MQ=255)
                                                                                                                              gTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAAGAGAGTACGGAATTCGGCa     <  5:11206/149‑1 (MQ=255)
                                                                                                                               tATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAAGAGAGTACGGAATTCGGCAg    <  7:265067/149‑1 (MQ=255)
                                                                                                                               tATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAAGAGAGTACGGAATTCGGCAg    >  2:176569/1‑149 (MQ=255)
                                                                                                                               tATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAAGAGAGTACGGAATTCGGCAg    >  5:198080/1‑149 (MQ=255)
                                                                                                                                 tgatgGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAAGAGAGTACGGAATTCGGCAGGc  <  8:16083/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                            
CGTTGGAAACAACCGATTTGAGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAAGAGAGTACGGAATTCGGCAGGC  >  NC_002947/5059732‑5060009

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