Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 780,559 A→C 9.5% intergenic (‑86/+36) PP_0670 ← / ← glyA‑II bile acid/Na+ symporter family transporter/serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947780,5590AC 9.5% 106.7 / 3.8 42intergenic (‑86/+36)PP_0670/glyA‑IIbile acid/Na+ symporter family transporter/serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (21/17);  new base C (2/2);  total (23/19)
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

GGCCAGTACGGCGGCGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCAGCCGGCCA  >  NC_002947/780414‑780707
                                                                                                                                                 |                                                                                                                                                    
ggCCAGTACGGCGGCGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTgg                                                                                                                                                  >  2:1225398/1‑150 (MQ=255)
ggCCAGTACGGCGGCGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTgg                                                                                                                                                  <  1:387999/150‑1 (MQ=255)
    aGTACGGCGGCGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGAATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAGAGAGCCCGGCACTGGGCCg                                                                                                                                              >  1:370416/1‑150 (MQ=255)
          gcggcgAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGCAGGAAACGAAAGAGCCCGGCACTGGGCCGGGctct                                                                                                                                        <  1:622373/150‑1 (MQ=255)
           cggcgAAGTAGATGGGGATCTGGCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATCCGAGGCGCAGTCTTTGTACGAGGGGCGTTGCGCCTGTATTAAACGATAGAGCCCGGCACTGGGCCGGGCTCtt                                                                                                                                       <  1:1053679/150‑1 (MQ=255)
                 aGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTCTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACtt                                                                                                                                 <  1:1225398/150‑1 (MQ=255)
                          ggATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAgc                                                                                                                        >  1:838200/1‑150 (MQ=255)
                                gCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAAt                                                                                                                  <  1:1062215/150‑1 (MQ=255)
                                            gTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTc                                                                                                                                     >  2:495052/1‑119 (MQ=255)
                                            gTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAACAGCCCGGCACTGGGCCGGGCTCTTTc                                                                                                                                     <  1:495052/119‑1 (MQ=255)
                                                  ctctCGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGc                                                                                                <  1:149899/150‑1 (MQ=255)
                                                   tctGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCa                                                                                               >  2:340125/1‑150 (MQ=255)
                                                       gTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGctt                                                                                           >  2:1190384/1‑150 (MQ=255)
                                                        tCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGcttc                                                                                          <  2:838200/150‑1 (MQ=255)
                                                            ggTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGc                                                                                      >  2:201285/1‑150 (MQ=255)
                                                                  gCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATgg                                                                                >  2:255508/1‑150 (MQ=255)
                                                                   ccTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGc                                                                               >  2:266858/1‑150 (MQ=255)
                                                                    cTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGcc                                                                              <  1:340125/150‑1 (MQ=255)
                                                                           aaGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCt                                                                       >  1:938663/1‑150 (MQ=255)
                                                                                  ggTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCg                                                                                                                                              <  2:724398/72‑1 (MQ=255)
                                                                                  ggTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCg                                                                                                                                              >  1:724398/1‑72 (MQ=255)
                                                                                      ttGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCacggacg                                                            >  2:580559/1‑150 (MQ=255)
                                                                                      ttGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTGAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCacggacg                                                            <  2:370416/150‑1 (MQ=255)
                                                                                          gggCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGt                                                        >  2:785495/1‑150 (MQ=255)
                                                                                                 gTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCAc                                                 >  1:206636/1‑150 (MQ=255)
                                                                                                           gAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCgttg                                       >  2:1091826/1‑150 (MQ=255)
                                                                                                                gCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTCCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAg                                  <  1:266858/150‑1 (MQ=255)
                                                                                                                       gcCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGcc                           >  1:1268026/1‑150 (MQ=255)
                                                                                                                        ccTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGTAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCa                          >  1:176289/1‑150 (MQ=255)
                                                                                                                              ttAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGt                                                        <  2:879069/114‑1 (MQ=255)
                                                                                                                              ttAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGt                                                        >  1:879069/1‑114 (MQ=255)
                                                                                                                                aaaCGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCANGGACGGNGTCGATCACCGCNTCGTTGTTCAGGTCGGCCAGGATGTCg                  <  1:1091826/150‑1 (MQ=255)
                                                                                                                                  gtgAAAGAGCCCGGCCCAGTGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGa                                                                          >  2:331088/3‑92 (MQ=37)
                                                                                                                                  gtgAAAGAGCCCGGCCCAGTGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGa                                                                          <  1:331088/90‑1 (MQ=37)
                                                                                                                                    gAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTgcc                                                                         <  1:185582/91‑3 (MQ=255)
                                                                                                                                    gAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGAc                                                                         >  2:185582/1‑91 (MQ=255)
                                                                                                                                      aaGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATc            <  1:201285/150‑1 (MQ=255)
                                                                                                                                        gagCCCGGCCCAGTGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGTAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCgtt                                        <  2:829891/120‑1 (MQ=255)
                                                                                                                                        gagCCCGGCCCAGTGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCgtt                                        >  1:829891/1‑120 (MQ=255)
                                                                                                                                        gagCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCa          >  1:723930/1‑150 (MQ=255)
                                                                                                                                            ccGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCAGCCg      >  2:663634/1‑150 (MQ=255)
                                                                                                                                                cACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCAGCCGGCCa  <  1:255508/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GGCCAGTACGGCGGCGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCAGCCGGCCA  >  NC_002947/780414‑780707

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