Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,405,237	A→C	10.2%	Q78H (CAA→CAC)	PP_2108 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,405,237	0	A	C	10.2%	143.7 / 7.3	49	Q78H (CAA→CAC)	PP_2108	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (20/24);  new base C (2/3);  total (22/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 = 9.95e-01

GCCAGGGCCAACCTGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACG  >  NC_002947/2405094‑2405375                                                                                                                                                 |                                                                                                                                           gCCAGGGCCAACCTGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGAGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGc                                                                                                                                        <  1:76387/149‑1 (MQ=255) gCCAGGGCCAACCTGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCt                                                                                                                                       >  7:39510/1‑149 (MQ=255)        ccAACCTGCAATTACTGTGTGGTGCGCAGCAACAGGTACGTGCGTTGTTTGACTCGCCCATCAACATCAGCAGCGGGTATCTCAGCGAAAGG‑GTCAATCGCTTGATCGGTGCAGCCTCCGATACCCAGCACGTGCAAGGTCTGgcggcg                                                                                                                                >  8:169020/1‑149 (MQ=255)         cAACCTGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCCCTTGCACGTGCTggctatcg                                                                                                                               <  5:70301/149‑6 (MQ=255)         cAACCTGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGgcggcg                                                                                                                                >  8:246069/1‑148 (MQ=255)                              tgCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCa                                                                                                          <  7:246069/149‑1 (MQ=255)                               gcgcTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAg                                                                                                         >  6:151646/1‑149 (MQ=255)                                     gAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGt                                                                                                            <  6:27293/140‑1 (MQ=255)                                     gAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGt                                                                                                            >  5:27293/1‑140 (MQ=255)                                          ggTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCgaga                                                                                              <  5:151646/149‑1 (MQ=255)                                                tGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTcc                                                                                        >  6:184777/1‑149 (MQ=255)                                                 gCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCg                                                                                       >  5:75164/1‑149 (MQ=255)                                                   gttgttTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAgg                                                                                                             <  7:2264/125‑1 (MQ=255)                                                   gttgttTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAgg                                                                                                             >  8:2264/1‑125 (MQ=255)                                                    ttgttTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGc                                                                                     >  5:236423/1‑148 (MQ=255)                                                                  catcatcatCAGCAGCGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTg                                                                      <  2:209798/149‑1 (MQ=255)                                                                      atcctcaGCTGCGGGTATCGCAGCTAAAGG‑GTCAAGAGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCACGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGt                                                                  <  8:268182/149‑1 (MQ=255)                                                                         atcaGCAGCGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAg                                                                        <  6:52291/140‑1 (MQ=255)                                                                         atcaGCAGCGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAg                                                                        >  5:52291/1‑140 (MQ=255)                                                                          tcaGCAGCGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGccc                                                              >  4:265581/1‑149 (MQ=255)                                                                            agcagcGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCtt                                                            <  2:71537/149‑1 (MQ=255)                                                                            agcagcGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCt                                                             <  7:265007/148‑1 (MQ=255)                                                                              cagcGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTc                                                           >  7:290612/1‑148 (MQ=255)                                                                              cagcGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCg                                                          <  3:169882/149‑1 (MQ=255)                                                                                gcGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGAt                                                        <  6:75164/149‑1 (MQ=255)                                                                                 cGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCga                                                                                                <  3:61157/108‑1 (MQ=255)                                                                                 cGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCga                                                                                                >  4:61157/1‑108 (MQ=255)                                                                                    gTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCg                                                                    <  6:246630/133‑1 (MQ=255)                                                                                    gTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCg                                                                    >  5:246630/1‑133 (MQ=255)                                                                                    gTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGc                                                    <  6:236423/149‑1 (MQ=255)                                                                                                 aGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAAtt                                       <  8:39510/149‑1 (MQ=255)                                                                                                  gg‑gTCAATCGGTTGTTCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACCTCACGGTGATCGAGGTCAGCCCTAGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTc                                      <  7:169020/149‑1 (MQ=255)                                                                                                   g‑gTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCg                                     <  3:265581/149‑1 (MQ=255)                                                                                                     gTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCg                                     >  7:106231/1‑148 (MQ=255)                                                                                                      tCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGAc                                   <  5:184777/149‑1 (MQ=255)                                                                                                          tCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCg                                                                                                                >  8:104236/1‑68 (MQ=255)                                                                                                          tCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCg                                                                                                                <  7:104236/68‑1 (MQ=255)                                                                                                             gTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTc                                      >  6:156889/1‑139 (MQ=255)                                                                                                             gTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTc                                      <  5:156889/139‑1 (MQ=255)                                                                                                                  tCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGc                                                                                                        >  4:287703/1‑68 (MQ=255)                                                                                                                  tCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGc                                                                                                        <  3:287703/68‑1 (MQ=255)                                                                                                                     gTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTg                    <  8:290612/149‑1 (MQ=255)                                                                                                                                 gccgCCAGCCCTTGCACGTGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGAcgcgc        >  7:6636/4‑149 (MQ=255)                                                                                                                                 aTAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGAcgcgc        >  2:90052/1‑149 (MQ=255)                                                                                                                                    gCCAGCCCTTGCACGTGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGc     <  8:6636/149‑1 (MQ=255)                                                                                                                                    gCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCgg      <  1:90052/148‑1 (MQ=255)                                                                                                                                     ccAGCCCTTGCACGTGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCa    >  8:58267/1‑149 (MQ=255)                                                                                                                                      cAGCCCTTGCACGTGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCAc   <  7:58267/149‑1 (MQ=255)                                                                                                                                       aGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCAcg  >  1:176192/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                           GCCAGGGCCAACCTGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGG‑GTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACG  >  NC_002947/2405094‑2405375

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