Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,295,827	G→A	12.7%	intergenic (+73/+60)	PP_5496 → / ← PP_2023	hypothetical protein/glutathione S‑transferase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,295,827	0	G	A	12.7%	116.2 / 8.2	48	intergenic (+73/+60)	PP_5496/PP_2023	hypothetical protein/glutathione S‑transferase family protein
Reads supporting (aligned to +/- strand):  ref base G (20/21);  new base A (3/3);  total (23/25)
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.37e-01

CACAGCTACACCGAGCAGTACCTGCGCCAGAGCGCCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAAGCCAGCCACCCCTCCAGCCAGCTCCGCACACGCGGATAAGCCGCCTCGGCGAAC  >  NC_002947/2295680‑2295974                                                                                                                                                     |                                                                                                                                                    cacaGCTACACCGAGCAGTACCTGCGCCAGAGCGCCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCgg                                                                                                                                                    >  5:243280/1‑149 (MQ=255)  acaGCTACACCGAGCAGTACCTGCGCCAGAGCGCCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGt                                                                                                                                                   <  6:243280/149‑1 (MQ=255)      cTACACCGAGCAGTACCTGCGCCAGAGCGCCAATTTCCATGCGGCCGTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGccc                                                                                                                                               <  4:269180/149‑1 (MQ=255)            cGAGCAGTACCTGCGCCAGAGCGCCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTcg                                                                                                                                         <  4:282833/149‑1 (MQ=255)                  gTACCTGCGCCAGAGCGCCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCa                                                                                                                                   >  3:107985/1‑149 (MQ=255)                             agagCGCCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGGGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTccc                                                                                                                        <  7:90466/149‑1 (MQ=255)                               agCGCCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCcac                                                                                                                      <  2:353707/149‑1 (MQ=255)                                gccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttaccctgtgggagcgggctTGCCCGCGAACACCGGCAGAGCCGGTGCCATGCACCGCAGTGCTTGCTTcgcg                                                                                                                                       <  7:166011‑M1/53‑1 (MQ=255)                                gccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttaccctgtgggagcgggctTGCCCGCGAACACCGGCAGAGCCGGTGCCATGCACCGCAGTGCTTGCTTcgcg                                                                                                                                       >  8:166011‑M1/80‑132 (MQ=255)                                 cgcCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCacac                                                                                                                    <  2:307554/149‑1 (MQ=255)                                 cgcCAAATTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCacac                                                                                                                    <  8:172319/149‑1 (MQ=255)                                  gcCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACg                                                                                                                   >  5:292263/1‑149 (MQ=255)                                  cgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcacctcaccctgtgggagcgggctTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCgg                                                                                                                                      >  5:124007‑M1/78‑131 (MQ=255)                                  cgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcacctcaccctgtgggagcgggctTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCgg                                                                                                                                      <  6:124007‑M1/54‑1 (MQ=255)                                            aTGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGATGGGGCCGGGCATGCCCCGCACCCCCGGCAAGGCCGGGCCCTACCCCCCCGTTCCCCGCTTGCCGGGCATCCCCCGGCCCCCCCGGCCCCAtgt                                                                                                          >  1:107633/1‑149 (MQ=255)                                            aTGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCAtgtg                                                                                                         >  1:306082/1‑149 (MQ=255)                                             tGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGGATGCCCACTCCCACAAGGACCCCTGTGc                                                                                                        >  8:372759/1‑149 (MQ=255)                                                ggCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCAGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGAc                                                                                                                <  8:301658/138‑1 (MQ=255)                                                ggCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCAGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGAc                                                                                                                >  7:301658/1‑138 (MQ=255)                                                  cccTCAGGGGCAAGCACCACCACTGTGCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCTAACACCGGCAGCGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGccc                                                                                                   <  2:370195/149‑1 (MQ=37)                                                      gAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCGTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGAt                                                                                                <  6:292263/148‑1 (MQ=255)                                                       aGTGGCAAGCACCACCAGTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGTTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGAt                                                                                                <  7:372759/147‑1 (MQ=255)                                                       aGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGt                                                                                              >  1:265008/1‑149 (MQ=255)                                                        gTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGtt                                                                                             <  2:265008/149‑1 (MQ=255)                                                           gCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTa                                                                                           >  4:349038/1‑148 (MQ=255)                                                           gCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTa                                                                                           <  5:195438/148‑1 (MQ=255)                                                           gCAAGAACCACCACTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTaa                                                                                          <  8:4799/149‑1 (MQ=255)                                                               gCACCCCCATTTTTCTGGCTTGTGTGTGGGCCTTG‑TGGGTGCGGGCTTGCCCGCTAAAACCGGCATAGCGGTTGCAATCCACCGCGGTGCCCGCTTCGCGGGAATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAg                                                                                      <  2:157203/149‑1 (MQ=255)                                                               gCACCACCACTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGtt                                                                                             >  7:2285/1‑142 (MQ=255)                                                                   caccacTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCCTGTGCAGCCCGCTGTGTAATCAGCggc                                                                                  >  4:258318/1‑148 (MQ=255)                                                                    accacTGATCTGGGTTGGGGGAGTGCCTTG‑GGGGCGCGGGCATCCCCGCTAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGAc                                                                                 <  3:196464/149‑1 (MQ=37)                                                                           aTCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCCCATAGCCGGCCTGCCCGCGCCCCCACGGACCCCAGTGCAGCCCCATGTTTAATCAGCGGACCATGAtt                                                                          >  6:38934/1‑149 (MQ=37)                                                                               gtcTTGTTGGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCTAACACCGGCAGAGCCGGTGCAATCCACCGCTGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGcc                                                                       <  7:149213/146‑1 (MQ=37)                                                                                  ttGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTg                                                                    >  1:6226/1‑148 (MQ=255)                                                                                           gTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCCCCCGGACCCATGTGCAGCCCGATGTTTAATAAGCCGACCATGATTTCCTTGAACAGCTCCg                                                          >  4:254909/1‑149 (MQ=37)                                                                                               ttG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAagc                                                      >  1:110264/1‑149 (MQ=255)                                                                                               ttG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAagc                                                      >  6:394847/1‑149 (MQ=255)                                                                                                 tgggggCGCGGGCATGCCCGCTAACACCGGCAGAGCGGGTGCCATCCACGGGGGTGCCCGCTTGGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAagc                                                      <  5:38934/145‑1 (MQ=38)                                                                                                   tgGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAagcca                                                    <  4:107985/148‑1 (MQ=255)                                                                                                    gggAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCAACCACCCCGACCCCTGTGCAGCCCCATGTTTAATAAGCGGACCATTATTTCCTTGAACACCAACGAagccagc                                                  >  4:233906/1‑149 (MQ=25)                                                                                                                 CCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAAc                                                                 >  7:5641‑M1/1‑121 (MQ=255)                                                                                                                 CCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAAc                                                                 <  8:5641‑M1/121‑1 (MQ=255)                                                                                                                          acCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAAGCCAGCCACCCCTCCAGCCAGCTCCGca                            <  3:233906/149‑1 (MQ=255)                                                                                                                            cGGCAGAGCCAGTGCCATCCACCGCGGTGTCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAAGCCAGCCACCCCTCCAGCCAGCTCCGcaca                          >  6:345442/1‑149 (MQ=255)                                                                                                                               cAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAAGCCAGCCACCCCTCCAGCCAGCTCCGCACAcgc                       <  2:6226/149‑1 (MQ=255)                                                                                                                                           cATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAAGCCAGCCACCCCTCCAGCCAGCTCCGCACACGCGGATAAgccgc            >  1:229474/1‑148 (MQ=255)                                                                                                                                              ccACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAAGCCAGCCACCCCTCCAGCCAGCTCCGCACACGCGGATAAGCCGCCTCg        >  5:147490/1‑149 (MQ=255)                                                                                                                                                aCCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCg                                                                                    >  6:123085/1‑71 (MQ=255)                                                                                                                                                aCCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCg                                                                                    <  5:123085/71‑1 (MQ=255)                                                                                                                                                   gcgGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGAGCATGATTTCCTTGAACAGCTCCGAAGCCAGCCACCCCTCCAGCCAGCTCCGCAGACGCTGATAAGCCGCCTCGGCGaa   <  3:349038/149‑1 (MQ=255)                                                                                                                                                    cgGGGCCCGCTTCGCGGGCATGCCCGCGCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAAGCCAGCCACCCCTCCAGCCAGCTCCGCACACGCGGATAAGCCGCCTCGGCGAAc  <  3:258318/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CACAGCTACACCGAGCAGTACCTGCGCCAGAGCGCCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTG‑TGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAAGCCAGCCACCCCTCCAGCCAGCTCCGCACACGCGGATAAGCCGCCTCGGCGAAC  >  NC_002947/2295680‑2295974

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