Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,442,673	C→G	12.5%	intergenic (‑99/+39)	nicP‑II ← / ← nicT	porin‑like protein/metabolite transport protein NicT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,442,673	0	C	G	12.5%	130.4 / 10.8	48	intergenic (‑99/+39)	nicP‑II/nicT	porin‑like protein/metabolite transport protein NicT
Reads supporting (aligned to +/- strand):  ref base C (23/19);  new base G (4/2);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

AGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCGGT  >  NC_002947/4442531‑4442819                                                                                                                                                 |                                                                                                                                                   aGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGc                                                                                                                                             <  2:125830/150‑1 (MQ=255)      tGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGTGCCCGAACAGCCAGGCCAACCGGC‑AGCGTGGTTG‑ATCTAGAGAGGTCGGCAGGCGGGCCCCGGCCCGcctgcc                                                                                                                                        <  1:21247/150‑1 (MQ=255)               agaaAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGc                                                                                                                                                 >  1:484976/1‑132 (MQ=255)               agaaAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGc                                                                                                                                                 <  2:484976/132‑1 (MQ=255)                 aaaaaCGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCg                                                                                                                             >  1:920127/1‑150 (MQ=255)                     aCGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGaa                                                                                                                         <  2:969847/150‑1 (MQ=255)                      cgcgGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAg                                                                                                                        <  1:707158/150‑1 (MQ=255)                           cgcAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACCGATTCCAGTGCCATCGTCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTg                                                                                                                   <  2:870544/150‑1 (MQ=255)                            gcAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTgg                                                                                                                  <  2:209975/150‑1 (MQ=255)                                  aaaCCTTGCATAGGGCTACTCCTGGATTTGGGCAAAAGCAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCt                                                                                                            >  1:694370/1‑150 (MQ=255)                                      cTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCg                                                                                                        >  2:1026898/1‑150 (MQ=255)                                           aTAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGa                                                                                                   >  2:740272/1‑150 (MQ=255)                                           aTAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGa                                                                                                   >  1:450493/1‑150 (MQ=255)                                            tAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGAc                                                                                                  <  1:740272/150‑1 (MQ=255)                                                 cTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTg                                                                                                                   >  2:814115/1‑128 (MQ=255)                                                 cTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTg                                                                                                                   <  1:814115/128‑1 (MQ=255)                                                     tCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGc                                                                                         >  2:5773/1‑150 (MQ=255)                                                        tGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCCGGCCCCGGCCCGCCTGCCCTGCCGGTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCgg                                                                                      >  1:23931/1‑150 (MQ=255)                                                            tttGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGACTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAg                                                                                  >  2:337886/1‑150 (MQ=255)                                                             tctggCAACAGCAACCCGGCGCGGCGGCACTGCATGGAC‑TCCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAgg                                                                                 <  1:10183/147‑1 (MQ=255)                                                                 gCAACAGCAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCAGGGCCCCTACCGGTTGACCAGCTTGGCCGGCAAGGcaga                                                                             >  2:10183/1‑147 (MQ=255)                                                                   aaCAGCAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGa                                                                           >  2:447063/1‑150 (MQ=255)                                                                   aaCAGCAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGa                                                                           >  1:610267/1‑150 (MQ=255)                                                                        cAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAgg                                                                      >  1:1125543/1‑150 (MQ=255)                                                                           cccGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAag                                                                   <  1:337886/150‑1 (MQ=255)                                                                            ccGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCCAGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAagc                                                                  >  1:618862/1‑150 (MQ=255)                                                                              ggcgcggcgGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAagcag                                                                <  2:921800/150‑1 (MQ=255)                                                                                  cggcggCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGa                                                            <  1:541341/150‑1 (MQ=255)                                                                                         aCTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCtt                                                                                            >  1:1125963/1‑111 (MQ=255)                                                                                         aCTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCtt                                                                                            <  2:1125963/111‑1 (MQ=255)                                                                                           tGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCAATACCTACAGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAAAAGCTGAGGAACACGc                                                   >  2:261982/1‑150 (MQ=255)                                                                                              aTGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGc                                                <  2:333794/150‑1 (MQ=255)                                                                                               tGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGcc                                               <  2:920127/150‑1 (MQ=255)                                                                                                      ccGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCaaa                                         >  1:149111/1‑150 (MQ=255)                                                                                                             gCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCcac                                  >  1:802469/1‑150 (MQ=255)                                                                                                                   ttGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCCGGCAGCCGGCAAACACCCAcagcgc                             <  1:297890/150‑1 (MQ=255)                                                                                                                    tGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCAcagcgca                            <  1:5773/150‑1 (MQ=255)                                                                                                                       aTCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCAcagcgcagcg                         >  2:171992/1‑150 (MQ=255)                                                                                                                                 gTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGc               >  2:1084886/1‑150 (MQ=255)                                                                                                                                 gTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGc               >  1:906800/1‑150 (MQ=255)                                                                                                                                       aGGCGGGCCGGGGCCCTCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAacac                                      >  2:839736/1‑121 (MQ=255)                                                                                                                                       aGGCGGGCCGGGGCCCTCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAacac                                      <  1:839736/121‑1 (MQ=255)                                                                                                                                       aGGCGGGCCGGGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCAcagcgcagc                          <  2:1010458/133‑1 (MQ=255)                                                                                                                                       aGGCGGGCCGGGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCAcagcgcagc                          >  1:1010458/1‑133 (MQ=255)                                                                                                                                       aGGCGGGCCGGGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGa         >  2:1100846/1‑150 (MQ=255)                                                                                                                                       aGGCGGGCCGGGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGa         >  2:1028752/1‑150 (MQ=255)                                                                                                                                       aGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATc                                                       <  2:675899/104‑1 (MQ=255)                                                                                                                                       aGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATc                                                       >  1:675899/1‑104 (MQ=255)                                                                                                                                              ccccGGCCCGCCGGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCGGt  <  2:610267/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   AGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCC‑ACCGGCGTGCG‑GGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCGGT  >  NC_002947/4442531‑4442819

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