Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	718,460	C→T	10.7%	intergenic (+30/+47)	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,460	0	C	T	10.7%	115.7 / 3.9	48	intergenic (+30/+47)	thiO/PP_0613	FAD‑dependent glycine/D‑amino acid oxidase/amidase family protein
Reads supporting (aligned to +/- strand):  ref base C (18/23);  new base T (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 = 6.13e-01

ATTATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCG  >  NC_002947/718312‑718600                                                                                                                                                     |                                                                                                                                              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                                                                                                                           >  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)                                                                                                                                                     |                                                                                                                                              ATTATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGCGAGCGGGTTT‑ACCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCG  >  NC_002947/718312‑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: ‑