Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	616,451	T→A	7.1%	noncoding (71/142 nt)	PP_mr11 ←	FMN regulatory region

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	616,451	0	T	A	7.1%	157.1 / 3.0	56	noncoding (71/142 nt)	PP_mr11	FMN regulatory region
Reads supporting (aligned to +/- strand):  ref base T (29/23);  new base A (2/2);  total (31/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 = 3.33e-01

GCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCGGCGAGACTTGCAGCAGT  >  NC_002947/616307‑616594                                                                                                                                                 |                                                                                                                                                gCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTc                                                                                                                                             <  3:312870/149‑1 (MQ=255)     ctcGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACCATGTCgggg                                                                                                                                         <  3:258017/149‑1 (MQ=255)     ctcGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCgg                                                                                                                                           <  1:10230/147‑1 (MQ=255)      tcGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGgc                                                                                                                                        >  2:58706/1‑149 (MQ=255)             tGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTcgcg                                                                                                                                 >  1:225872/1‑149 (MQ=255)               aTCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCgg                                                                                                                                <  7:282251/148‑1 (MQ=255)                      aaCGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACCATGTCgggggtcag                                                                                                                                    >  6:315511/1‑132 (MQ=255)                      aaCGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACCATGTCgggggtcag                                                                                                                                    <  5:315511/137‑6 (MQ=255)                               gggCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGgcgc                                                                                                                                      <  8:64879/126‑1 (MQ=255)                                   gAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGCCCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTAccgc                                                                                                           >  3:169650/1‑149 (MQ=255)                                     aCAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTAccgccg                                                                                                         <  6:25153/149‑1 (MQ=255)                                             cgcgcgcAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGa                                                                                                  <  7:243287/148‑1 (MQ=255)                                                gcgcAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGCGCTCATCTATCGATTTACCGCCGGTGGGGACTTa                                                                                              >  8:247032/1‑149 (MQ=255)                                                    aaGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACAcc                                                                                          >  7:12143/1‑149 (MQ=255)                                                      gggCGCCCGGATGGCCCCGGCATGACGCCTTCTTTCATCCGGACTATGCCCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACAcccc                                                                                        <  2:322710/149‑1 (MQ=255)                                                       ggCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCg                                                                                       >  6:230553/1‑149 (MQ=255)                                                         cgcCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGcc                                                                                     <  2:371380/149‑1 (MQ=255)                                                              ggATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGaa                                                                                <  2:225872/149‑1 (MQ=255)                                                              ggATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGaa                                                                                <  1:58706/149‑1 (MQ=255)                                                                   gCCCCTTCATGTCGCCTTATTTCATCCGGACTATTAGCGTCGGCTCTGGAGTTTCACCAGATCTTCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAgacc                                                                            <  2:369484/148‑1 (MQ=255)                                                                     cccTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAgaccgac                                                                         >  4:374782/1‑149 (MQ=255)                                                                     cccTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAgac                                                                             <  2:351469/145‑1 (MQ=255)                                                                     cccTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAgac                                                                             >  1:351469/1‑145 (MQ=255)                                                                      ccTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAgaccgacc                                                                        >  8:236105/1‑149 (MQ=255)                                                                           aTGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACAGACCCAGc                                                                    >  7:360024/1‑148 (MQ=255)                                                                                gCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAg                                                              >  4:294722/1‑149 (MQ=255)                                                                                 ccTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAg                                                              >  3:15754/1‑148 (MQ=255)                                                                                     cttTCATCCGGACTATGACCGTCGTCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGcaca                                                          <  1:66365/148‑1 (MQ=255)                                                                                      tttCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGt                                                        >  4:378856/1‑149 (MQ=255)                                                                                           tCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCAc                                                   >  7:193765/1‑149 (MQ=255)                                                                                                  tATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTg                                            >  8:350047/1‑149 (MQ=255)                                                                                                     gACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTAACGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCAAGCATACAGCACAGTGCCACCCCAGTGGc                                          >  6:217301/1‑148 (MQ=255)                                                                                                       ccGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTC‑GGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAAcc                                      >  4:29861/1‑149 (MQ=255)                                                                                                              ctctGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGt                                <  7:236105/149‑1 (MQ=255)                                                                                                                   gagaTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTctac                           <  3:50567/149‑1 (MQ=255)                                                                                                                   gAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCTTTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCA‑CCCAGTGGCAACCCGCGGTTctac                           <  4:169650/148‑1 (MQ=255)                                                                                                                    aGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTctact                          <  3:378856/149‑1 (MQ=255)                                                                                                                         cACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCtt                     >  2:196824/1‑149 (MQ=255)                                                                                                                           ccAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCg                   >  6:264771/1‑149 (MQ=255)                                                                                                                           ccAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCg                   >  1:67584/1‑149 (MQ=255)                                                                                                                            cAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCgg                  >  1:384006/1‑149 (MQ=255)                                                                                                                             aGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCGGc                 <  8:360024/149‑1 (MQ=255)                                                                                                                             aGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCGGc                 <  3:294722/149‑1 (MQ=255)                                                                                                                                   gctgACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCGGCGAGACt           >  4:133324/1‑149 (MQ=255)                                                                                                                                    ctgACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGaa                                                                                >  4:245500/1‑79 (MQ=255)                                                                                                                                    ctgACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGaa                                                                                <  3:245500/79‑1 (MQ=255)                                                                                                                                    ctgACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCGGCgaga             <  7:247032/146‑1 (MQ=255)                                                                                                                                       aCCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCGGCGAGACTTgca       <  5:230553/149‑1 (MQ=255)                                                                                                                                        gccccGACCATGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCGGCGAGACTTgcag      >  1:2387/2‑149 (MQ=255)                                                                                                                                         ccccGACATGGTCGGGGCGCTCGCGGGCTCATGTTTCGATTTACCGCCGGTGGGGACGTACACCCCGCCCTGCAGACCGACCCAGCATACAGCACAGTGCCACCGCAGTGGCAACCCGCGGTTCTACTACCTTCGGCGAGACTTgcag      >  3:295999/1‑148 (MQ=255)                                                                                                                                         ccccGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCGGCGAGACTTgcag      >  8:15004/1‑148 (MQ=255)                                                                                                                                         ccccGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCGGCGAGACTTgcag      >  3:366892/1‑148 (MQ=255)                                                                                                                                         ccccGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCGGCGAGACTTgcag      >  2:19905/1‑148 (MQ=255)                                                                                                                                           ccGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCGGCGAGACTTGcagcag   >  8:147217/1‑149 (MQ=255)                                                                                                                                           ccGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCGGCGAGACTTGcagcag   >  1:260750/1‑149 (MQ=255)                                                                                                                                            cGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCGGCGAGACTTGCAGCAGt  <  4:15754/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                GCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCGGCGAGACTTGCAGCAGT  >  NC_002947/616307‑616594

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