Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	718,461	G→A	10.5%	intergenic (+31/+46)	thiO → / ← PP_0613	FAD‑dependent glycine/D‑amino acid oxidase/amidase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	718,461	0	G	A	10.5%	116.8 / 3.3	48	intergenic (+31/+46)	thiO/PP_0613	FAD‑dependent glycine/D‑amino acid oxidase/amidase family protein
Reads supporting (aligned to +/- strand):  ref base G (19/23);  new base A (2/3);  total (22/26)
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.61e-01

TTATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACC  >  NC_002947/718313‑718609                                                                                                                                                     |                                                                                                                                                      ttATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTAcg                                                                                                                                                       <  6:3501/149‑1 (MQ=255)        aaTGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCgg                                                                                                                                                <  6:311693/149‑1 (MQ=255)         aTGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCggg                                                                                                                                               <  8:296875/149‑1 (MQ=255)          tGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGt                                                                                                                                              <  5:32681/149‑1 (MQ=255)              cTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCGCTCCTAAGCGAGCCGGGTT‑Ac                                                                                                                                         >  2:71353/1‑149 (MQ=255)              cTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCAGAGCCGATTATCGACCCGGCGCCGTGTGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCGCGTAGGAGCGGGTTT‑Ac                                                                                                                                         <  6:78070/149‑1 (MQ=38)               tGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCGACGCGAGCGGGTTT‑Acc                                                                                                                                        >  5:274252/1‑149 (MQ=255)                       gcgcCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCGCGTAGGAGCGGGTTT‑ACCcgcg                                                                                                                                    >  5:102853/1‑145 (MQ=38)                       gcgcCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCGCGTAGGAGCGGGTTT‑ACCcgcg                                                                                                                                    <  6:102853/145‑1 (MQ=38)                       gcgcCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTT‑ACCCGCGaaa                                                                                                                                 <  6:274252/148‑1 (MQ=255)                           cGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCgg                                                                                                                                                <  1:41661/130‑1 (MQ=255)                           cGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCgg                                                                                                                                                >  2:41661/1‑130 (MQ=255)                           cGGCTTCGTGCATGCTGTTTACCGATTTGCTGACAGGGGCCGAGCGTATTATCGACCCGGCGCCGTATGCGCCAGATGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGcc                                                                                                                            <  8:241886/149‑1 (MQ=39)                                  gTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTTAACCCCCTCCAACGCGAGCGGGTTT‑ACCCCCGAAAGGGCCGCAGaaa                                                                                                                     >  2:218522/1‑148 (MQ=255)                                   tGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGc                                                                                                                    <  6:86703/149‑1 (MQ=255)                                   tGCAAGCTGTTTAACGTTTTGCTGACAGGGGCCGAGCCTTTTATCGACCCGGCGCCGTATGCGCCATAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCGGCTCCTACGCTAGCGGTTTT‑ACCCGCGAAAGGGCCGCAGAAGc                                                                                                                    <  6:93894/149‑1 (MQ=255)                                     ccagagttttACCGATTTGCTTACAGGGGCCGTGCCGATTATCGACCCGGCGCCGTTTGCGCCTGAGGGGGGTTTGGGGTGTGGACCTCTTCGCGGGTAACCCCGCTCCTACGCGAGCGGTTTT‑ACCCGCGAAAGGGCCGCAGAAGCCa                                                                                                                  <  6:236444/142‑1 (MQ=255)                                      cAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCCCAAAAACCCCAGCCGGTTT‑TACCCCGAAAAGGCCCCAGAAACCa                                                                                                                  >  2:289420/1‑148 (MQ=255)                                                  gATTTGCTGACCGGGGCCGAGCCGATTATCGACCCGGCGCCGTAGGCGCCAGAAGGGCGTTTGGGCTGCGCACCTCTTCGCGGGAAAACCCGCTCCTACGCGCGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAg                                                                                                     <  1:218522/149‑1 (MQ=255)                                                        cTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTTAAACCGCTCCTACGCGAGCGGGTTT‑TACCCCGAAAAGGCCGCCCAAGCCAAAACAGAATTAAGCCAGt                                                                                                >  4:303451/1‑148 (MQ=255)                                                           aCAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGCACCTCTTCGCGGGTAAACCCGATCCTACGCGAGCGGGTTTCACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAAt                                                                                             <  1:71353/149‑1 (MQ=255)                                                           aCAGGGGCCGACCCTTTTATGGACCCGGCGCCGTTTGCGCCATAAGGGGGTTTGGGCTGAGGACCTCTGCGCGGTTAAACCCGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATc                                                                                            <  4:23575/149‑1 (MQ=255)                                                             cggggCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGc                                                                                                    <  7:55950/138‑1 (MQ=255)                                                             cgagctcgagccgaTTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGTACCTATTCGCGGGTAAACCCGCTCCTACGCGCGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCg                                                                                           <  7:12645/142‑1 (MQ=39)                                                             aGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCGCTCCAACGCCAGCGGGTTT‑AACCGCGAAAGGGCCCCAGAAGCCAATACACAATTCAGCCAgaacgga                                                                                           >  8:55950/1‑142 (MQ=255)                                                                                 aCCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGc                                                                             <  3:193215/142‑1 (MQ=255)                                                                                 aCCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCGACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCCATACAGATTTAAGGCAGTAATCGGTGGTTGCAGGTGCCggggcg                                                                      >  4:83393/1‑145 (MQ=39)                                                                                 aCCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCGCGCCGACGCGAGCGGGGTT‑AACCGCGCAAGGGGCGCAGAAGCCAATAAAGAATTGAGCCAGTAATAGGTGCTTTCAGGTgg                                                                             >  4:193215/1‑141 (MQ=255)                                                                                        gcCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAAtg                                                               <  1:97830/149‑1 (MQ=255)                                                                                          cGTATGCGCCTGATGGGGGTTTGGGCTTCGGACCGCTTCGCGGGTAAACCCGCTCCTACGCGCGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAAtgct                                                             <  3:94181/149‑1 (MQ=38)                                                                                                  ccAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTTAAAACCCTCCTACGCGAGCGGGTTT‑ACCCGCGAAATGGCCGCa                                                                                                                         >  1:276086/1‑81 (MQ=38)                                                                                                       aGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCCCTACTAAGCGACCGGGTTT‑ACCCCCCAAAGGGCCGCAGAAACCAATACAGATTTAAGCCAAAAATAGGTGCTTTCAGGTGCCGGCTACGATAATGCAGCTCATAAACCtt                                                >  6:248008/1‑149 (MQ=255)                                                                                                         ggCGTTTGGGATGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTg                                               <  1:289420/148‑1 (MQ=255)                                                                                                          gcttttGGGCTGAGCACCTCTTCGCGGTTAAACCCGCTCCTACGCCTGCGGTTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATgctgct                                                          <  7:260321/133‑1 (MQ=255)                                                                                                          gCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCGACGCGCGCGGGTTT‑ACCCGCGGAAGGGCCGCCGAAGCCAATACAGATTTCAGCCAGGAATCGGTGCTTGCCGGTGCCGGGGCGGATTATGCTGCTCATACACCTTGcc                                             >  5:304083/1‑148 (MQ=38)                                                                                                                          ccTCTTGGCGGGTAAACCCGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGAAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGc                             <  5:248008/149‑1 (MQ=255)                                                                                                                           ctctTCGCGGGTAAACCCGCTCCTACGCGCGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTGAGCCAGTAATCGGTGCTTGCAGGTGGCGGGGGGGATAAtg                                                               >  2:87252/1‑114 (MQ=37)                                                                                                                           ctctTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCt                            >  7:241784/1‑149 (MQ=255)                                                                                                                           ctctTCGCGGGTAAACCCCCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGGAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCt                            >  6:90116/1‑149 (MQ=255)                                                                                                                             tttCGCGGGTAAACCCGCTCGCGTAGGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGcc                                                                            <  1:101246/98‑1 (MQ=21)                                                                                                                             tttCGCGGGTAAACCCGCTCGCGTAGGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGcc                                                                            >  2:101246/2‑99 (MQ=21)                                                                                                                                   gggTAAACCCGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCa                    <  8:171975/149‑1 (MQ=255)                                                                                                                                     gTAAACCCGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCAc                                 <  3:214960/134‑1 (MQ=255)                                                                                                                                     gTAAACCCGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCAc                                 >  4:214960/1‑134 (MQ=255)                                                                                                                                         aCCCGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAgcgc              >  5:17273/1‑149 (MQ=255)                                                                                                                                          cccGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCACCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGcc             >  4:158192/1‑149 (MQ=255)                                                                                                                                           ccGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAGGCGAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGGTCAGACACGTGGCAGACCCGCAGCACGTGCGCATCAGCACAGCGccc            >  6:26630/1‑149 (MQ=255)                                                                                                                                           ccGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTGATGAGCaagtgctgc            >  3:77214/1‑142 (MQ=255)                                                                                                                                            cGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCg           <  6:17273/149‑1 (MQ=255)                                                                                                                                                     gcgAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCaccacc  >  7:260362/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                      TTATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACC  >  NC_002947/718313‑718609

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