Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,061,134 G→A 8.5% intergenic (‑301/‑190) PP_4458 ← / → PP_4459 opine ABC transporter substrate‑binding protein/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,061,1340GA 8.5% 127.6 / 3.0 47intergenic (‑301/‑190)PP_4458/PP_4459opine ABC transporter substrate‑binding protein/transposase
Reads supporting (aligned to +/- strand):  ref base G (19/24);  new base A (2/2);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

CTTTTCCATAGCTTCAAAAACATCGGTCAATGTCTTTATCCTATAGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTGGTAGGTCTCCATCTCACATTGAGGCCGCACTGACTCTGCTAGTGTTCTGTCTCAAAAATAAGTATTCCAT  >  NC_002947/5061001‑5061280
                                                                                                                                     |                                                                                                                                                  
cttTTCCATAGCTTCAAAAACATCGGTCAATGTCTTTATCCTATAGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAAc                                                                                                                                      <  4:63943/148‑1 (MQ=255)
cttTTCCATAGCTTCAAAAACATCGGTCAATGTCTTTATCCTATAGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACa                                                                                                                                     >  3:361395/1‑149 (MQ=255)
 ttttCCATAGCTTCAAAAACATCGGTCAATGTCTTTATCCTATAGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACAt                                                                                                                                    >  3:298755/1‑149 (MQ=255)
              cAAAAACATCGGTCAATGTCTTTATCCTATAGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTaca                                                                                                                        <  4:361395/148‑1 (MQ=255)
               aaaaaCATCGGTCAATGTCTTTATCCTATAGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTacac                                                                                                                       >  3:113173/1‑148 (MQ=255)
                     aTCGGTCAATGTCTTTATCCTATAGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAgagc                                                                                                                <  1:350534/149‑1 (MQ=255)
                         gTCAATGTCTTTATCCTATAGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAgagcgagc                                                                                                            <  4:231054/149‑1 (MQ=255)
                          tCAATGTCTTTATCCTATAGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGagctagcg                                                                                                           >  3:84937/1‑149 (MQ=255)
                                 cTTTATCCTATAGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTc                                                                                                                                 >  6:191324/1‑120 (MQ=255)
                                 cTTTATCCTATAGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTc                                                                                                                                 <  5:191324/120‑1 (MQ=255)
                                    tATCCTATAGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCACTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGa                                                                                                 >  7:134550/1‑149 (MQ=255)
                                    tATCCTATAGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCACGGCTATGGAACGCTCCAAAATGAAAAATCGTAGTTACGACCAATACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGa                                                                                                 >  7:129414/1‑149 (MQ=255)
                                     aTCCTATAGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTAGTGAGTACACTAAGAGAGAGCGTGAATACGGAt                                                                                                >  8:318454/1‑149 (MQ=255)
                                      tCCTATAGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGAt                                                                                                <  4:46880/148‑1 (MQ=255)
                                      tCCTATAGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGAt                                                                                                <  4:113173/148‑1 (MQ=255)
                                         tataGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATgaga                                                                                            <  6:244384/149‑1 (MQ=255)
                                         tataGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATgaga                                                                                            >  1:254908/1‑149 (MQ=255)
                                         tataGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCACGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATgaga                                                                                            >  8:217353/1‑149 (MQ=255)
                                          ataGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAgagcgag                                                                                                             >  2:875/1‑131 (MQ=255)
                                          ataGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAgagcgag                                                                                                             <  1:875/131‑1 (MQ=255)
                                           taGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATgaga                                                                                            <  1:94870/147‑1 (MQ=255)
                                                    aGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCgtgt                                                                                 <  5:80754/149‑1 (MQ=255)
                                                      aaTATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTc                                                                                >  1:202803/1‑148 (MQ=255)
                                                             cTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTg                                                                        >  6:324069/1‑149 (MQ=255)
                                                                ttCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTgg                                                                       <  2:202803/147‑1 (MQ=255)
                                                                  cGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTGGTAg                                                                    <  2:138033/148‑1 (MQ=255)
                                                                             cTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTGGTAGGTCTCCAtct                                                          <  1:52978/147‑1 (MQ=255)
                                                                             cTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTGGTAGGTCTCCATCTca                                                        <  4:370499/149‑1 (MQ=255)
                                                                                tAGTGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTTCGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTGGTAGGTCTCCATCTcac                                                       <  8:285040/147‑1 (MQ=255)
                                                                                          ccGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCGCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTGGTAGGTCTCCATCTCACATTGAGGCCGCa                                           <  5:324069/149‑1 (MQ=255)
                                                                                              cTATGGAACGCTCCAAAATGAAACCTCGTAGTTACGCCCACTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTGGTAGGTCTCCATCTCACATTGAGGCCGCACTGa                                       <  8:134550/149‑1 (MQ=255)
                                                                                              cTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTGGTAGGTCTCCATCTCACATTGAGGCCGCACTGa                                       <  2:254908/149‑1 (MQ=255)
                                                                                                aTGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAAt                                                                             >  2:368370/1‑109 (MQ=255)
                                                                                                aTGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAAt                                                                             <  1:368370/109‑1 (MQ=255)
                                                                                                 tGGAACGCTCCAAAATGAAACATCTTAGTTACTCCCACTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTGGTAGGTCTCCATCTCACATTGAGGCCGCACTGa                                       <  8:129414/146‑1 (MQ=255)
                                                                                                   gAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTGGTAGGTCTCCATCTCACATTGAGGCCGCACTGActct                                   <  3:233367/148‑1 (MQ=255)
                                                                                                   gAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTGGTAGGTCTCCATCTCACATTGAGGCCGCACTGACTCTg                                  >  3:162489/1‑149 (MQ=255)
                                                                                                         tCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTGGTAGGTCTCCATCTCACATTGAGGCCGCACTGACTCTGCTAgtg                            >  6:32069/1‑149 (MQ=255)
                                                                                                           cAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTGGTAGGTCTCCATCTCACATTg                                                   >  7:184871/1‑124 (MQ=255)
                                                                                                           cAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTGGTAGGTCTCCATCTCACATTg                                                   <  8:184871/124‑1 (MQ=255)
                                                                                                               aTGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTGGTAGGTCTCCATCTCACATTGAGGCCGCACTGACTCTGCTAGTGTTCTGt                      >  7:295129/1‑149 (MQ=255)
                                                                                                                  aaaCATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTGGTAGGTCTCCATCTCACATTGAGGCCGCACTGACTCTGCTAGTGTTCTGtct                    >  1:296971/1‑148 (MQ=255)
                                                                                                                       tCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTGGTAGGTCTCCATCTCACATTGAGGCCGCACTGACTCTGCTAGTGTTCTGTCTCaaaaa              >  5:275990/1‑149 (MQ=255)
                                                                                                                         gTAGTTACGCCAGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTGGTAGGTCTCCATCTCACATTGAGGCCGCACTGACTCTGCTAGTGTTCTGTCTCAAAAAt             <  7:318454/148‑1 (MQ=255)
                                                                                                                            gTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTGGTAGGTCTCCATCTCACATTGAGGCCGCACTGACTCTGCTAGTGTTCTGTCTCAAAAATAAg          <  4:84937/148‑1 (MQ=255)
                                                                                                                                   ccGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTGGTAGGTCTCCATCTCACATTGAGGCCGCACTGACTCTGCTAGTGTTCTGTCTCAAAAATAAGTATTCCAt  <  7:217353/149‑1 (MQ=255)
                                                                                                                                   ccGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTGGTAGGTCTCCATCTCACATTGAGGCCGCACTGACTCTGCTAGTGTTCTGTCTCAAAAATAAGTATTCCAt  >  3:203755/1‑149 (MQ=255)
                                                                                                                                     |                                                                                                                                                  
CTTTTCCATAGCTTCAAAAACATCGGTCAATGTCTTTATCCTATAGAAAATAAGAATATTTCTTTTCGAGAGCGCTACTGTAGGGCGTGCCCGGCTATGGAACGCTCCAAAATGAAACATCGTAGTTACGCCCGCTACCTTCAGTAACATTTCGTGCGTACACTAAGAGCGAGCGTGAATACGGATGAGACGCAGTCGTGTCAATAGCTGGTAGGTCTCCATCTCACATTGAGGCCGCACTGACTCTGCTAGTGTTCTGTCTCAAAAATAAGTATTCCAT  >  NC_002947/5061001‑5061280

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