Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	707,910	A→G	intergenic (‑1/‑121)	mntB ← / → sirR	ABC superfamily ATP binding cassette transporter, ABC protein/possible metal‑dependent transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	707,910	0	A	G	96.5%	179.4 / ‑4.7	58	intergenic (‑1/‑121)	mntB/sirR	ABC superfamily ATP binding cassette transporter, ABC protein/possible metal‑dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (30/26);  total (32/26)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

AAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAACTGAAGAAAAAGAGG  >  CP000730/707768‑708052                                                                                                                                               |                                                                                                                                               aaGACTTGATAGGGGGGGATTTACCAGCACCATTCGGGCCCATGATACCAATCATTTAGCAGCGTCATGGTGTTGATTTGGAACAGAGGTTAAGTACATCCTGATTAATTTAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGttt                                                                                                                                           >  2:245819/1‑148 (MQ=255) aaGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGttt                                                                                                                                           >  1:439419/1‑148 (MQ=255) aaGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGttt                                                                                                                                           >  1:75516/1‑148 (MQ=255) aaGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGttt                                                                                                                                           >  1:136479/1‑148 (MQ=255)   gACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACAGGTATCGATAAGGCAATGTTTTTATGTACATGCTTATTATCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGGTtat                                                                                                                                         >  2:92243/1‑148 (MQ=255)    aCTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTtata                                                                                                                                        >  1:294937/1‑148 (MQ=255)      ttGAGAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATAcc                                                                                                                                      >  2:560769/1‑148 (MQ=255)        gATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATAcctc                                                                                                                                    <  2:564844/148‑1 (MQ=255)         gtgaGGCAAGATTTCCCAGCACCTTTGGGGCCCAGGATACCAATTATTTCGCCGCGGACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCCTATTACGTAAAGACAGATTTAAATCTGTGGTTTCTAACAGGCGTTTATAcctcc                                                                                                                                   <  2:576031/145‑1 (MQ=255)         aTAAGGGAACGTGTACCCGCACCATTCGGGCCCAGGNTACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTANCTAAAAACAGATTTAAATCGTTTGTTTCTAACAAGCGTTTATAcctcc                                                                                                                                   <  2:373602/148‑1 (MQ=255)           aaGGCAAGATTTACCAGCACCATTCGGGCCCATGATTCCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTa                                                                                                                                 <  1:153814/148‑1 (MQ=255)             gggAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTACATCTTTTGTTTCTAACAAGAGTTTATACCTCCTAAt                                                                                                                               >  1:35385/1‑148 (MQ=255)             gggAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAAt                                                                                                                               >  1:269020/1‑148 (MQ=255)                aaGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTCATTaa                                                                                                                            >  1:257957/1‑148 (MQ=255)                    tttCCCAGCACCAGTCGGGCCCATGATACCAATTATTTCTCCTCGTACTCGTATCGATAATGAATTTTTTTGAAGTACATGCTTATTACCTAAAAACAGATTTAATTCTTTTGTTTCTAACAAGCGTTTATACCTCCTAAGTAAAAGt                                                                                                                        <  2:23263/148‑1 (MQ=255)                           gCACCATTCGGGCCCATGATACCAATTATTTAGCCGCGTACTGGTATGGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCGAACAAACGTTTATACCTCCTAGTTAATAGTTTAGGCt                                                                                                                 >  2:598232/1‑148 (MQ=255)                             aCCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAACGTTTAGGCTaa                                                                                                               >  2:243296/1‑148 (MQ=255)                                   cGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAACTGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACATGCGTTTATACCTCCTGATTACAAGTTTAGGCTAAcatgct                                                                                                         >  2:155508/1‑143 (MQ=255)                                        ccATGATACCAATTATTTCGGCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTaattaatt                                                                                                    <  2:69323/148‑1 (MQ=255)                                        ccATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTaattaatt                                                                                                    <  2:60733/148‑1 (MQ=255)                                        ccATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTaattaatt                                                                                                    <  2:60731/148‑1 (MQ=255)                                              tACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTataa                                                                                              <  1:455918/148‑1 (MQ=255)                                               aCCAATTATTTCACCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAACACAGATTTAAATCTTTTGTTTTTAACAAGCGATTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGtgtcat                                                                                             >  2:381877/1‑143 (MQ=255)                                                      attTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTga                                                                                      <  2:525247/148‑1 (MQ=255)                                                          cgccgcGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAat                                                                                  >  1:632033/1‑148 (MQ=255)                                                                 tACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATc                                                                           <  1:91918/148‑1 (MQ=255)                                                                         tCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAg                                                                   <  1:624230/148‑1 (MQ=255)                                                                         tCGATAAGGAAATGTTTTTAAGTACATGCTTATTACATAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAg                                                                   <  2:33741/148‑1 (MQ=255)                                                                           gATCAGGAAATGTTTTTCAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTAATACCTCCTACTTAAAAGTTTAGGCTAACCTAATTAATTGTATAATCAACTGGAAATCTTTATCTTGTCAAGTa                                                                 >  2:179145/1‑148 (MQ=255)                                                                                ggAAATGTTTTTAAGTACATGCTTATTAACTAAAAACAGATTTAAATCTTTTTTTTCTAATAAGCGGTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAAGAAACTGAGAATATTTATCGTGTCAAGTAAATTc                                                            >  1:383506/1‑148 (MQ=255)                                                                                 taaaTGTTTTTAAGTACATGTTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCg                                                           <  1:537059/147‑1 (MQ=255)                                                                                 gAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCg                                                           <  1:466186/148‑1 (MQ=255)                                                                                   aaTTTTTTTAAGTACATTCTTATTATCTAAAAACAGATTTAAATCTTTTGGTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTg                                                         <  2:147291/148‑1 (MQ=255)                                                                                   aaTGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTg                                                         <  1:45295/148‑1 (MQ=255)                                                                                   aaTGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTg                                                         >  1:160896/1‑148 (MQ=255)                                                                                    aTGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTACTTGGATAATAATCTGAGAATATTTATCATGTCAAGTAAATTCGTGa                                                        <  1:560769/148‑1 (MQ=255)                                                                                            aaGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGatataatat                                                <  1:324850/148‑1 (MQ=255)                                                                                               tACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGa                                             >  1:545223/1‑148 (MQ=255)                                                                                               tACATGCTGATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGa                                             <  2:452915/148‑1 (MQ=255)                                                                                                      ttattaCCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAatgta                                      <  1:245819/148‑1 (MQ=255)                                                                                                          taCCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTAtgtg                                  <  1:414726/148‑1 (MQ=255)                                                                                                              tAAAAACAGATTTATATCTTTTGTTTCTAACATGCGTTTATACCTCCTAATTAAAAGTTTAGACTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGGTATAATATAGTCATTGTATGTGAGGt                              >  2:450876/1‑148 (MQ=255)                                                                                                                aaaaCAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGa                            <  2:489157/148‑1 (MQ=255)                                                                                                                  aaCAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGaaa                          <  2:437587/148‑1 (MQ=255)                                                                                                                  aaCAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAAGGTATGTGAGGTGaag                          >  1:23040/1‑147 (MQ=255)                                                                                                                  aaCAGATTTAAATCTTTTGTTTCCAACAAGCGTTAATACCTCCTAATAAAAAGTTTAGTCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGCTATAATATAGACCATGTATGTGAGGTGaga                          <  2:294666/148‑3 (MQ=255)                                                                                                                         ttAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGtt                   <  2:517326/148‑1 (MQ=255)                                                                                                                            aaTCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACACTGTATGTGAGGTGAAAGTATGTTAAc                >  2:128729/1‑148 (MQ=255)                                                                                                                            aaTCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAAc                >  2:231275/1‑148 (MQ=255)                                                                                                                             aTCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATACTAAACTGAGAATATTTATCATTTCAAGTAATTTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTACCt               >  1:39084/1‑148 (MQ=255)                                                                                                                             aTCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAACt               >  1:505367/1‑148 (MQ=255)                                                                                                                              tCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAAGTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAACTg              >  1:491086/1‑148 (MQ=255)                                                                                                                                 tttGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATCTTTATCATGTCATGTAAATTCGTGATCTAGTATAAACAATGTATGTGAGGTGAAAGTATGTTAACTgaag           >  2:191747/1‑148 (MQ=255)                                                                                                                                 tttGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAACTgaag           >  2:497489/1‑148 (MQ=255)                                                                                                                                  ttGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAACTgaaga          >  1:185386/1‑148 (MQ=255)                                                                                                                                   tttttCTAACAAGCGTTTATACCTCCTAATTAACGGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAACTgaagaa         >  1:183419/3‑148 (MQ=255)                                                                                                                                       tCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAACGTATGTTAATTGAAGAAAAag     >  2:98767/1‑148 (MQ=255)                                                                                                                                          aaCAAGCGTTTATACCTCCTAATTAAAAGTTTATGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAACTGAAGAAAAAGAgg  >  1:446222/1‑148 (MQ=255)                                                                                                                                               |                                                                                                                                               AAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAACTGAAGAAAAAGAGG  >  CP000730/707768‑708052

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑