Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,118,728	G→C	8.3%	G802R (GGG→CGG)	PP_1884 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,118,728	0	G	C	8.3%	139.3 / 4.7	48	G802R (GGG→CGG)	PP_1884	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (20/24);  new base C (2/2);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

ATATCGATGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGCCTACCAATTACCTCTCGG  >  NC_002947/2118588‑2118871                                                                                                                                             |                                                                                                                                                atatcgatGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCCCAGCGGCAGCGGCAGcac                                                                                                                                         <  7:335706/149‑1 (MQ=255)    tcgatGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGt                                                                                                                                      >  5:309161/1‑149 (MQ=255)       atGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTc                                                                                                                                     >  4:255400/1‑147 (MQ=255)       atGAAGACGGCAAGACTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCaa                                                                                                                                   >  7:150110/1‑149 (MQ=255)         gAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGc                                                                                                                                 >  6:233262/1‑149 (MQ=255)           aGACGGCAAGCCTTACTACCAGGCCGCCGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCAATCACGCTCGAGCGTTTCATCGTGCTGCCGCTGCCCCTGGGGCAGCACAGTCAAGcc                                                                                                                                >  5:338923/1‑148 (MQ=255)               ggCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCa                                                                                                                           <  3:255400/149‑1 (MQ=255)                gCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAg                                                                                                                          <  4:22164/149‑1 (MQ=255)                 cAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAgg                                                                                                                         <  6:309161/149‑1 (MQ=255)                    gCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATc                                                                                                                            <  3:324246/143‑1 (MQ=255)                    gCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATc                                                                                                                            <  3:329545/143‑1 (MQ=255)                    gCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATc                                                                                                                            >  4:324246/1‑143 (MQ=255)                    gCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATc                                                                                                                            >  4:329545/1‑143 (MQ=255)                      cTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCtt                                                                                                                    <  2:136457/149‑1 (MQ=255)                          ctacCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAg                                                                                                                 <  2:26740/148‑1 (MQ=255)                                     aCGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGccc                                                                                                                      <  8:171175/132‑1 (MQ=255)                                     aCGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGccc                                                                                                                      >  7:171175/1‑132 (MQ=255)                                                  cGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAAc                                                                                                    >  1:242250/1‑137 (MQ=255)                                                  cGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAAc                                                                                                    <  2:242250/137‑1 (MQ=255)                                                   gCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTGGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAg                                                                                       <  5:272225/149‑1 (MQ=255)                                                              cGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCg                                                                            >  7:45918/1‑149 (MQ=255)                                                                                gAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAAcc                                                             <  8:88473/146‑1 (MQ=255)                                                                                gAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAAcc                                                             <  8:150110/146‑1 (MQ=255)                                                                                gAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCa                                                           >  7:298550/1‑148 (MQ=255)                                                                                 aGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGtt                                                                       <  4:332532/135‑1 (MQ=255)                                                                                 aGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGtt                                                                       >  3:332532/1‑135 (MQ=255)                                                                                      ctcCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCttt                                                    <  7:95753/149‑1 (MQ=255)                                                                                         ccccGACTTTCCCATCCCGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCTGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGTCGTTCACCTAAACCCATACCTTTCa                                                  <  6:338923/148‑1 (MQ=255)                                                                                         ccccGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCTGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCa                                                  <  5:233262/148‑1 (MQ=255)                                                                                         ccccGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGt                                                                                     <  6:9331/113‑1 (MQ=255)                                                                                         ccccGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGt                                                                                     >  5:9331/1‑113 (MQ=255)                                                                                               cTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTcc                                           <  5:360831/149‑1 (MQ=255)                                                                                                  aCCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCa                                          <  8:45918/147‑1 (MQ=255)                                                                                                          cGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCggg                                <  8:298550/149‑1 (MQ=255)                                                                                                                 gggTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCCCTCCAGTTGAGCGGGGCCAcat                         >  3:296005/1‑149 (MQ=255)                                                                                                                    tttCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCAcatcat                      >  6:142646/1‑149 (MQ=255)                                                                                                                       gactgTGCTGCCCCAGCGGCAGCGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGcc                   >  7:73209/5‑149 (MQ=255)                                                                                                                        aTCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGc                    >  6:306323/1‑147 (MQ=255)                                                                                                                        aTCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGc                    <  5:306323/147‑1 (MQ=255)                                                                                                                         tCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAg                                                 >  6:288402/1‑117 (MQ=255)                                                                                                                         tCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAg                                                 <  5:288402/117‑1 (MQ=255)                                                                                                                         tCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACACATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGCCTa                 >  2:298689/1‑149 (MQ=255)                                                                                                                             gctgccgctgccGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGCCTACCaa             >  7:204333/1‑149 (MQ=255)                                                                                                                                  tgctgccccagCGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCa                                          >  5:12028/2‑115 (MQ=255)                                                                                                                                  tgctgccccagCGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCa                                          <  6:12028/114‑1 (MQ=255)                                                                                                                                   gctgccgctgGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAAc                                                              <  8:164370/94‑1 (MQ=255)                                                                                                                                   gctgccgctgGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAAc                                                              >  7:164370/1‑94 (MQ=255)                                                                                                                                        cgctgGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCATCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACACATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGCCTACCAATTACCTCTCgg  <  1:298689/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                ATATCGATGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGCCTACCAATTACCTCTCGG  >  NC_002947/2118588‑2118871

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑