Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	718,457	C→G	9.2%	intergenic (+27/+50)	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,457	0	C	G	9.2%	124.6 / 3.4	52	intergenic (+27/+50)	thiO/PP_0613	FAD‑dependent glycine/D‑amino acid oxidase/amidase family protein
Reads supporting (aligned to +/- strand):  ref base C (21/25);  new base G (2/3);  total (24/28)
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.19e-01

CATTATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCG  >  NC_002947/718311‑718600                                                                                                                                                   |                                                                                                                                                 cATTATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTa                                                                                                                                                <  4:47183/149‑1 (MQ=255) cATTATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCCCTCCTa                                                                                                                                                >  7:241886/1‑149 (MQ=255)  aTTATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTTAACCCGCTCCTAc                                                                                                                                               >  3:23575/1‑149 (MQ=255)   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                                                                                                                           <  6:102853/145‑1 (MQ=38)                         gcgcCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCGCGTAGGAGCGGGTTT‑ACCcgcg                                                                                                                           >  5:102853/1‑145 (MQ=38)                         gcgcCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTT‑ACCCGCGaaa                                                                                                                        <  6:274252/148‑1 (MQ=255)                             cGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCgg                                                                                                                                       >  2:41661/1‑130 (MQ=255)                             cGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCgg                                                                                                                                       <  1:41661/130‑1 (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)                                                                                                                                                   |                                                                                                                                                 CATTATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCG  >  NC_002947/718311‑718600

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