Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,059,548 Δ1 bp 7.8% coding (1334/1581 nt) 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,5480C. 7.8% 153.7 / 7.9 51coding (1334/1581 nt)PP_4458opine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (24/23);  new base . (2/2);  total (26/25)
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.49e-01

CCCAGACCTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTCGCGT  >  NC_002947/5059402‑5059695
                                                                                                                                                  |                                                                                                                                                    
cccAGACCTGTTTCAGACCAAACAAAGACATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGa                                                                                                                                                   >  1:870725/1‑150 (MQ=255)
  cAGACCTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGaac                                                                                                                                                 <  1:360383/150‑2 (MQ=255)
   aGACCTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAAt                                                                                                                                                >  2:1101888/1‑150 (MQ=255)
       cTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGa                                                                                                                                           >  2:572900/1‑150 (MQ=255)
       cTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCGAATTTCgga                                                                                                                                          <  2:743643/150‑3 (MQ=255)
                aCCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCGAATTTCGGACATTGctca                                                                                                                                  >  1:604219/1‑147 (MQ=255)
                aCCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCGAATTTCGGACATTGctca                                                                                                                                  <  2:604219/150‑4 (MQ=255)
                      caaaGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTgc                                                                                                                            <  2:335559/150‑1 (MQ=255)
                       aaaGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTgcg                                                                                                                           >  1:993151/1‑150 (MQ=255)
                          gggATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCg                                                                                                                        <  1:659747/150‑1 (MQ=255)
                              tGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATgg                                                                                                                    <  2:22694/150‑1 (MQ=255)
                                    aTACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTc                                                                                                                                              >  1:440174/1‑119 (MQ=255)
                                    aTACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTc                                                                                                                                              <  2:440174/119‑1 (MQ=255)
                                           gtcgtGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATc                                                                                                       <  2:185288/150‑1 (MQ=255)
                                                             cgcCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTc                                                                                          >  2:873811/1‑145 (MQ=255)
                                                             cgcCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTc                                                                                          <  1:873811/145‑1 (MQ=255)
                                                                cTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGc                                                                                                                                  >  1:528711/1‑102 (MQ=255)
                                                                cTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGc                                                                                                                                  <  2:528711/102‑1 (MQ=255)
                                                                     tCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCgg                                                                             <  1:1103437/150‑1 (MQ=255)
                                                                        acaAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGgcgc                                                                          >  1:1110049/1‑150 (MQ=255)
                                                                        acaAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGATACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGt                                                                                           <  2:878372/133‑1 (MQ=255)
                                                                        acaAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGATACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGt                                                                                           >  1:878372/1‑133 (MQ=255)
                                                                          aaaTCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGccc                                                                        <  2:440586/150‑1 (MQ=255)
                                                                           aaTCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGcc                                                                                                            >  2:1027257/1‑113 (MQ=255)
                                                                           aaTCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGcc                                                                                                            <  1:1027257/113‑1 (MQ=255)
                                                                             tCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTgcgc                                                                                                                          >  2:287479/1‑97 (MQ=255)
                                                                             tCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTgcgc                                                                                                                          <  1:287479/97‑1 (MQ=255)
                                                                                       tCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCa                                                                                     >  2:700309/1‑124 (MQ=255)
                                                                                       tCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCa                                                                                     <  1:700309/124‑1 (MQ=255)
                                                                                         tGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAgcagca                                                         >  2:897402/1‑150 (MQ=255)
                                                                                                  gCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAac                                                >  2:610128/1‑150 (MQ=255)
                                                                                                  gCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGGGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAac                                                >  2:617178/1‑150 (MQ=255)
                                                                                                   ccTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAaca                                               >  2:277194/1‑150 (MQ=255)
                                                                                                     tCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACg                                             >  1:430467/1‑150 (MQ=255)
                                                                                                        aCCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGCCATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTATGTAAGAAAGGACTCGGGGTTCGGATATCGAGCAGCACCGTAAAACACGATc                                          <  1:1101888/150‑1 (MQ=255)
                                                                                                        aCCCTGGCGTTCAGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGACGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATc                                          >  1:294049/1‑150 (MQ=255)
                                                                                                          ccTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGACATTTGAAGTAGCAGCGGGTGCGCCGAGTGCAGAAGCAGAATCGTAAAACTCTGTGATCCGTGAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATcgc                                        >  1:570489/1‑150 (MQ=255)
                                                                                                          ccTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATcgc                                        >  1:147698/1‑150 (MQ=255)
                                                                                                            tGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCgcgc                                      >  2:1017053/1‑150 (MQ=255)
                                                                                                             ggCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACTTTGCAGCGGGTGCGCCGATGGCAGAAGCCGCATCGTAAAACTCTGTCAGCCAATAGTCGGCGCCCGGATATCGAGCAGCACCGCAAAACACGATCGCGct                                     <  2:570489/150‑1 (MQ=255)
                                                                                                                     ggATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCtt                             >  1:1118567/1‑150 (MQ=255)
                                                                                                                            atcatcAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCt                      <  1:610128/150‑1 (MQ=255)
                                                                                                                                tcaGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAActct                                                                                             <  1:241576/75‑1 (MQ=255)
                                                                                                                                tcaGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAActct                                                                                             >  2:241576/1‑75 (MQ=255)
                                                                                                                                tcaGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTg                  <  2:506782/150‑1 (MQ=255)
                                                                                                                                   gCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGcc               <  2:430467/150‑1 (MQ=255)
                                                                                                                                      aCTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTgg            >  1:385757/1‑150 (MQ=255)
                                                                                                                                        gctgCAATGTCGAATTTCGGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCCTGACGGCTTTTGGCCTGGTa          >  2:885871/4‑150 (MQ=255)
                                                                                                                                              aTGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTcgc    <  2:385757/150‑1 (MQ=255)
                                                                                                                                               tGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTcgcg   <  1:617178/150‑1 (MQ=255)
                                                                                                                                                gTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTCGCGt  <  2:147698/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCCAGACCTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTCGCGT  >  NC_002947/5059402‑5059695

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