Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,072,519	G→C	8.9%	P330A (CCG→GCG)	mdtA ←	multidrug transporter membrane fusion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,072,519	0	G	C	8.9%	134.8 / 5.4	45	P330A (CCG→GCG)	mdtA	multidrug transporter membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base G (16/25);  new base C (2/2);  total (18/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

GACGGTCGCCGGCCTTCAAGCCGTCGAGAATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCTTTGAACTTCAGCGTGCCGGTGGTGGTGTCGATCTGGTTGTCGATGCTGCCCA  >  NC_002947/4072372‑4072665                                                                                                                                                    |                                                                                                                                                   gACGGTCGCCGGCCTTCAAGCCGTCGAGAATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGTCCGCGGCCggcg                                                                                                                                                  >  1:156120/1‑148 (MQ=255)    ggTCGCCGGCCTTCAAGCCGTCGAGAATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCTGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGccg                                                                                                                                               <  2:548385/150‑2 (MQ=255)             ccTTCAAGCCGTCGAGAATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATCGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCt                                                                                                                                      <  2:372415/150‑1 (MQ=255)              cTTCAAGCCGTCGAGAATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGa                                                                                                                                            >  2:85551/1‑143 (MQ=255)              cTTCAAGCCGTCGAGAATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGa                                                                                                                                            <  1:85551/143‑1 (MQ=255)                 cAAGCCGTCGAGAATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGttt                                                                                                                                  <  2:674173/150‑1 (MQ=255)                     ccGTCGAGAATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAgg                                                                                                                              <  2:306438/150‑1 (MQ=255)                          gagaATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGaccac                                                                                                                                        <  1:101901/135‑1 (MQ=255)                          gagaATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGaccac                                                                                                                                        >  2:101901/1‑135 (MQ=255)                             aaTGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGc                                                                                                                      <  2:162367/150‑1 (MQ=255)                               tGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCa                                                                                                                    >  1:55803/1‑150 (MQ=255)                                gACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGt                                                                                                                                    <  1:152717/133‑1 (MQ=255)                                gACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGt                                                                                                                                    >  2:152717/1‑133 (MQ=255)                                          tcGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCCCGGCCGCGGCCAtcca                                                                                                                                           <  1:198486/116‑4 (MQ=255)                                          tcGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCCCGGCCGCGGCCAtcca                                                                                                                                           >  2:198486/1‑113 (MQ=255)                                            gCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGtt                                                                                                       >  1:786008/1‑150 (MQ=255)                                              cATCGCTGGCACGCACCTTCCGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGCCCATGACCACCTGTTTCAGGGTGTCGGCCAGCAAGCGCACCTTGa                                                                                                     <  2:336172/150‑1 (MQ=255)                                                   cTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCCGCAAGCGCACGTTGAGGAAc                                                                                                >  2:95646/1‑150 (MQ=255)                                                   cTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAAc                                                                                                <  2:786008/150‑1 (MQ=255)                                                    tGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACt                                                                                               >  1:485725/1‑150 (MQ=255)                                                       cacccaccTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGt                                                                                            >  1:61649/1‑150 (MQ=255)                                                         cccaccTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTg                                                                                          >  2:604322/1‑150 (MQ=255)                                                           caccTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTggg                                                                                        <  2:549427/150‑1 (MQ=255)                                                           caccTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTggg                                                                                        <  2:485725/150‑1 (MQ=255)                                                           caccTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGNTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTggg                                                                                        <  1:604322/150‑1 (MQ=255)                                                                    cTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGc                                                                               <  2:55803/150‑1 (MQ=255)                                                                              ttGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGAc                                                                                                    >  2:329192/1‑119 (MQ=255)                                                                              ttGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGAc                                                                                                    <  1:329192/119‑1 (MQ=255)                                                                                gACGGTGCTCTGCTCATTCACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGGGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGtt                                                                   <  2:343830/150‑1 (MQ=255)                                                                                 aCGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTtc                                                                  >  2:163972/1‑150 (MQ=255)                                                                                   ggTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTtctc                                                                <  1:785097/150‑1 (MQ=255)                                                                                        tctgctcATTGACAACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAAc                                                           >  2:631756/1‑150 (MQ=255)                                                                                              cATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCtt                                                     >  1:117792/1‑150 (MQ=255)                                                                                                 tGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCTTTGa                                                  <  2:117792/150‑1 (MQ=255)                                                                                                         tAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCTTTGAACTTCAGc                                          <  1:631756/150‑1 (MQ=255)                                                                                                                  gTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCTTTGAACTTCAGCGTGCCggtg                                 <  1:163972/150‑1 (MQ=255)                                                                                                                     ccATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCTTTGAACTTCAGCGTGCCggtggtg                              <  2:61649/150‑1 (MQ=255)                                                                                                                      cATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTggg                                                                                        <  1:334154/91‑1 (MQ=255)                                                                                                                      cATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTggg                                                                                        >  2:334154/1‑91 (MQ=255)                                                                                                                      cATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGaaa                                                            <  2:21729/119‑1 (MQ=255)                                                                                                                      cATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGaaa                                                            >  1:21729/1‑119 (MQ=255)                                                                                                                             gCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCTTTGAACTTCAGCGTGCCGGTGGTGGTGTCGAt                      <  1:396303/150‑1 (MQ=255)                                                                                                                                 aaCTGGATGGCCGCGGCCGGGGCCATGACCGCCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCTTTGAACTTCAGCGTGCCGGTGGTGGTGTCGATCTgg                  <  1:592672/150‑1 (MQ=255)                                                                                                                                       aTGGCCCCGGCCGCGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCttt                                                    <  2:347256/110‑1 (MQ=255)                                                                                                                                       aTGGCCCCGGCCGCGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCttt                                                    >  1:347256/1‑110 (MQ=255)                                                                                                                                                 ccGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCTTTGAACTTCAGCGTGCCGGTGGTGGTGTCGATCTGGTTGTCGATGCTGCCCa  >  1:538249/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GACGGTCGCCGGCCTTCAAGCCGTCGAGAATGACGCTGTTCTCGCCATCGCTGGCACCCACCTTCAGCTTGCGGACATTGACGGTGCTCTGCTCATTGACCACATAGGCGAAGGTGCCATCGTTGCCAAACTGGATGGCCGCGGCCGGGGCCATGACCACCTGTTTGAGGGTGTCGGCCAGCAAGCGCACGTTGACGAACTGGTTGGGGAACAGGGCCAGGTCCTTGTTCTCGAAACGGCCTTTGAACTTCAGCGTGCCGGTGGTGGTGTCGATCTGGTTGTCGATGCTGCCCA  >  NC_002947/4072372‑4072665

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