Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 780,561 T→A 8.8% intergenic (‑88/+34) 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,5610TA 8.8% 134.2 / 3.8 46intergenic (‑88/+34)PP_0670/glyA‑IIbile acid/Na+ symporter family transporter/serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base T (22/19);  new base A (2/2);  total (25/21)
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.83e-01

CGGCCAGTACGGCGGCGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCAGCCG  >  NC_002947/780413‑780703
                                                                                                                                                    |                                                                                                                                              
cGGCCAGTACGGCGGCGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTg                                                                                                                                               <  1:920701/150‑1 (MQ=255)
             ggcgAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGCGGCGGTGCGCCTGTATTAACCGCAAGAGCCCGGCCTAGTGCCGGGCTCttt                                                                                                                                  <  1:629327/150‑1 (MQ=255)
                      gATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAgg                                                                                                                         <  2:212786/150‑1 (MQ=255)
                           ggATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAgc                                                                                                                    <  1:315645/150‑1 (MQ=255)
                           ggATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAgc                                                                                                                    <  1:586878/150‑1 (MQ=255)
                                       tttCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTg                                                                                                        <  1:693247/150‑1 (MQ=255)
                                             gTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCATGCTTTAAACTAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTaa                                                                                                  >  1:74813/1‑150 (MQ=255)
                                                  tctctGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCgg                                                                                             >  2:961593/1‑150 (MQ=255)
                                                    tctGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCa                                                                                           >  1:206875/1‑150 (MQ=255)
                                                    tctGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCa                                                                                           <  1:736494/150‑1 (MQ=255)
                                                               tGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAg                                                                                >  2:703604/1‑150 (MQ=255)
                                                               ggCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAg                                                                                <  1:474323/149‑1 (MQ=255)
                                                                   gCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGcttc                                                                                      <  2:1035109/140‑1 (MQ=255)
                                                                   gCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGcttc                                                                                      >  1:1035109/1‑140 (MQ=255)
                                                                                 gAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCa                                                              >  1:852403/1‑150 (MQ=255)
                                                                                      aTTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCacggac                                                         <  1:805346/150‑1 (MQ=255)
                                                                                                   tCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTg                                                                                                        <  1:1043326/90‑1 (MQ=255)
                                                                                                   tCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTg                                                                                                        >  2:1043326/1‑90 (MQ=255)
                                                                                                     tttGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAg                                                                                                                          <  1:697232/70‑1 (MQ=255)
                                                                                                     tttGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAg                                                                                                                          >  2:697232/1‑70 (MQ=255)
                                                                                                     tttGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTGGCAGAGGGCCTTGACCTTCGCACGTACGGCGGCTATCACCGc                                          >  2:665346/1‑150 (MQ=255)
                                                                                                            gAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCgttg                                   >  1:1138078/1‑150 (MQ=255)
                                                                                                                ggCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCa                               >  1:131199/1‑150 (MQ=255)
                                                                                                                 gCGGTGCGGCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGACAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAg                              >  1:387485/1‑150 (MQ=255)
                                                                                                                     tGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCg                          <  2:206875/150‑1 (MQ=255)
                                                                                                                         ccTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCa                      >  1:800546/1‑150 (MQ=255)
                                                                                                                          cTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAg                                                                                          <  2:642223/81‑1 (MQ=255)
                                                                                                                          cTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAg                                                                                          >  1:642223/1‑81 (MQ=255)
                                                                                                                              aTTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATg                 >  2:755487/1‑150 (MQ=255)
                                                                                                                               ttAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGAt                                                                              <  2:218185/88‑1 (MQ=255)
                                                                                                                               ttAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGAt                                                                              >  1:218185/1‑88 (MQ=255)
                                                                                                                               ttAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGt                >  2:428726/1‑150 (MQ=255)
                                                                                                                                 aaaCGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTc                                                   >  2:535551/1‑113 (MQ=255)
                                                                                                                                 aaaCGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTc                                                   <  1:535551/113‑1 (MQ=255)
                                                                                                                                 aaaCGAAAGAGCCCGGCACTGGGCCGGGCTCTTTAACTGGAGTCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCg              <  1:665346/150‑1 (MQ=255)
                                                                                                                                  aaCGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGc             >  2:707385/1‑150 (MQ=255)
                                                                                                                                  aaCGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGc             >  2:826486/1‑150 (MQ=255)
                                                                                                                                   aCGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTc               >  2:391743/1‑147 (MQ=255)
                                                                                                                                   aCGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTc               <  1:391743/147‑1 (MQ=255)
                                                                                                                                      aaaGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCATAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGAt         >  1:1063124/1‑150 (MQ=255)
                                                                                                                                        agagCCCGGCCCAGTGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCgtt                                    <  2:490180/121‑1 (MQ=255)
                                                                                                                                        agagCCCGGCCCAGTGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCgtt                                    >  1:490180/1‑121 (MQ=255)
                                                                                                                                         gagCCCGGCCCAGTGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCa      >  1:469156/1‑150 (MQ=255)
                                                                                                                                           gCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCAGc    <  2:1138078/150‑1 (MQ=255)
                                                                                                                                            aacGGCACGGGGCCGGTCTCTTTCACTTTACGCAAGCGGCACACAGTTGCCGTAAACCGGCAGCTTCTTGCGGATGTCCTTGACCTTCCCACGGACGGCGTCGATGACCGGCTCGTTTTTCAGGTCGGGCAGGTTGTCGCCGATCCCGcc   >  2:75576/3‑150 (MQ=255)
                                                                                                                                             ccGGCACTGGGCCGGGCTCTTTACCTTTAGGCAACAGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCAGCCg  <  2:800546/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CGGCCAGTACGGCGGCGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCAGCCG  >  NC_002947/780413‑780703

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