Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,897,778	C→G	10.4%	intergenic (‑178/+56)	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_002947	5,897,778	0	C	G	10.4%	111.6 / 3.3	46	intergenic (‑178/+56)	PP_5172/PP_5173	hypothetical protein/RND family transporter
Reads supporting (aligned to +/- strand):  ref base C (19/20);  new base G (3/2);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-01

GAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGA‑GAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGAT  >  NC_002947/5897640‑5897920                                                                                                                                            |                                                                                                                                               gAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGA‑GAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCAt                                                                                                                                      >  1:15795/1‑149 (MQ=255)    ggggTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGA‑GAACACGCCTCTGGTAGGAGCGGGTTTTCCCGCGAATACGGTGGTCGATGCACCAtag                                                                                                                                    >  7:116429/1‑146 (MQ=255)             cgcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGA‑GAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                           <  8:42139/148‑1 (MQ=25)             cgcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGA‑GAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCggg                                                                                                                          <  6:121460/149‑1 (MQ=255)             cgacTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACAGTTTAAGAATTTTAATATAACCATTTAAGA‑GAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                           <  6:272685/145‑1 (MQ=25)                cTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGA‑GAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTa                                                                                                                        >  5:214510/1‑148 (MQ=255)                     ttGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGA‑GAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCa                                                                                                                                       <  6:233373/128‑1 (MQ=255)                     ttGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGA‑GAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCa                                                                                                                                       >  5:233373/1‑128 (MQ=255)                         gggCAATAGTAGTGAGGCGGTATTATTATTAAAAATACTGTTTAATAATTTTTATATAACCATTTACGA‑GACCACGCCTCTGGTCGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTc                                                                                                               <  6:214510/148‑1 (MQ=38)                          ggCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGATAAAAACCCCTATGGTATGAGCCGGTTTAACCGCGAATACGGTTGTGGTTGCACAAAACCATTACCGGGTAAAccccaaa                                                                                                               >  1:107278/1‑144 (MQ=255)                                 gtagtGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGA‑GAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAgg                                                                                                       <  2:71183/148‑1 (MQ=255)                                       aGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGA‑GAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTAGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCata                                                                                                >  7:143303/1‑149 (MQ=255)                                                     ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGA‑GAACTCGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCg                                                                                  >  1:254915/1‑149 (MQ=255)                                                     ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGA‑GAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCg                                                                                  >  3:256093/1‑149 (MQ=255)                                                         aaaTACTGTTTAAGAATTTTTATATAACCATTTAAGA‑GAACACGCCTCTGGTAGGAGCGGGTTTAACCGCGAATAAGGTGGTAGATTGACCATAGCATTAGCGGGTAAAACCCCACCTACCAGGGGGCCAATGGGCCTTAAGGCGTGtt                                                                              >  7:117716/1‑149 (MQ=255)                                                                gTTTAAGAATTTTTATATAACCATTTAAGA‑GAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCgg                                                                       >  1:15606/1‑149 (MQ=255)                                                                 tttAAGAATTTTTATATAACCATTTAAGA‑GAACACGCCTCTGGTAGGAGCGGGTTTTACCGCGAATACGGGGGTGGCTGGCCCCTCGCCTTGGCGGGGTAACCCGCGCCGGCCCGGGGGGAGATGGGGCGTGAGGGGGGTTTGTggggg                                                                      >  8:146775/1‑148 (MQ=255)                                                                  ttAAGAATTTTTATATAACCATTTAAGA‑GAACACGCCTCTGGTAGGAGCGGGTTTAACCGCGAATACGGTGGTCGATTCACCATCGCATTCGCgg                                                                                                                           >  5:90293/1‑95 (MQ=255)                                                                   tAAGAATTTTTATATAACCATTTAAGA‑GAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCt                                                                     <  2:107278/148‑1 (MQ=255)                                                                            ttATATAACCATTTAAGA‑GAACACGCCTCTGGTAGGAGCGGGTTTACCCGCTAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCt                                                            <  8:143303/148‑1 (MQ=255)                                                                                  aaCCATTTAAGA‑GAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTAc                                                       >  4:122414/1‑147 (MQ=255)                                                                                          aaGA‑GAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACg                                              >  2:217810/1‑148 (MQ=255)                                                                                          aaGA‑GAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCACCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACg                                              >  1:270683/1‑148 (MQ=255)                                                                                           aga‑gAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGt                                             <  1:217810/148‑1 (MQ=255)                                                                                            ga‑gaACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTAc                                           >  7:108596/1‑149 (MQ=255)                                                                                                   ccGCCTCTGGTCGGAGGGGGTTTACCCGCTAATACGGTGGTCGTTGCACCATCGCTTTCGCGGTTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGcc                                     <  8:143472/148‑1 (MQ=255)                                                                                                     gCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTcgc                                                                  <  6:96423/118‑1 (MQ=255)                                                                                                     gCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTcgc                                                                  >  5:96423/1‑118 (MQ=255)                                                                                                     gCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCGTGAGGCGTGGTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCgg                                   >  7:210470/1‑149 (MQ=255)                                                                                                        cctGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAg                                >  5:165133/2‑149 (MQ=21)                                                                                                         ctGGCAGGAGCGGGTTTACCCGCGAATACGGGGGCCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAgg                               <  8:210470/149‑1 (MQ=255)                                                                                                             tAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAca                           <  4:256093/149‑1 (MQ=255)                                                                                                                gAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAcagca                        <  8:108596/149‑1 (MQ=255)                                                                                                                 aGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAcagcag                       >  5:103363/1‑149 (MQ=255)                                                                                                                  gCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACagcagc                      <  4:282366/149‑1 (MQ=255)                                                                                                                         tccccGCTAATAGGGTGGGCGCTGCACCATCGCATTCGCGGTTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGc                <  2:176374/146‑1 (MQ=255)                                                                                                                           cccGCTAATACGGTGGTCGATGCACATTCGCATTCGCGGGAAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAg             <  8:117716/149‑1 (MQ=255)                                                                                                                            ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACagcagc                      <  6:274626/139‑1 (MQ=255)                                                                                                                            ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACagcagc                      >  5:274626/1‑139 (MQ=255)                                                                                                                             cgcgAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt           >  3:130077/1‑149 (MQ=255)                                                                                                                             cgcgAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCACATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCgg                                   >  7:200974/1‑125 (MQ=255)                                                                                                                             cgcgAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCACATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCgg                                   <  8:200974/125‑1 (MQ=255)                                                                                                                                gAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGc        <  2:15795/149‑1 (MQ=255)                                                                                                                                   tACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAg                                                    <  3:34321/102‑1 (MQ=255)                                                                                                                                   tACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAg                                                    >  4:34321/1‑102 (MQ=255)                                                                                                                                      ggtggtCGATGCACCATCGCATTCGCGGGTAAACCCGCACCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACgat  <  2:270683/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                               GAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGA‑GAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGAT  >  NC_002947/5897640‑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: ‑