Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	152,157	A→T	11.3%	intergenic (‑51/+43)	PP_0143 ← / ← PP_0144	membrane protein/insulinase family metalloprotease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	152,157	0	A	T	11.3%	107.8 / 4.3	44	intergenic (‑51/+43)	PP_0143/PP_0144	membrane protein/insulinase family metalloprotease
Reads supporting (aligned to +/- strand):  ref base A (20/18);  new base T (3/2);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTA  >  NC_002947/152011‑152304                                                                                                                                                   |                                                                                                                                                     gCCGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAg                                                                                                                                                    >  3:16816/1‑149 (MQ=255)  ccGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCaa                                                                                                                                                     <  4:140222/147‑1 (MQ=255)    gATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTAGCTGGCTTGccc                                                                                                                                                 >  1:63073/1‑149 (MQ=255)         cGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                            <  4:16816/149‑1 (MQ=255)            agagGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAgg                                                                                                                                         <  4:131398/149‑1 (MQ=255)             gcggACACCGCCTAGCGCTATGATCCCCGGGTAGCGCTTGAGTAGTTGCATTCTTCGCTGGACTACATGGGCGATTTGGGGCATTGGGGCCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCTTGCCCGCGaaga                                                                                                                                          <  2:63073/145‑3 (MQ=255)                gAAACCGCCGAGCGCAAGGATGCCCGGGTAGCGCTTGAGTTGTTCGCTGATGCGCTGCACGATCTGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCg                                                                                                                                     <  2:286064/149‑1 (MQ=255)                cataCCGCCGAGCGCAAGGCTGCCCGGGGAGCGTTTGAGTAGTTCGATGCGCCGAGGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCg                                                                                                                                     <  6:325196/146‑1 (MQ=255)                     cgccgAGCGCAATGTTGCCCGGGTTGCGCTTGTGTTGTGCGATGATGCGCTGGACTATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTc                                                                                                                                      <  8:411920/143‑1 (MQ=255)                     cgccgAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTc                                                                                                                                      >  7:411920/1‑143 (MQ=255)                     cgccgAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACt                                                                                                                                >  5:393326/1‑149 (MQ=255)                             gcAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAg                                                                                                                        >  2:78937/1‑149 (MQ=255)                             gcAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAg                                                                                                                        <  8:281098/149‑1 (MQ=255)                             gcAATGATGCCCGGGTAGCGCGTGACTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAg                                                                                                                        <  6:257482/149‑1 (MQ=255)                                             aGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTAc                                                                                                       >  5:215966/1‑149 (MQ=255)                                                         gtTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCAAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGAAGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTa                                                                                           >  1:360345/1‑149 (MQ=255)                                                                  cgcggTGGACGTTCTGGGCGAGTCGCGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTg                                                                                  <  6:221548/148‑1 (MQ=39)                                                                     gcTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGc                                                                               >  4:380362/1‑149 (MQ=255)                                                                         gACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCCAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGCTACCGCC‑GATTTACCGCTGCCGTGAAGCACCACTGCGCCGtt                                                                           >  8:366042/1‑149 (MQ=255)                                                                              tatGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGt                                                                      >  2:36210/1‑149 (MQ=255)                                                                                  ggggATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCAt                                                                  >  2:169590/1‑149 (MQ=255)                                                                                  ggggAGTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAgg                                                                       <  6:393326/144‑1 (MQ=255)                                                                                                   ggggCTACGACATAAACGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCaa                                                  <  7:157749/148‑1 (MQ=255)                                                                                                    gggCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGGACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCaa                                                  >  5:372045/1‑147 (MQ=255)                                                                                                                  aaGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTcc                                                                    >  1:28384/1‑115 (MQ=255)                                                                                                                  aaGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTcc                                                                    <  2:28384/115‑1 (MQ=255)                                                                                                                  aaGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAgg                                  >  4:135893/1‑149 (MQ=255)                                                                                                                   aGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGc                                 <  6:372045/149‑1 (MQ=255)                                                                                                                           ggCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGcccc                         <  6:215966/149‑1 (MQ=255)                                                                                                                            gCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCa                        >  7:412114/1‑149 (MQ=255)                                                                                                                                  ccAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGa                     <  8:259695/146‑1 (MQ=255)                                                                                                                                  ccAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGa                     >  7:259695/1‑146 (MQ=255)                                                                                                                                  ccAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGc                   >  3:416109/1‑148 (MQ=255)                                                                                                                                  ccAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGc                   <  3:60623/148‑1 (MQ=255)                                                                                                                                   cAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGc                   >  3:329628/1‑147 (MQ=255)                                                                                                                                    aGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAg                >  6:388562/1‑149 (MQ=255)                                                                                                                                        tctTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGTGCCGTTCAGGTCCATGCCGCCGCCAGCGCAAc                                                 <  5:258887/112‑1 (MQ=255)                                                                                                                                        tctTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAAc                                                 >  6:258887/1‑112 (MQ=255)                                                                                                                                         tttCGCGGGCTTGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCAc                       >  3:410457/2‑137 (MQ=255)                                                                                                                                         tttCGCGGGCTTGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCAc                       <  4:410457/136‑1 (MQ=255)                                                                                                                                         tttCGCGGGCTTGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCAt           >  1:224248/2‑149 (MQ=255)                                                                                                                                          ttCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCAt           >  6:258619/1‑148 (MQ=255)                                                                                                                                           tCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCTGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCCCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCAt           >  1:101543/1‑148 (MQ=255)                                                                                                                                                 gCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGt   <  5:388562/149‑1 (MQ=255)                                                                                                                                                  cAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTa  <  4:416109/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GCCGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCC‑GATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTA  >  NC_002947/152011‑152304

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