Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,897,773 G→A 32.3% intergenic (‑173/+61) 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,7730GA32.3% 45.2 / 22.6 34intergenic (‑173/+61)PP_5172/PP_5173hypothetical protein/RND family transporter
Reads supporting (aligned to +/- strand):  ref base G (11/12);  new base A (4/7);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GCCGGGCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGAT  >  NC_002947/5897625‑5897920
                                                                                                                                                    |                                                                                                                                                   
gcCGGGCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACgg                                                                                                                                                     >  8:106091/1‑149 (MQ=255)
   cccatggcctacgccatgattggcgggatcatcgtggcgaccttgctgacgctgctgttcctgccggcgctgtacgtggcctggtacaagatccgcgagccgcagaAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACggtgg                                                                                                                                                  <  8:214293‑M1/43‑1 (MQ=255)
            tttGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCAc                                                                                                                                         <  7:38271/149‑1 (MQ=255)
                aaGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCAt                                                                                                                                      <  8:341227/148‑1 (MQ=255)
                aaGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATc                                                                                                                                     <  7:77442/149‑1 (MQ=255)
                 atgggTTTTGCGCCTGGCTTTAGGGCAATAGTAGTGATGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATc                                                                                                                                     <  8:183576/146‑1 (MQ=255)
                        ggcgggatcatcgtggcgaccttgctgacgctgctgttcctgccggcgctgtacgtggcctggtacaagatccgcgagccgcagaAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTcg                                                                                                                              >  3:396461‑M1/86‑148 (MQ=255)
                         tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCCCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTcgc                                                                                                                             <  5:247847/148‑1 (MQ=25)
                         tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCCCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTcgc                                                                                                                             <  7:5900/148‑1 (MQ=25)
                         gcgggatcatcgtggcgaccttgctgacgctgctgttcctgccggcgctgtacgtggcctggtacaagatccgcgagccgcagaAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTcgcg                                                                                                                            >  2:80655‑M1/85‑149 (MQ=255)
                         gcgggatcatcgtggcgaccttgctgacgctgctgttcctgccggcgctgtacgtggcctggtacaagatccgcgagccgcagaAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTcgc                                                                                                                             >  7:214293‑M1/85‑148 (MQ=255)
                           cgcCTGGCTTGCTGGGAATAGTAGTGAGGGGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATTTAACCATTTAAGAGATCACGCCCCTGGTAGGAGCGGGTTTACCCGCGAATGCTATGGTGCATCGACCACCGTATTCGCgg                                                                                                                           <  3:48100/148‑1 (MQ=18)
                           cgcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                           <  1:229525/148‑1 (MQ=25)
                           cgcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                           <  1:139263/148‑1 (MQ=255)
                           cgcCTGGCTCGCCGGCAATAGTACTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                           <  8:411993/148‑1 (MQ=255)
                           cgacTGGCTTGAGGGCAATACTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                           <  8:203284/145‑1 (MQ=255)
                                 gCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAAc                                                                                                                     <  1:158946/148‑1 (MQ=255)
                                                tagtGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                           <  6:43601/127‑1 (MQ=21)
                                                tagtGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                           >  5:43601/1‑127 (MQ=21)
                                                 agtGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCg                                                                                                                  <  8:283325/135‑1 (MQ=255)
                                                 agtGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCg                                                                                                                  >  7:283325/1‑135 (MQ=255)
                                                       gCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCCCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                           <  1:428630/120‑1 (MQ=18)
                                                       gCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCCCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                           >  2:428630/1‑120 (MQ=18)
                                                         gggATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATggg                                                                                            >  8:7716/1‑149 (MQ=255)
                                                                                   aGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTcg                                                                   <  7:106091/148‑1 (MQ=255)
                                                                                                          gagaACACGCCTCTGGTAGGAGCGGGTTTAACCGCGAATACGGTGGTAGATGCACCATAGCATTAGCGGGTAAAACCGCTCCTACCAGGGGCCATATGGGCCTTAAGGCGTTTTTATGAGGCTCGGGGGTATTGTAACAGGCCACGTa                                            >  8:266114/1‑148 (MQ=255)
                                                                                                               CACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAgc                                        >  2:310640‑M1/1‑147 (MQ=255)
                                                                                                               CACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAgc                                        <  1:310640‑M1/147‑1 (MQ=255)
                                                                                                                    CTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTTAACCCGCGCCTACCAGGGGCCATATGGGGCGTCAGGCGTGTTTCTGCGGCTCGCGGATCGTGTACCCGGCCACGTACCGCGCCgg                                   <  4:260590‑M1/147‑1 (MQ=255)
                                                                                                                                   gTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTGTGCGGCTCGAGGATCTTGTCCCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGt                    >  7:201584/1‑147 (MQ=255)
                                                                                                                                         ccGCGAATGCGATGGTGCATCGACCACCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCa                  >  8:119914/1‑143 (MQ=37)
                                                                                                                                         ccGCGAATGCGATGGTGCATCGACCACCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCa                  <  7:119914/143‑1 (MQ=37)
                                                                                                                                          cgcgAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt           >  4:110787/1‑149 (MQ=255)
                                                                                                                                           gcgAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTc          >  3:832/1‑149 (MQ=37)
                                                                                                                                           gcgAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTc          >  4:160098/1‑149 (MQ=255)
                                                                                                                                             gAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGc        >  3:6056/1‑149 (MQ=255)
                                                                                                                                                    gtggtCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACgat  >  3:86888/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCCGGGCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGAT  >  NC_002947/5897625‑5897920

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: