Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,897,779 G→C 9.8% intergenic (‑179/+55) PP_5172 ← / ← PP_5173 hypothetical protein/RND family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,897,7790GC 9.8% 125.8 / 3.9 50intergenic (‑179/+55)PP_5172/PP_5173hypothetical protein/RND family transporter
Reads supporting (aligned to +/- strand):  ref base G (24/21);  new base C (2/3);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

AAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCAC  >  NC_002947/5897641‑5897917
                                                                                                                                          |                                                                                                                                           
aaGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATc                                                                                                                                   <  6:221831/149‑1 (MQ=255)
aaGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCA‑TTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCg                                                                                                                                  <  4:377476/149‑1 (MQ=255)
 aGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTTCCCGCGAATACGGTGGTGGGTGGCCCATCg                                                                                                                                  >  8:152355/1‑149 (MQ=255)
 aGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGGATTAGGTGGTCGATGCACCATCg                                                                                                                                  >  3:20888/1‑149 (MQ=255)
              cTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTa                                                                                                                      >  6:90641/1‑148 (MQ=255)
              cTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTAGCGGGTa                                                                                                                      >  3:151700/1‑148 (MQ=255)
                 gCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTTACCGCGAATACGGTGGTCGATGGACCATCGCATTCGCGGGGaaaac                                                                                                                  >  6:50440/1‑147 (MQ=255)
                   ttGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCg                                                                                                                >  2:270054/1‑149 (MQ=255)
                      agtgCAATAGTAGTTATGCGGCATTATTATTAAAAATACTGTTTAAGAATTTTTTTAAAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCt                                                                                                              <  6:216739/145‑1 (MQ=38)
                      aGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTc                                                                                                             <  3:291268/149‑1 (MQ=255)
                      aGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCATCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTc                                                                                                             >  5:268167/1‑149 (MQ=255)
                          cAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTAc                                                                                                         >  6:157749/1‑149 (MQ=255)
                           aaTAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGGAAACCCGCTCCTAc                                                                                                         >  8:101793/1‑148 (MQ=255)
                              agaagtCAGGCGGGATAATTCTTAAAAATACTGTTTAAGAATTTTTAAAAAACCATTAAAGACACCACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAgg                                                                                                     <  7:152355/146‑1 (MQ=37)
                               gtagtGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAgg‑g                                                                                                   <  7:101793/149‑1 (MQ=255)
                                        cGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATgg                                                                                          <  5:157749/149‑1 (MQ=255)
                                         gggATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGAGGGTTTACCCGCGAAAACCGTGGTCGAAGCACAATCGCATTAGCCGGTAAACACCATCCTACCAGG‑GGCCATATggg                                                                                         >  6:136174/1‑149 (MQ=255)
                                                 tCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAgg                                                                                 >  7:332924/1‑149 (MQ=255)
                                                 tCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAg                                                                                  <  5:90641/148‑1 (MQ=255)
                                                   ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCg                                                                               >  4:284234/1‑149 (MQ=255)
                                                   ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGGCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCGACCAGG‑GGCCATAGGGGCCTGCAGGCg                                                                               >  8:249084/1‑149 (MQ=255)
                                                     aaaaaTACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCgtg                                                                             >  2:56772/1‑149 (MQ=255)
                                                             tGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCg                                                                     >  3:329940/1‑149 (MQ=255)
                                                                 tAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCTAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTAc                                                                                                         <  4:255459/110‑1 (MQ=255)
                                                                 tAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTAGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTAc                                                                                                         >  3:255459/1‑110 (MQ=255)
                                                                          ttATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGCGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATc                                                          <  1:270054/148‑1 (MQ=255)
                                                                           tatataACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCtt                                                        <  1:56772/148‑1 (MQ=255)
                                                                               taACCATTTAAGTGAACACGCCTCTGGTAGCAGCGGGTTTACCCGCTAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTAcc                                                   <  5:136174/149‑1 (MQ=255)
                                                                                 aCCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAg                                                 >  1:322453/1‑149 (MQ=255)
                                                                                     tttAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCGGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAg                                                 <  6:147379/145‑1 (MQ=255)
                                                                                     tttAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAg                                                 >  5:147379/1‑145 (MQ=255)
                                                                                       tAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCAc                                            <  3:289317/148‑1 (MQ=255)
                                                                                         agagAACACGCCGCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTa                                         <  4:329940/149‑1 (MQ=255)
                                                                                              acccGCCTCTGGTAGGTGCGGGTTTACCCGCGATTAGGGTGGTCGATGCACCTTCGCATTCGCGGTTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAgcg                                    <  4:151700/146‑1 (MQ=255)
                                                                                                   ccTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCgg                                >  2:238943/1‑148 (MQ=255)
                                                                                                     tctGGTAGGAGCGGGTTTTACCGCGAATACGGTGGTCGATGGACCATCGCATTAGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGGCTTCAGGCGTGTTTCTGCGGGTCGCGGATCTTGTACCAGGCCACGTACAGAGCCGGGAg                             >  8:207145/1‑149 (MQ=37)
                                                                                                       tGGTAGGAGCGGTTTACCCCGCTAATACGGTGGTCGAGCCACCATCGAATTCGCGGGAAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAg                             <  5:39645/147‑1 (MQ=255)
                                                                                                             gcgCGGGTTTACCCGCTAATACGGGGGTCCATCCCCCATCGCTTTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAcagca                     <  2:322453/147‑1 (MQ=38)
                                                                                                             gAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAcag                       >  5:14386/1‑147 (MQ=255)
                                                                                                              aGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAcagcag                    >  8:101849/1‑149 (MQ=21)
                                                                                                               gCGGTTTAACCCGCTAATAGGGTGGCGCATGACCCACGGCTTTGCCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACagcagc                   <  3:284234/149‑1 (MQ=255)
                                                                                                               gCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCCTTCGCGGGTTAACCCGCTGCTACCAGG‑GGCCATATGGGGCGTAAGGCGTGTTTATGCGGCGCGCGGATCGTGTACCAGGCCACGTTCAGCGCCGGCAGGCACagcagc                   >  6:56212/1‑149 (MQ=255)
                                                                                                               gCGGGTTTACCCGCGAATACGGTGGGCGATGAACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACagcagc                   <  4:20888/149‑1 (MQ=255)
                                                                                                                        cccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAg          <  8:332924/149‑1 (MQ=255)
                                                                                                                         ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAgg         >  8:114280/1‑149 (MQ=255)
                                                                                                                          cgcgAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAcag                       >  3:220214/1‑134 (MQ=25)
                                                                                                                          cgcgAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAcag                       <  4:220214/134‑1 (MQ=25)
                                                                                                                              aagaCGGTGGTCGACGCACCATCTCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGcc    <  6:268167/146‑1 (MQ=255)
                                                                                                                              aaTACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGcc    <  7:207145/149‑1 (MQ=255)
                                                                                                                                 aCGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCAc  >  2:365383/1‑148 (MQ=255)
                                                                                                                                          |                                                                                                                                           
AAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCAC  >  NC_002947/5897641‑5897917

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: ‑