Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,803,264	G→A	6.9%	A218V (GCC→GTC)	PP_3361 ←	aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,803,264	0	G	A	6.9%	171.7 / 2.5	58	A218V (GCC→GTC)	PP_3361	aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (22/32);  new base A (2/2);  total (24/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

AGCCGGCCACCAGCGCAGCAGCAGCGCCGATCGGGTCGTCCATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTAGGCTGTACAGTCTGTTCGAACCGGGCGATGAAGCGCTCGGCAACGGCCCGCCGATCAGGGTCCTTGATGCA  >  NC_002947/3803118‑3803403                                                                                                                                                   |                                                                                                                                            aGCCGGCCACCAGCCCAGCAGCAGCGCCGATCGGGTCGTCAATGTCGAGAATGGCGTACGAGCTGGCCATAGCGACCTCGGCAATCAGAACGCTGAGCACCGCAACGCCAAAATCGATCAAGCCAGTCACGCTACGCGGCGTTACGGc                                                                                                                                            >  5:369031/1‑148 (MQ=255)    cGGCCACCAGCGCAGCAGCAGCGCCGATCGGGTCGTCCATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCg                                                                                                                                        <  8:207674/149‑1 (MQ=255)       ccaccaGCGCAGCAGCAGCGCCGATCGGGTCGTCCATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGAcc                                                                                                                                     <  5:276539/149‑1 (MQ=255)           caGCGCAGCAGCAGCGCCGATCGGGTCGTCCATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGc                                                                                                                                            >  4:41260/1‑138 (MQ=255)           caGCGCAGCAGCAGCGCCGATCGGGTCGTCCATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGc                                                                                                                                            <  3:41260/138‑1 (MQ=255)             gcgcAGCAGCAGCGCCGATCGGGTCGTCCATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATg                                                                                                                               <  3:98306/149‑1 (MQ=255)                cagcagcagCGCCGATCGGGTCGTCCATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTc                                                                                                                            <  5:353439/149‑1 (MQ=255)                cagcagcagCGCCGATCGGGTCGTCCATGTCGAGAATGGCGTACGCGCGGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCTTTTCGGCAGCGACCAGGATGtcc                                                                                                                            <  6:371953/149‑3 (MQ=255)                  gcagcagcGCCGATCGGGTCGTCCATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCa                                                                                                                          <  5:305174/149‑1 (MQ=255)                      cagcGCCGATCGGGTCGTCCATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCAt                                                                                                                                          <  1:134525/129‑1 (MQ=255)                      cagcGCCGATCGGGTCGTCCATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCAt                                                                                                                                          >  2:134525/1‑129 (MQ=255)                         cgcCGATCGGGTCGTCCATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTg                                                                                                                   <  3:312612/149‑1 (MQ=255)                           ccGATCGGGTCGTCCATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTgg                                                                                                                  >  7:261270/1‑148 (MQ=255)                             gATCGGGTCGTCCATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCAt                                                                                                               >  3:147086/1‑149 (MQ=255)                              aTCGGGTCGTCCATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATc                                                                                                              <  3:97604/149‑1 (MQ=255)                               tCGGGTCGTCCATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGCAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCa                                                                                                                          >  7:76258/1‑136 (MQ=255)                                   gtcgtcCATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGAc                                                                                                                                      >  7:252016/1‑120 (MQ=255)                                   gtcgtcCATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGAc                                                                                                                                      <  8:252016/120‑1 (MQ=255)                                   gtcgtcCATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTtgt                                                                                                         >  3:129731/1‑149 (MQ=255)                                   gtcgtcCATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTtgt                                                                                                         <  4:294160/149‑1 (MQ=255)                                    tcgtcCATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTtgtg                                                                                                        <  4:56983/149‑1 (MQ=255)                                       tcCATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGAt                                                                                                     >  5:118432/1‑149 (MQ=255)                                         cATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCa                                                                                                   >  2:12574/1‑149 (MQ=255)                                                      gCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTAAGGGGCGTTTCGACATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAg                                                                                      >  3:64439/1‑149 (MQ=255)                                                        gTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAggg                                                                                    <  1:12574/149‑1 (MQ=255)                                                           cgcgCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGaa                                                                                  <  3:7687/148‑1 (MQ=255)                                                           cgcgCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGaa                                                                                  <  4:129731/148‑1 (MQ=255)                                                           cgcgCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAg                                                                                 <  6:26822/149‑1 (MQ=255)                                                              gcTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGtt                                                                                                                    >  6:8228/1‑111 (MQ=255)                                                              gcTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGtt                                                                                                                    <  5:8228/111‑1 (MQ=255)                                                                     aTCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTAgg                                                                       <  1:77440/149‑1 (MQ=255)                                                                      tCGCGACCTCGGCAATCAGAACGGTGTGCAGCGGATCTCCAAAATCGATGATTCCAGTCACGCTACGGGGCGTTTCGGGATCGACGATGATGTTAGTGTAGTTGGCATCGTTGTGGATCACCTGTGCGCGAAGGGAAGCCAGTTTAgg                                                                       <  6:369031/148‑1 (MQ=255)                                                                        gcgACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTAGGCt                                                                     <  4:174425/148‑1 (MQ=255)                                                                             cTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTtgtg                                                                                                        >  4:237496/1‑108 (MQ=255)                                                                             cTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTtgtg                                                                                                        <  3:237496/108‑1 (MQ=255)                                                                                ggCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTAGGCTGTACAGTCt                                                            >  4:181755/1‑149 (MQ=255)                                                                                        gAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTAGGCTGTACAGTCTGTTCGAAc                                                    <  3:181755/149‑1 (MQ=255)                                                                                         aaCGGTGTGCACCTGATCTCCAAAATCGATGACGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGGCGTTGGCATCGTTGTGGAGCACCTGGGCGCGAAGGGAAGCCAGTTTAGGCTGTACAGTCTGTTCGAAcc                                                   <  8:261270/149‑1 (MQ=255)                                                                                            ggTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTAGGCTGTACAGTCTGTTCGAACCggg                                                <  1:135902/149‑1 (MQ=255)                                                                                            ggTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCAGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTAGGCTGTACAGTCTGTTCGAACCggg                                                <  1:176717/149‑1 (MQ=255)                                                                                             gtgtGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTa                                                                         <  3:132802/123‑1 (MQ=255)                                                                                             gtgtGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTa                                                                         >  4:132802/1‑123 (MQ=255)                                                                                                   aCCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCTTCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTAGGCTGTACAGTCTGTTCGAACCGGGCGATg                                           >  1:271683/1‑147 (MQ=255)                                                                                                   aCCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTAGGCTGTACAGTCTGTTCGAACCGGGCGATGaa                                         <  4:147086/149‑1 (MQ=255)                                                                                                           gCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTAGGCTGTACAGTCTGTTCGAACCGGGCGATGAAGCGCTCgg                                 <  6:118432/149‑1 (MQ=255)                                                                                                           gCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCAAGTTTAGGCGTTACAGTCTGTTCGAACCGGGCGATGAACCGCTAgg                                 >  4:270625/1‑149 (MQ=255)                                                                                                           gCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGACATCGACCAGGATGTTCCAGTCGTTGGCATCGTtgt                                                                                                         >  8:6352/1‑77 (MQ=255)                                                                                                           gCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGACATCGACCAGGATGTTCCAGTCGTTGGCATCGTtgt                                                                                                         <  7:6352/77‑1 (MQ=255)                                                                                                                aaTCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTAGGCTGTACAGTCTGTTCGAACCGGGCGATGAAGCGCTCGGCAACg                            >  5:61039/1‑149 (MQ=255)                                                                                                                     aTGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTAGGCTGTACAGTCTGTTCGAACCGGGCGATGAAGCGCTCGGCAACGGCccg                       >  6:264879/1‑149 (MQ=255)                                                                                                                           ccAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTAGTTGGCATAGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTAGGCTGTACAGTCTGTTCGAACCGGGCGATGAAGCGCTCGGCAACGGCCCGCCGAt                  >  1:262656/1‑148 (MQ=255)                                                                                                                             aGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTAGGCTGTACAGTCTGTTCGAACCGGGCGATGAAGCGCTCGGCAACgg                           <  7:186308/137‑1 (MQ=255)                                                                                                                             aGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTAGGCTGTACAGTCTGTTCGAACCGGGCGATGAAGCGCTCGGCAACGGCCCGCCGATCAg               <  2:271683/149‑1 (MQ=255)                                                                                                                             aGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTAGGATGTACAGTCTGTTCGAAACGTGCGATGAATCGCTCGGCAAAgg                           >  8:186308/1‑137 (MQ=255)                                                                                                                                   gCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTAGGCTGTACAGTCTGTTCGAACCGGGCGATGAAGCGCTCGGCAACGGCCCGCCGATCAGGGTCCt         >  6:363917/1‑149 (MQ=255)                                                                                                                                          ggCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTAGGCTGTACAGTCTGTTCGAACCGGGCGATGAAGCGCTCGGCAACgg                           >  4:292596/1‑124 (MQ=255)                                                                                                                                          ggCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTAGGCTGTACAGTCTGTTCGAACCGGGCGATGAAGCGCTCGGCAACgg                           <  3:292596/124‑1 (MQ=255)                                                                                                                                          ggCGTTTCGACATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTAGGCTGTACAGGCTGTTCGAACCGGGCGATGAAGCGCTCGGCAACGGCCCGCCGATCAGGGTCCTTGAtgca  <  4:64439/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                            AGCCGGCCACCAGCGCAGCAGCAGCGCCGATCGGGTCGTCCATGTCGAGAATGGCGTACGCGCTGGCGATCGCGACCTCGGCAATCAGAACGGTGTGCACCGCATCGCCAAAATCGATGAAGCCAGTCACGCTACGGGGCGTTTCGGCATCGACCAGGATGTTCCAGTCGTTGGCATCGTTGTGGATCACCTGGGCGCGAAGGGAAGCCAGTTTAGGCTGTACAGTCTGTTCGAACCGGGCGATGAAGCGCTCGGCAACGGCCCGCCGATCAGGGTCCTTGATGCA  >  NC_002947/3803118‑3803403

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