Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	103,821	A→T	8.6%	intergenic (+119/‑21)	PP_0098 → / → PP_0099	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	103,821	0	A	T	8.6%	151.5 / 4.9	58	intergenic (+119/‑21)	PP_0098/PP_0099	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/27);  new base T (3/2);  total (29/29)
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.68e-01

GTTTGGTGCAGTGCTGCAGCCCTAAGCGCTGTCTGTCTGGGATGGGCGGCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGCCCCTGGTCATTCGGTGGCGGAAGATGATGGCTGGTACCAGGAAGTGCCGCAAACCCAGGGTAC  >  NC_002947/103678‑103966                                                                                                                                                 |                                                                                                                                                  gttTGGTGCAGTGCTGCAGCCCTAAGCGCTGTCTGTCTGGGATGGGCGTCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TTAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATAc                                                                                                                                               <  5:254252/148‑1 (MQ=255)  tttGGTGCAGTGCTGCAGCCCTAAGCGCTGTCTGTCTGGGATGGGCGGCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACa                                                                                                                                              <  5:234393/148‑1 (MQ=255)     ggTGCAGTGCTGCAGCCCTAAGCGCTGTCTGTCTGGGATGGGCGGCATTCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTa                                                                                                                                           <  5:3189/148‑1 (MQ=255)      gTGCAGTGCTGCAGCCCTAAGCGCTGTCTGTCTGGGATGGGCGGCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATTAATATACAGTa                                                                                                                                           <  1:16300/147‑1 (MQ=255)      gTGCAGTGCTGCAGCCCTAAGCGCTGTCTGTCTGGGATGGGCGGCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATTAATATACAGTa                                                                                                                                           >  2:16300/1‑147 (MQ=255)       tGCAGTGCTGCAGCCCTAAGCGCTGTCTGTCTGGGATGGGCGGCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTc                                                                                                                                        <  6:254036/149‑1 (MQ=255)        gCAGTGCTGCAGCCCTAAGCGCTGTCTGTCTGGGATGGGCGGCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTAtt                                                                                                                                         >  2:72967/1‑147 (MQ=255)           gTGCTGCAGCCCTAAGCGCTGTCTGTCTGGGATGGGCGGCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGt                                                                                                                                    <  1:229383/149‑1 (MQ=255)                gcAGCCCTAAGCGCTGTCTGTCTGGGATGGGCGGCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGc                                                                                                                                >  5:27474/1‑148 (MQ=255)                 cAGCCCTAAGCGCTGTCTGTCTGGGATGGGCGGCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCa                                                                                                                              >  8:87990/1‑149 (MQ=255)                    cccTAAGCGCTGTCTGTCTGGGATGGGCGGCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGTTACTGTATATTAATACAGTATTCAGGt                                                                                                                                    <  4:98729/140‑1 (MQ=255)                    cccTAAGCGCTGTCTGTCTGGGATGGGCGGCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGTTACTGTATATTAATACAGTATTCAGGt                                                                                                                                    >  3:98729/1‑140 (MQ=255)                     ccTAAGCGCTGTCTGTCTGGGATGGGCGGCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGAt                                                                                                                          <  7:105502/149‑1 (MQ=255)                         aGCGCTGTCTGTCTGGGATGGGCGGCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATtc                                                                                                                        <  7:240619/147‑1 (MQ=255)                          gcgcTGTCTGTCTGGGATGTGCGGCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGTATTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATtcctc                                                                                                                     <  8:281857/149‑1 (MQ=255)                           cgcTGTCTGTCTGGGATGGGCGGCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATtcctcc                                                                                                                    <  1:72967/149‑1 (MQ=255)                                  tgtctgGGATGGGCGGCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGt                                                                                                                                    >  4:7111/1‑126 (MQ=255)                                  tgtctgGGATGGGCGGCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGt                                                                                                                                    <  3:7111/126‑1 (MQ=255)                                  tgtctgGGATGGGCGGCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATtcctc                                                                                                                     >  2:76635/1‑141 (MQ=255)                                  tgtctgGGATGGGCGGCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCg                                                                                                             >  2:44196/1‑149 (MQ=255)                                  tgtctgGGATGGGCGGCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTCCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATtcctc                                                                                                                     <  1:76635/141‑1 (MQ=255)                                    tctgGGATGGGCGGCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGc                                                                                                            >  2:69685/1‑148 (MQ=255)                                        ggATCGGCGTCAATCGAGACTGGCAGGGTATGTTCGCTGGCATGGCC‑TGAAGTAGTCAGAACAGATGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAg                                                                                                       <  4:255775/149‑1 (MQ=255)                                               cggcAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTg                                                                                                >  4:200871/1‑149 (MQ=255)                                                ggcAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATtcctcc                                                                                                                    >  4:30534/1‑128 (MQ=255)                                                ggcAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATtcctcc                                                                                                                    <  3:30534/128‑1 (MQ=255)                                                ggcAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTgg                                                                                               <  6:27474/149‑1 (MQ=255)                                                    aTCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCCTTGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCAc                                                                                            >  1:230222/1‑149 (MQ=255)                                                                aGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCTCAAAAGGGCCGTGGCACGGCGCACAATCca                                                                               <  3:323930/149‑1 (MQ=255)                                                                aGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCca                                                                               <  2:230222/149‑1 (MQ=255)                                                                 gCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCcac                                                                              <  1:69685/149‑1 (MQ=255)                                                                             tGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATc                                                                                 <  2:103894/134‑1 (MQ=255)                                                                             tGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATc                                                                                 >  1:103894/1‑134 (MQ=255)                                                                               gCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGt                                                                                                                                    <  6:336133/81‑1 (MQ=255)                                                                               gCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGt                                                                                                                                    >  5:336133/1‑81 (MQ=255)                                                                                  tGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGCCCCt                                                             >  3:193721/1‑149 (MQ=255)                                                                                            gtagtCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGCCCCTGGTCATTCg                                                    <  4:193721/149‑1 (MQ=255)                                                                                              agtCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGccc                                                               >  6:142414/1‑136 (MQ=255)                                                                                              agtCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGccc                                                               <  5:142414/136‑1 (MQ=255)                                                                                                tCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGCCCCTGGTCATTCGGTgg                                                >  8:176533/1‑149 (MQ=255)                                                                                                tCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACAGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGCCCCTGGTCATTCGGTgg                                                >  8:44711/1‑149 (MQ=255)                                                                                                  aGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGCCCCTGGTCATTCGGTGGCg                                              >  8:160595/1‑149 (MQ=255)                                                                                                         agGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGCCCCTGGTCATTCGGTGGCGGAAgatg                                       <  1:44196/149‑1 (MQ=255)                                                                                                               tGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGCCCCTGGTCATTCGGTGGCGGAAGATGATGGCt                                 >  2:57819/1‑149 (MQ=255)                                                                                                                    aGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGCCCCTGGTCATTCGGTGGCGGAAGATGATGGCTGGTAc                            >  8:78276/1‑149 (MQ=255)                                                                                                                     gATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGCCCCTGGTCATTCGGTGGCGGAAGATGATGGCTGGTAcc                           <  7:87990/149‑1 (MQ=255)                                                                                                                            gtgtGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGCCCCTGGTCATTCGGTGGCGGAAGATGATGGCTGGTACCAGGAAGt                    <  3:200871/149‑1 (MQ=255)                                                                                                                                 ctgaaTACTGTATTAATATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGCCCCTGGTCATTCGGTGGCGGAAGATGATGGCTgg                               <  3:174175/129‑1 (MQ=255)                                                                                                                                 ctgaaTACTGTATTAATATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGCCCCTGGTCATTCGGTGGCGGAAGATGATGGCTgg                               >  4:174175/5‑133 (MQ=255)                                                                                                                                 cACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGCCCCTGGTGATGCGGTGGCGGAAGGTGATGGCTGGTACCAGGAAGTGCCg                >  6:300903/1‑148 (MQ=255)                                                                                                                                 cACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGCCCCTGGTCATTCGGTGGCGGAAGATGATGGCTGGTACCAGGAAGTGCCg                <  7:78276/148‑1 (MQ=255)                                                                                                                                  aCGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGCCCCTGGTCATTCGGTGGCGGAAGATGATGGCTGGTACCAGGAAGTGCCGc               >  8:24466/1‑148 (MQ=255)                                                                                                                                   gaaTACTGTATTAATATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGAACAATACACACAACCGCTTTGCCCCTGGTCATTCGGTGGCGGAAGATGATGGCTGGTACCAGGAAGTGCCGCa              >  3:320308/3‑148 (MQ=255)                                                                                                                                         tGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGCCCCTGGTCATTCGGTGGCGGAAgatgat                                     >  6:311330/1‑119 (MQ=255)                                                                                                                                         tGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGCCCCTGGTCATTCGGTGGCGGAAgatgat                                     <  5:311330/119‑1 (MQ=255)                                                                                                                                           tataTTAATACAGTATTCAGGTGAGCCATGATTCCCCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGCCCCTGGTCATTCGGTGGCGGAAGATGATGGCTGGTACCAGGAAGTGCCGCAAACCCAggg     >  2:279186/1‑149 (MQ=255)                                                                                                                                            atatTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGCCCCTGGTCATTCGGTGGCGGAAGATGATGGCTGGTACCAGGAAGTGCCGCAAACCCAggg     >  1:143237/1‑148 (MQ=255)                                                                                                                                              atTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGCCCCTGGTCATTCGGTGGCGGAAGATGATGGCTGGTACCAGGAAGTGCCGCAAACCCAGGGTAc  >  4:110780/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                  GTTTGGTGCAGTGCTGCAGCCCTAAGCGCTGTCTGTCTGGGATGGGCGGCAATCGAGCCAGGCAGCGTATGTGCGCTGGCATGGCC‑TGAAGTAGTCAGACCAGAGGTGCTGGTAAGATCCCGGTGTGCACGATACTGTATATTAATACAGTATTCAGGTGAGCCATGATTCCTCCAGTACCGCAAAAGGGCCGTGGCACGGCGCACAATCCACACAACCGCTTTGCCCCTGGTCATTCGGTGGCGGAAGATGATGGCTGGTACCAGGAAGTGCCGCAAACCCAGGGTAC  >  NC_002947/103678‑103966

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