Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,948,683 T→A 8.5% intergenic (‑64/+18) PP_mr62 ← / ← trxA Pseudomon‑Rho regulatory region/thioredoxin I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,948,6830TA 8.5% 137.7 / 3.4 47intergenic (‑64/+18)PP_mr62/trxAPseudomon‑Rho regulatory region/thioredoxin I
Reads supporting (aligned to +/- strand):  ref base T (20/23);  new base A (2/2);  total (22/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.67e-01

GATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCGGG  >  NC_002947/5948541‑5948827
                                                                                                                                               |                                                                                                                                                
gATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCggg                                                                                                                                            <  5:235724/149‑1 (MQ=255)
 aTGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCgggg                                                                                                                                           >  7:146001/1‑149 (MQ=255)
       aCAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCtttt                                                                                                                                     <  4:245877/149‑1 (MQ=255)
        cAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGTCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCtttt                                                                                                                                     >  1:58490/1‑148 (MQ=255)
        cAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTCTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCtttt                                                                                                                                     <  1:277848/149‑1 (MQ=255)
           gATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTttaca                                                                                                                                 >  2:344956/1‑149 (MQ=255)
           gATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCAAATGCGGGGCttttctttt                                                                                                                                 <  6:223366/149‑6 (MQ=255)
            aTCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTttacat                                                                                                                                <  1:344956/149‑1 (MQ=255)
                    cTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAGGCCCCGCAAATGCGGGGCttttct                                                                                                                                    >  5:31728/1‑135 (MQ=255)
                    cTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAGGCCCCGCAAATGCGGGGCttttct                                                                                                                                    <  6:31728/137‑3 (MQ=255)
                       aTTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTgg                                                                                                                       <  6:354214/147‑1 (MQ=255)
                                  ggATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATg                                                                                                          <  8:146001/149‑1 (MQ=255)
                                   gaggcgtCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTCCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGc                                                                                                         <  8:83648/143‑1 (MQ=255)
                                    aTGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCg                                                                                                        <  3:229499/149‑1 (MQ=255)
                                    aTGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCg                                                                                                        <  6:94865/149‑1 (MQ=255)
                                                     tGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGa                                                                                       <  6:54510/149‑1 (MQ=255)
                                                               aCTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGcc                                                                              >  8:215442/1‑148 (MQ=255)
                                                                   gAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGTGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCt                                                                         >  5:72484/1‑149 (MQ=255)
                                                                     aaaGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtg                                                                       >  7:319391/1‑149 (MQ=255)
                                                                       aGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtggt                                                                     <  8:319391/149‑1 (MQ=255)
                                                                                  gCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGt                                                          >  4:181543/1‑149 (MQ=255)
                                                                                    cGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTg                                                        <  7:215442/149‑1 (MQ=255)
                                                                                    cGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTg                                                        <  2:58490/149‑1 (MQ=255)
                                                                                    cGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTg                                                        <  3:181543/149‑1 (MQ=255)
                                                                                         gTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTc                                                               >  3:77404/1‑137 (MQ=255)
                                                                                         gTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTc                                                               <  4:77404/137‑1 (MQ=255)
                                                                                             cacTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtg                                                                       >  3:85268/1‑125 (MQ=255)
                                                                                             cacTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtg                                                                       <  4:85268/125‑1 (MQ=255)
                                                                                                       ttttCTGCGTCTAGTGCTCAGAAAAAGAAAAGCCCCGCATTTGCCGGGCCTTTTCACATCACAGCTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCAAACCTTGGTGGCATCGACGTAGCCGTACTTGAACAGCAt                                       >  4:91831/1‑148 (MQ=255)
                                                                                                           cTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGc                                   <  6:292641/148‑1 (MQ=255)
                                                                                                                   aGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATa                          >  8:114533/1‑149 (MQ=255)
                                                                                                                    gTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATa                          >  3:272884/1‑148 (MQ=255)
                                                                                                                     tGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCtt                                                <  4:93151/125‑1 (MQ=255)
                                                                                                                     tGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCtt                                                >  3:93151/1‑125 (MQ=255)
                                                                                                                       ctctGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCa                       <  7:332229/148‑1 (MQ=255)
                                                                                                                           gAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGca                   >  5:6673/1‑148 (MQ=255)
                                                                                                                            aaaaaGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGGTTCGACGTTGCCGTTCTTGAACAGc                                         >  8:226160/1‑125 (MQ=255)
                                                                                                                            aaaaaGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGc                                         <  7:226160/125‑1 (MQ=255)
                                                                                                                             aaaagaaaagCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGaa                                             <  2:87269/120‑1 (MQ=255)
                                                                                                                             aaaagaaaagCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGaa                                             >  1:87269/1‑120 (MQ=255)
                                                                                                                               aagaaaagACCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGAGTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACACAGCACAcct              >  8:155960/1‑148 (MQ=255)
                                                                                                                                 aaaaggCCCCGCAAATGCGGGGCTTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTg            >  3:208922/2‑149 (MQ=255)
                                                                                                                                  aaaagCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGc           >  7:113760/1‑149 (MQ=255)
                                                                                                                                         ccGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCg    >  8:95484/1‑149 (MQ=255)
                                                                                                                                           gCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACa                                           >  3:8535/1‑108 (MQ=255)
                                                                                                                                           gCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACa                                           <  4:8535/108‑1 (MQ=255)
                                                                                                                                           gCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCggg  <  8:113760/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCGGG  >  NC_002947/5948541‑5948827

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: