Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,151,982	A→T	9.3%	intergenic (+128/+23)	PP_3655 → / ← PP_3656	cytosine/purine/uracil/thiamine/allantoin permease family protein/aromatic compound‑specific porin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,151,982	0	A	T	9.3%	131.6 / 4.7	43	intergenic (+128/+23)	PP_3655/PP_3656	cytosine/purine/uracil/thiamine/allantoin permease family protein/aromatic compound‑specific porin
Reads supporting (aligned to +/- strand):  ref base A (18/21);  new base T (2/2);  total (20/23)
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.92e-01

GGCAGCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTCTGAA  >  NC_002947/4151839‑4152123                                                                                                                                                 |                                                                                                                                              ggCAGCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAggc                                                                                                                                         >  1:229266/1‑150 (MQ=255)    aGCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAggccgg                                                                                                                                      >  1:169400/1‑150 (MQ=255)               aaTCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCt                                                                                                                           >  2:1254253/1‑150 (MQ=255)                aTCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGGCGGCCCGTTTCGCTa                                                                                                                          >  1:967105/1‑150 (MQ=255)                aTCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTa                                                                                                                          <  1:1267959/150‑1 (MQ=255)                      aCTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGa                                                                                                                    <  2:229266/150‑1 (MQ=255)                              ttCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCAGTTTCGCTATTTAGAAGTTGAAg                                                                                                            <  1:432047/150‑1 (MQ=255)                                       aaaGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAggccggcc                                                                                                                                    <  1:1121720/117‑1 (MQ=255)                                       aaaGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAggccggcc                                                                                                                                    >  2:1121720/1‑117 (MQ=255)                                        aaGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAACGGGCCGGCCTTCTGAA‑GCCGGcccatt                                                                                                                                >  2:1276221/1‑117 (MQ=38)                                        aaGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAACGGGCCGGCCTTCTGAA‑GCCGGCCCAtt                                                                                                                                <  1:1276221/120‑1 (MQ=38)                                            tGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGaa                                                                                                             >  2:329509/1‑135 (MQ=255)                                            tGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGaa                                                                                                             <  1:329509/135‑1 (MQ=255)                                               aaCGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAg                                                                                                            >  2:211954/1‑133 (MQ=255)                                               aaCGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAg                                                                                                            <  1:211954/133‑1 (MQ=255)                                                aCGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAgg                                                                                          >  2:341095/1‑150 (MQ=255)                                                  ggCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATc                                                                                             <  2:10604/145‑1 (MQ=255)                                                  ggCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATc                                                                                             >  1:10604/1‑145 (MQ=255)                                                  ggCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATc                                                                                             >  1:724073/1‑145 (MQ=255)                                                  ggCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATc                                                                                             <  2:724073/145‑1 (MQ=255)                                                     tGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGAtt                                                                                     >  2:834931/1‑150 (MQ=255)                                                      gggTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGAttt                                                                                    <  1:1413630/150‑1 (MQ=255)                                                            aGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTa                                                                              >  1:282473/1‑150 (MQ=255)                                                            aGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTa                                                                              >  2:1299941/1‑150 (MQ=255)                                                                        gCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAg                                                                                                            <  1:80196/108‑1 (MQ=255)                                                                        gCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAg                                                                                                            >  2:80196/1‑108 (MQ=255)                                                                           tGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAg                                                               <  1:1299941/150‑1 (MQ=255)                                                                           tGGTGCAATTCAGGCGAAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAg                                                               <  1:1343577/150‑1 (MQ=255)                                                                            ggTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAgg                                                              <  2:967105/150‑1 (MQ=255)                                                                                cAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGAtttt                                                                                   >  1:707937/1‑125 (MQ=255)                                                                                cAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGAtttt                                                                                   <  2:707937/125‑1 (MQ=255)                                                                                 aaTTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCAt                                                                                <  1:1086813/127‑1 (MQ=255)                                                                                 aaTTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCAt                                                                                >  2:1086813/1‑127 (MQ=255)                                                                                            aaaaCAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGa                                              <  1:964850/150‑1 (MQ=255)                                                                                                 aGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTAc                                         <  1:834931/150‑1 (MQ=255)                                                                                                        ccTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCt                                                                         >  2:251316/1‑111 (MQ=255)                                                                                                        ccTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCt                                                                         <  1:251316/111‑1 (MQ=255)                                                                                                                                      ggccggccTTCTGAA‑GCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTc                                                                                  >  2:472036/1‑72 (MQ=37)                                                                                                                                      ggccggccTTCTGAA‑GCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTc                                                                                  <  1:472036/72‑1 (MQ=37)                                                                                                                                        ccGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTCTGaa  <  2:169400/150‑1 (MQ=255)                                                                                                                                        ccGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTCTGaa  <  1:341095/150‑1 (MQ=255)                                                                                                                                          gg‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTc                                                                                  >  2:455577/1‑68 (MQ=255)                                                                                                                                          gg‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTc                                                                                  <  1:455577/68‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                              GGCAGCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGG‑CTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTCTGAA  >  NC_002947/4151839‑4152123

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