Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	915,857	C→T	9.0%	A5V (GCA→GTA)	sufB →	ABC superfamily ATP binding cassette transporter, membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	915,857	0	C	T	9.0%	132.8 / 2.7	44	A5V (GCA→GTA)	sufB	ABC superfamily ATP binding cassette transporter, membrane protein
Reads supporting (aligned to +/- strand):  ref base C (17/23);  new base T (2/2);  total (19/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 = 9.95e-01

GATATTAGACATATAAAAGTATAAAAATTTTTATAAGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAGAGAAATTTCTAACATGAAAAATGAGCCGGAATGGATGTTAGATTTCCGT  >  CP000730/915719‑915996                                                                                                                                           |                                                                                                                                            gATATTAGACATATAAAAGTATAAAAATTTTTATAAGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCAc                                                                                                                                           <  2:258849/141‑1 (MQ=255)   tatTAGACATATAAAAGTATAAAAATTTTTATAAGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCt                                                                                                                                         <  2:64515/141‑1 (MQ=255)             ataAAAGTATAAAAATTTTTATAAGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTgggg                                                                                                                               <  1:203833/141‑1 (MQ=255)                      taAAAATTTTTATAAGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAat                                                                                                                      <  1:267580/141‑1 (MQ=255)                           aTTTTTATAAGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGa                                                                                                                 <  1:238822/141‑1 (MQ=255)                                tataAGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCa                                                                                                            >  2:398100/1‑141 (MQ=255)                                tataAGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCa                                                                                                            >  2:180584/1‑141 (MQ=255)                                  taAGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCacg                                                                                                          >  2:446917/1‑141 (MQ=255)                                   aaGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGTACCTGATGTTGGGGATTATAAATATGGATTCCacga                                                                                                         >  1:201020/1‑141 (MQ=255)                                   aaGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGTACCTGATGTTGGGGATTATAAATATGGATTCCacga                                                                                                         >  1:293493/1‑141 (MQ=255)                                   aaGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCacga                                                                                                         <  2:315014/141‑1 (MQ=255)                                   aaGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCacga                                                                                                         <  2:408871/141‑1 (MQ=255)                                             gtcatTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGAtgtg                                                                                               <  1:7580/141‑1 (MQ=255)                                                  tGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCAt                                                                                          <  1:107454/141‑1 (MQ=255)                                                  tGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCAt                                                                                          <  1:185976/141‑1 (MQ=255)                                                   gTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCAtt                                                                                         <  1:6942/141‑1 (MQ=255)                                                      ataataTGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTtc                                                                                      >  1:286996/1‑141 (MQ=255)                                                         atatGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTtcaga                                                                                   >  1:219818/1‑141 (MQ=255)                                                           atGGTTTACATCATGAATTAAAAACTTACGCACGTCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTtcagatc                                                                                 <  2:136581/141‑1 (MQ=255)                                                           atGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTtcagatc                                                                                 <  1:309131/141‑1 (MQ=255)                                                           atGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTtcagatc                                                                                 <  2:6513/141‑1 (MQ=255)                                                             ggTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTtcagatcag                                                                               >  1:301083/1‑141 (MQ=255)                                                                 tACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACg                                                                           >  2:425658/1‑141 (MQ=255)                                                                 tACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACg                                                                           >  2:266811/1‑141 (MQ=255)                                                                  aCATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGt                                                                          >  1:32604/1‑141 (MQ=255)                                                                     tcatGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGt                                                                       >  1:141458/1‑141 (MQ=255)                                                                           aTTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACt                                                                 <  2:227123/141‑1 (MQ=255)                                                                               aaaaCTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTgaga                                                             >  2:279169/1‑141 (MQ=255)                                                                               aaaaCTTACGCACGCCGTGGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTgaga                                                             <  2:113991/141‑1 (MQ=255)                                                                                      acgcacgcCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGt                                                      <  1:446917/141‑1 (MQ=255)                                                                                         cacgcCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTag                                                   <  2:6188/141‑1 (MQ=255)                                                                                         cacgcCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTag                                                   <  2:398844/141‑1 (MQ=255)                                                                                            gcCGTTGTAAATATATTTTTAAGGGGTGATTGAAATGGCTAAAAAAGTACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAgaga                                                <  2:201020/141‑1 (MQ=255)                                                                                            gcCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGTACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAgaga                                                <  2:293493/141‑1 (MQ=255)                                                                                                 tGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAGAGAAAttt                                           >  2:248148/1‑141 (MQ=255)                                                                                                           tttttAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTCGGGGACGACAAAGCTGGATTCCACGACGATGATGTGTCCATTTTCAGACCCGAACGGGGGTCAACTGAGAATATCGTTAGAGAAACTTCTAAcgcgaa                                 >  2:69249/1‑136 (MQ=255)                                                                                                                aaGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAGAGAAATTTCTAACATGAAAAATg                            <  1:40312/141‑1 (MQ=255)                                                                                                                        aTTGAAATGACTAAAAAAGCAACTGATGTTGGGGATTATAAATATGGATTCCACGACAATGATGTGTCCATTTTCAGATCAGAACGTGGTTTATCTGTGAAGATCGTTAGAGAAATAAAAAAAAAAAAAAATGAGCCGGaa                    <  2:376426/141‑1 (MQ=255)                                                                                                                            aaaTGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAGAGAAATTTCTAACATGAAAAATGAGCCGGAATGGa                >  2:182734/1‑141 (MQ=255)                                                                                                                             aaTGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAGAGAAATTTCTAACATGAAAAATGAGCCGGAATGGAt               >  2:25825/1‑141 (MQ=255)                                                                                                                             aaTGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAGAGAAATTTCTAACATGAAAAATGAGCCGGAATGGAt               >  1:323173/1‑141 (MQ=255)                                                                                                                             aaTGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAGAGAAATTTCTAACATGAAAAATGAGCCGGAATGGAt               >  1:237713/1‑141 (MQ=255)                                                                                                                                      aaaaGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAGAGAAATTTCTAACATGAAAAATGAGCCGGAATGGATGTTAGAttt      <  2:296909/141‑1 (MQ=255)                                                                                                                                          gCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAGAGAAATTTCTAACATGAAAAATGAGCCGGAATGGATGTTAGATTTCCGt  <  2:219818/141‑1 (MQ=255)                                                                                                                                           |                                                                                                                                            GATATTAGACATATAAAAGTATAAAAATTTTTATAAGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAGAGAAATTTCTAACATGAAAAATGAGCCGGAATGGATGTTAGATTTCCGT  >  CP000730/915719‑915996

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