Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,897,770 A→G 30.6% intergenic (‑170/+64) 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,7700AG30.6% 45.0 / 25.1 36intergenic (‑170/+64)PP_5172/PP_5173hypothetical protein/RND family transporter
Reads supporting (aligned to +/- strand):  ref base A (12/13);  new base G (4/7);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

GCCGGGCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGC  >  NC_002947/5897625‑5897914
                                                                                                                                                 |                                                                                                                                                
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                                                                                                                       <  7:5900/148‑1 (MQ=25)
                         tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCCCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTcgc                                                                                                                       <  5:247847/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                                                                                                                     >  5:43601/1‑127 (MQ=21)
                                                tagtGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                     <  6:43601/127‑1 (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            <  7:119914/143‑1 (MQ=37)
                                                                                                                                         ccGCGAATGCGATGGTGCATCGACCACCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCa            >  8:119914/1‑143 (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)
                                                                                                                                                 |                                                                                                                                                
GCCGGGCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGC  >  NC_002947/5897625‑5897914

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: