Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,896,377 C→T 9.5% intergenic (+10/‑198) PP_4304 → / → sbp‑I VIC family cation transporter/sulfate ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,896,3770CT 9.5% 111.7 / 2.8 43intergenic (+10/‑198)PP_4304/sbp‑IVIC family cation transporter/sulfate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (18/20);  new base T (2/2);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TTTACCGCCGAGCTGGCCAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAC  >  NC_002947/4896233‑4896506
                                                                                                                                                |                                                                                                                                 
tttACCGCCGAGCTGGCCAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCG                                                                                                                               >  4:183366‑M1/1‑149 (MQ=255)
 ttACCGCCGAGCTGGCCAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACCAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAAGGGGTTGTTGACCAAAGCGTGATAGCGACTGTGCCGG                                                                                                                              >  8:318882‑M1/1‑149 (MQ=255)
            cTGGCCAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGTGTGATAGCGCCTGCGCCGGCCTCTTCGCGG                                                                                                                   <  1:214128‑M1/149‑1 (MQ=255)
             ttgccAATGCGATGCGTGGGGAGCAGCTGCAGCTTGATTGCCCGACCTGCCACAAGAAGACCCAGGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGATTGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGG                                                                                                                  <  5:92930‑M1/147‑1 (MQ=255)
                 cAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAA                                                                                                              >  4:153107‑M1/1‑149 (MQ=255)
                   aTGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATc                                                                                                            <  3:183366/149‑1 (MQ=255)
                                ggAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGgc                                                                                               >  1:236435/1‑149 (MQ=255)
                                 gAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGgcg                                                                                              >  6:223639/1‑149 (MQ=255)
                                              aTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGcgcg                                                                                 <  5:223639/149‑1 (MQ=255)
                                              aTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCCGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGcgcg                                                                                 <  5:76412/149‑1 (MQ=255)
                                                      ccGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGTGCCGGCCTCTTCGCGGGTAAAcccgctcctacgggggaccgcgcacggtccgcgcctg                                                                         >  8:66285‑M1/1‑112 (MQ=255)
                                                       cGACCTGCCACAAGAAGACCCCCGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAAGGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGt                                                                        <  3:14540/149‑1 (MQ=255)
                                                         aCCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAg                                                                      <  3:153107/149‑1 (MQ=255)
                                                              ccACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAcccgctcctacgggggaccgcgcacggtccgcgcctgtaggagcg                                                                 >  7:206633‑M1/1‑104 (MQ=255)
                                                                acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAtccgctcctacgggggaccgcgcacggtccgcgcctgtagg                                                                     <  3:240857‑M1/143‑42 (MQ=255)
                                                                acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAtccgctcctacgggggaccgcgcacggtccgcgcctgtagg                                                                     >  4:240857‑M1/1‑102 (MQ=255)
                                                                acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCggg                                                               >  7:290130/1‑149 (MQ=255)
                                                                acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGCGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCACTTAGCCGGTAAATCCGCACCTACACCCGCGGACCCAGCGCCGGCCCACCTAGCAGCcgg                                                               >  5:310614/1‑149 (MQ=255)
                                                                 caAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGTGCCGGCCTCTTCGCGGGTAAAcccgctcctacgggggaccgcgcacggtccgcgcctgtaggagcggat                                                              <  7:318882‑M1/149‑49 (MQ=255)
                                                                  aagaagACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGTGCCGGCCTCTTCGCGGGTAAAcccgctcctacgggggaccgcgcacggtccgcgcctgtaggagcggatt                                                             <  7:66285‑M1/149‑50 (MQ=255)
                                                                    gaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAcccgctcctacgggggaccgcgcacggtccgcgcctgtaggagcggattta                                                           <  8:206633‑M1/149‑52 (MQ=255)
                                                                    gaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATcc                                                                                                           >  5:363727/1‑101 (MQ=255)
                                                                    gaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATcc                                                                                                           <  6:363727/101‑1 (MQ=255)
                                                                          ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAGAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAcccgctcctacgggggaccgcgcacggtccgcgcctgtaggagcggat                                                              >  1:118247‑M1/1‑92 (MQ=255)
                                                                          ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAGAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAcccgctcctacgggggaccgcgcacggtccgcgcctgtaggagcggat                                                              <  2:118247‑M1/140‑49 (MQ=255)
                                                                           cccgcgCAGGGGGCGGCGTTCTGTTCTCGCTGTGGTAATGGGTTGTTTCCCAACGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGTTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTCGTAGCGGGTTTACCCGCGa                                                    <  8:290130/144‑1 (MQ=255)
                                                                                gCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGtttacccgcgaagagg                                               <  1:215820‑M2/149‑17 (MQ=255)
                                                                                            gttctgttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGtttacccgcgaagaggccg                                            <  6:256850‑M2/140‑20 (MQ=255)
                                                                                            gttctgttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGtttacccgcgaagaggccg                                            >  5:256850‑M2/1‑121 (MQ=255)
                                                                                             ttctgttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCCCGGGTAAATCCGCGCCTACAGGCGCGGACCGTGCGCGGTCCCCCGGAGGAGCGGGtttacccgcgaagaggccgtgacagg                                     <  8:157163‑M2/146‑27 (MQ=255)
                                                                                                tgttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAg                               >  6:264293/1‑149 (MQ=255)
                                                                                                  ttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAg                               <  2:215820/147‑1 (MQ=255)
                                                                                                       ggtgtgGGAATGGGTGTTTTCAGAAAGCGTTATAGCGCGTGGGGCGGCGTCTTCGCGGGTAATTGCGCTGGTACTGGCGGGGAGCGTGGGCGGTCCCCGGTTGGAGCGGGTTTACCGGCTATGAGGCCGGCACAGGCAGAAGCATAAc                         <  6:310614/146‑1 (MQ=255)
                                                                                                       gCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGAAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATaa                          <  7:157163/147‑1 (MQ=255)
                                                                                                             ggAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTTAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGGAGGAGCGGGTTTAACCGCGAAGCGGCCGGCACCGGCAGAAGCATAACCAAAGcc                  >  3:328453/1‑148 (MQ=255)
                                                                                                               aaTGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCATCGCGGGTTAATCCGCTCCTACAGGCGCGGCCCGTGGGCGGTCCCCCGTTGGAGCGGGTTTTCCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCAtt                >  1:122271/1‑149 (MQ=39)
                                                                                                                    gttgttTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGCGCGGGTTTACCCGCGCAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTg           >  5:260402/1‑149 (MQ=255)
                                                                                                                    gttgttTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCCCGGGTAAATCCGCTCCTACAGGCGCCGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCACAGGCAGAAGCATAACCAAAGCATTTTTTg           >  6:103088/1‑149 (MQ=255)
                                                                                                                     ttgttTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATACGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGt          >  8:166621/1‑149 (MQ=255)
                                                                                                                       aattCCCACAGCGTTATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCAACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGTTTTACCGGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTc        <  5:315431/147‑1 (MQ=38)
                                                                                                                             cAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAc  <  5:264293/149‑1 (MQ=255)
                                                                                                                              aaaGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTccccc                                                                          >  3:133932/1‑76 (MQ=255)
                                                                                                                              aaaGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTccccc                                                                          <  4:133932/76‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                 
TTTACCGCCGAGCTGGCCAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAC  >  NC_002947/4896233‑4896506

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: