Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,521,682 A→G 6.5% T49T (ACA→ACG U712_18050 → Ribose import ATP‑binding protein RbsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,521,6820AG 6.5% 189.0 / 2.3 61T49T (ACA→ACGU712_18050Ribose import ATP‑binding protein RbsA
Reads supporting (aligned to +/- strand):  ref base A (28/29);  new base G (2/2);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATGAAAGACATTCATAAAACATTCGGAAAAAATCAGGTGCTGTCAGGCGTTTCCTTTCAGCTCATGCCTGGCGAGGTTCACGCATTAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGAATAGCCTTTATCCATCAGGAATTGAATATCTGGCCGGAAATGACCGTTC  >  CP006881/3521546‑3521818
                                                                                                                                        |                                                                                                                                        
aaatgaaaGACATTCATAAAACATTCGGAAAAAATCAGGTGCTGTCAGGCGTTTCCTTTCAGCTCATGCCTGGCGAGGTTCACGCATTAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAgg                                                                                                                                        >  8:261812/1‑139 (MQ=255)
 aatgaaaGACATTCATAAAACATTCGGAAAAAATCAGGTGCTGTCAGGCGTTTCCTTTCAGCTCATGCCTGGCGAGGTTCACGCATTAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGc                                                                                                                                       <  2:239420/139‑1 (MQ=255)
                  aaaCATTCGGTAAAAATCCGGTGCGGTCAGGCGTTTCCTTTCTGCTCATGCGTGGCGAGGTTCACGCATTAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATa                                                                                                                       <  2:15000/138‑1 (MQ=255)
                       ttCGGAAAAAATCAGGTGCTGTCAGGCGTTTCCTTTCAGCTCATGCCTGGCGAGGTTCACGCATTAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGt                                                                                                                  >  1:188481/1‑138 (MQ=255)
                        tCGGAAAAAATCAGGTGCTGTCAGGCGTTTCCTTTCAGCTCATGCCTGGCGAGGTTCACGCATTAATGGGAGAAACCGGCGCCGGCAAGTCCACGCTTATGACCATTTTGACAGGCCTGCACAAAGCAGATAAAGGTc                                                                                                                 >  7:88696/1‑138 (MQ=255)
                         cGGAAAAAATCAGGTGCTGTCAGGCGTTTCCTTTCAGCTCATGCCTGGCGAGGTTCACGCATTAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCa                                                                                                                >  4:329000/1‑138 (MQ=255)
                                                   ttgccTTTCAGCTCATGCCTGGCGAGGGTCACGCATAAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTa                                                                                     <  4:293986/136‑1 (MQ=255)
                                                      ccTTTCAGCTCATGCCTGGCGAGGTTCACGCATTAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACGGGc                                                                                                                                       >  8:11567/1‑86 (MQ=255)
                                                      ccTTTCAGCTCATGCCTGGCGAGGTTCACGCATTAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACGGGc                                                                                                                                       <  7:11567/86‑1 (MQ=255)
                                                          tCAGCTCATGCCTGGCGAGGTTCACGCATTAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTc                                                                               >  8:189538/1‑138 (MQ=255)
                                                              cTCATGCCTGGAGAGGTTCACGCATTAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACGGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATc                                                                          >  6:270863/1‑139 (MQ=255)
                                                                 aTGCCTGGCGAGGTTCACGCATTAATGGGAGAAAACGGCGCCGGCAAGTCCCCGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCg                                                                        >  5:284010/1‑138 (MQ=255)
                                                                 aTGCCTGGCGAGGTTCACGCATTAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCg                                                                        >  7:183612/1‑138 (MQ=255)
                                                                  tGCCTGGCGAGGTTCACGCATTAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGaa                                                                      <  8:183612/139‑1 (MQ=255)
                                                                       ggCGAGGTTCACGCATTAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAaagaag                                                                 >  3:370200/1‑139 (MQ=255)
                                                                            ggTTCACGCATTAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGaa                                                            <  5:365624/139‑1 (MQ=255)
                                                                            ggTTAACGCATTAATGGGACAAACCGGCGCCTGCAAGTCCACGCTTATGAACATTTTTACAGGCCTGCACAAAGCATATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGaa                                                            <  7:318179/139‑1 (MQ=255)
                                                                             gTTCACGCATTAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGaa                                                            <  3:329000/138‑1 (MQ=255)
                                                                                cACGCATTAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGc                                                        <  4:313221/139‑1 (MQ=255)
                                                                                 aCGCATTAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCa                                                       <  7:189538/139‑1 (MQ=255)
                                                                                 aCGCATTAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCa                                                       <  5:187570/139‑1 (MQ=255)
                                                                                  cGCATTAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGaa                                                            >  6:116011/1‑133 (MQ=255)
                                                                                  cGCATTAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGaa                                                            <  5:116011/133‑1 (MQ=255)
                                                                                      ttAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGaa                                                  <  3:325184/139‑1 (MQ=255)
                                                                                       tAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGaa                                                  >  5:61241/1‑138 (MQ=255)
                                                                                        aaTGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGAAt                                                 <  7:261812/138‑1 (MQ=255)
                                                                                               gacaccGGCGCCGGCAAGTCCACGCTTATGAACATTATGACACGCCTGCACCAAGCAGATAAAGGTCAAATCAGCATCAACGGAAACGACACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGAATAGCCTTTa                                         <  2:90230/134‑1 (MQ=255)
                                                                                               gaAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGaa                                                  >  4:282116/1‑130 (MQ=255)
                                                                                               gaAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGaa                                                  <  3:282116/130‑1 (MQ=255)
                                                                                               gaAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGAATAGCCTTTa                                         >  2:97841/1‑139 (MQ=255)
                                                                                               gaAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGAATAGCCTTTa                                         >  4:15581/1‑139 (MQ=255)
                                                                                                    cgccgcCGGCAAGGCCACGCGTATGAACATTTTGACAGGCCGGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGAATAGCCTTTATCCAt                                    <  5:300119/136‑1 (MQ=255)
                                                                                                    cGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGAATAGCCTTTATCCAt                                    <  3:101062/139‑1 (MQ=255)
                                                                                                       cgcCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAg                                                                                                                             <  4:189003/47‑1 (MQ=255)
                                                                                                       cgcCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAg                                                                                                                             >  3:189003/1‑47 (MQ=255)
                                                                                                        gcCGGCAAGTCCACGCTTATGAACATTTTGACGGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGAATAGCCTTTATCCATCAgg                                <  5:270863/139‑1 (MQ=255)
                                                                                                        gcCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGAAAAAGGTCAAATCAGCATACACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGAATAGCCTTTATCCATCAgg                                <  2:263235/139‑1 (MQ=255)
                                                                                                            gCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATATTTCCAATCCGAAAGAAGCGGAAc                                                           >  4:263542/1‑108 (MQ=255)
                                                                                                            gCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATATTTCCAATCCGAAAGAAGCGGAAc                                                           <  3:263542/108‑1 (MQ=255)
                                                                                                                   cACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGATACGTATTTTTCCAATCCTTACGCAGCGGAACAGCATGGATTAGCCTTTATCCATCAGGAATTGAatat                      <  7:41986/138‑1 (MQ=255)
                                                                                                                      gCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGa                                                             <  5:86029/96‑1 (MQ=255)
                                                                                                                      gCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGa                                                             >  6:86029/1‑96 (MQ=255)
                                                                                                                      gCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGAATAGCCTTTATCCATCAGGAATTGAATATCTgg                  <  8:111785/139‑1 (MQ=255)
                                                                                                                           tGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACg                                                                                                  <  4:83510/54‑1 (MQ=255)
                                                                                                                           tGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACg                                                                                                  >  3:83510/1‑54 (MQ=255)
                                                                                                                           tGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGAATAGCCTTTATCCATCAGGAATTGAATATCTGGCCgg              >  4:310231/1‑138 (MQ=255)
                                                                                                                           tGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGCAATAGCCTGTATCCATCAGGAAGTGAATATCTCGCCgg              >  8:120635/1‑138 (MQ=255)
                                                                                                                                  tttGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGgaaacgaaa                                                                                         <  4:123175/56‑1 (MQ=255)
                                                                                                                                  tttGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGgaaacgaaa                                                                                         >  3:123175/1‑56 (MQ=255)
                                                                                                                                  tttGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGgaaacg                                                                                            <  5:123156/53‑1 (MQ=255)
                                                                                                                                  tttGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGgaaacg                                                                                            >  6:123156/1‑53 (MQ=255)
                                                                                                                                  tttGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGaa                                                            <  4:1531/85‑1 (MQ=255)
                                                                                                                                  tttGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGaa                                                            >  3:1531/1‑85 (MQ=255)
                                                                                                                                  tttGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGAATAGCCTTTATcc                                      >  7:176211/1‑107 (MQ=255)
                                                                                                                                  tttGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGAATAGCCTTTATcc                                      <  8:176211/107‑1 (MQ=255)
                                                                                                                                  tttGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGAATAGCCTTTATCCATCAGGAATTGAatat                      >  8:105646/1‑123 (MQ=255)
                                                                                                                                  tttGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGAATAGCCTTTATCCATCAGGAATTGAatat                      <  7:105646/123‑1 (MQ=255)
                                                                                                                                  tttGACACGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGCAATAGCCTTTATCCATCAGCAATAGAATATCTGGCCGGAAATGAcc      >  3:172376/1‑139 (MQ=255)
                                                                                                                                    tGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGAATAGCCTTTATCCATCAGGAATTGAATATCTGGCCGGAAATGACCGt    >  1:232463/1‑139 (MQ=255)
                                                                                                                                      aCAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGAATAGCCTTTATCCATCAGGAATTGAATATCTGGCCGGAAATGACCGTTc  >  8:217344/1‑139 (MQ=255)
                                                                                                                                      aCAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGAATAGCCTTTATCCATCAGGAATTGAATATCTGGCCGGAAATGACCGTTc  >  8:222686/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                        
AAATGAAAGACATTCATAAAACATTCGGAAAAAATCAGGTGCTGTCAGGCGTTTCCTTTCAGCTCATGCCTGGCGAGGTTCACGCATTAATGGGAGAAAACGGCGCCGGCAAGTCCACGCTTATGAACATTTTGACAGGCCTGCACAAAGCAGATAAAGGTCAAATCAGCATAAACGGAAACGAAACGTATTTTTCCAATCCGAAAGAAGCGGAACAGCATGGAATAGCCTTTATCCATCAGGAATTGAATATCTGGCCGGAAATGACCGTTC  >  CP006881/3521546‑3521818

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