Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	326,176	G→T	P358T (CCG→ACG)	betA ←	choline dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	326,176	0	G	T	97.7%	126.6 / NA	44	P358T (CCG→ACG)	betA	choline dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (21/22);  minor base C (0/1);  total (21/23)
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.96e-01

GTCGCGGGATTTAATCCGCACATGCCCACGGCTTGGCGAGCGCATTGAGCCGACGTGGCACTGGAAACCGTGCTCTTTCACTGCATTCGAGCCGTTATAGTTAATCGCTACTGGCAGGAAATGGTACTGAATATTCGGCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCAAAGTGGTTGCTGGCACCAACGCCAGTGCCGCCAAACAGCCACTCCGCACCGATTTTCGGCTGGTTCCACCACTGCAGGGCAGGGTAGAGGG  >  NC_000913/326039‑326310                                                                                                                                          |                                                                                                                                       gtggcgGGACTTGCTCCGGCCCTGCCCACGGCTTGGCGAGCGCATTGAGCCGACGTGGCACTGGAAACCGTGCTCTTTCACTGCATTCGATCCGTTATAGTTAATCGCTACTGGCAGGAAATGGTATTGCTTATTCGTCCACg                                                                                                                                   <  2:891154/140‑1 (MQ=255) ggcgccGGTTTTAGTCCCCCCCTGCCCCGGGCTTGGCGAGCGGATTGTGCCGACGTGGCAGTGGAAACCGTGCTCTTTCACTGCATTCGAGCCGTTATAGTCAACCGCTCCTGGCAGGAAATGGTACTGAATATTCGTCCACg                                                                                                                                   <  2:891391/141‑1 (MQ=255)     cgGGATGTAATCCGCACATGCCCACGGCTTGGCGAGCGCATTGAGCCGACGTGGCACTGGAAACCGTGCTCTTTCACTGCATTCGAGCCGTTATAGTTAATCGCTACTGGCAGGAAATGGTACTGAATATTCGTCCACGCaaa                                                                                                                               <  1:318620/143‑1 (MQ=255)       ggATTTAATCCGCACATGCCCACGGCTTGGCGAGCGCATTGAGCCGACGTGGCACTGGAAACCGTGCTCTTTCACTGCATTCGAGCCGGTATAGTTAGTCGCTACTGGCACGCAATGGTACTGAACACTCGTCCACGGAAGtt                                                                                                                             >  2:546219/1‑143 (MQ=255)          tttAATCCGCACATGCCCACGGCTTGGCGAGCGCATTGAGCCGACGTGGCACTGGAAACCGTGCTCTTTCACTGCATTCGAGCCGTTATAGTTAATCGCTACTGGCAGGAAATGGTACTCACTATCCGTCCCCGCCACTTCCt                                                                                                                          >  2:34588/1‑143 (MQ=255)                    acaTGCCCACGGCTTGGCGAGCGCATTGAGCCGACGTGGCACTGGAAACCGTGCTCTTTCACTGCATGCGAGCCGTTATAGTTAATCGCTACTGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATCCCTGACGGCTGCg                                                                                                                >  1:157437/1‑143 (MQ=255)                      gtGCCCTCGGCTGGGCGGGCGGATTGAGCCGACGTGGGACTGGAAACCGTGCTCTTTCACTGCATTCGAGCCGTTATAGTTAATCGCTACTGGCAGGAAAGGGTACTTAATATTCGTCCACGCAAATTCCTCACGGCTGCGaa                                                                                                              <  2:496605/142‑1 (MQ=255)                      aTGCCCACGGCTTGGCGAGCGCATTGAGCCGACGTGGCACTGGAAACCGTGCTCTTTCACTGCATTCGAGCCGTTATAGTTAATCGCTACTGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGaa                                                                                                              >  1:797449/1‑143 (MQ=255)                       tGCCCACGGCTTGGCGTGCGCATTGAGCCGTCGTGGCACTGGAAACCGTTCTCTTTCACTGCATTCGAGCCGTTATAGTTAATCGCTACTGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAAt                                                                                                             <  1:76282/143‑1 (MQ=255)                         cccACGGCTTGGCGCGCGCATTGGGCCGACGTGGCACTGGAAACCGTGCTCTTTCACTGCATTCGAGCCGTTATAGTTAATCGCTACTGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAATaa                                                                                                           <  2:289143/143‑1 (MQ=255)                          ccACGGCGCGGCGAGCGCATTGAGCCGACGTGGCACTGGAACCCGTGCTCTTTCACTGCATTCGAGCCGTTATAGTTAATCGCTACTGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAATaaa                                                                                                          <  1:269706/143‑1 (MQ=255)                             gggCTTGGCGAGCGCCTTGGGCCGCCGTGGCACTGGAAACCGTGCTCTTTCTCTGCATTCGAGCCGTTATAGTTAATCGCTACTGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAATAAATcc                                                                                                       <  2:814647/142‑1 (MQ=255)                                         cgcATTGGGCCGAGGTGCCTGTGGAACTCGTGGTCTTTCTCTGCATTCGTGCCGTTATAGTTAATCGCTATTGGCAGGAAATGGTACTGAATATTCTTCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCaca                                                                                           <  2:137246/143‑3 (MQ=255)                                         cgcATTGAGCCGACGTGGCACTGGAAACCGTGCTCTTTCACTGCATTCGAGCCGTTATAGTTAATCGCTACTGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCaaa                                                                                           >  1:884572/1‑143 (MQ=255)                                                  ccGACGTGGCACTGGAAACCGTGCTCTTTCACTGCATTCGAGCCGTTATAGTTAATCGCTACTGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCAAAGTGGTTGCt                                                                                  >  1:196778/1‑143 (MQ=255)                                                  ccGACGTCGCCCTGGAAACCGTGCGCTTTCACTGCATTCGAGCCGTTATAGTTAATCGCTACTGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCAAAGTGGTTGCt                                                                                  <  1:654559/143‑1 (MQ=255)                                                   cGACGTGGCACTGGAAACCATGCTCTTTCACTGCATTCGAGCCGTTATAGTTAATCGCTACTGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCAAAGTGGTTGCTg                                                                                 >  1:492956/1‑143 (MQ=255)                                                      ggTGGCACTGGAAACCTGGCTCTTTCACTGCATCCGAGCCGCTATAGTTAATCGCTACTGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCAAAGTGGTTGCTGGCa                                                                              <  2:131778/142‑1 (MQ=255)                                                      cGTGGCACTGGATACCGTGCTCTTTCACTGCATTCGAGCCGTTATAGTGAATCGCTACTGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCAAAGTGGTTGCTGGCa                                                                              <  2:764545/143‑1 (MQ=255)                                                           gagtGGAAACGGTGCTCTTTCACTGCATTCGAGCCGTTACAGTCAATCGTTTTTGGCAGGAAAGGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCAAAGTGGTTGCTGGCACCAAc                                                                         <  1:565997/140‑1 (MQ=255)                                                            aCTGGAAACCGTGCTCTTTCACTGCATTCGAGCCGTTATAGTTAATCGCTACTGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCAAAGTGGTTGCTGGCACCAACg                                                                        >  1:379179/1‑143 (MQ=255)                                                             cTGGAATCCCTGCTCGTTCACTGCATTCGAGCCGTTATAGTTAATCGCTACTGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCGCACGGCAGCGAATAAATCCACCTGCTTCAAAGTGGTTGCTGGCACCAACGc                                                                       <  2:726773/143‑1 (MQ=255)                                                                    gcgggCTCCTTCCCCGCATTCGAGCCGTTATAGTGAATCGCTACTGGCAGGAAATGGTACTGAATATTCGCCCACGCCAATTCCCCACGGCTGCGAATAAATCCACATGCTTCAAAGTGGTTGCTGGCACCAACGCCAGTgcc                                                                <  1:90673/139‑1 (MQ=255)                                                                      gtcctctTTCACTGCATTCGTGCCGTTATAGTTAATCGCGACTGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCAAAGTGGTTGCTGGCACCAACGCCAGTgccgc                                                              <  2:229359/140‑1 (MQ=255)                                                                        gCTCTTTCACTGCATTCGAGCCGTTATAGTTAATCGCTACTGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCAGGGCTGCGACTAAAACCACCTGCTTCAAAGTGGTTGCTGGCACCAACGCCAGTGCCGCCa                                                            >  1:380140/1‑143 (MQ=255)                                                                          tctTTCACTGCATTCGAGCCGTTATAGTTAATCGCTACTGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCAAAGTGGTTGCTGGCACCAACGCCCGTGCCGCCaac                                                          >  1:28655/1‑142 (MQ=255)                                                                                 cTGCATTCGAGCCGTTATAGTTAATCGCTACTGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCAAAGTGGTTGCTGGCACCAACGCCAGTGCCGCCAAACACCCAc                                                   >  1:695834/1‑143 (MQ=255)                                                                                   gCATTCGAGCCGTTATATTTAATCGCTATTGGCCGGACCTGGTAATGAATATTCGTTATCGCTCTTTCCTCACGGCTGCGAATAAATCCACCTGCTTCAATGTGTTTACTCACACCAACGCCAGTGCCGCCAAACAGCCACTc                                                 <  2:16874/143‑1 (MQ=255)                                                                                      ttCGAGCCGTTATAGTTAATCGCTACTGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCAAAGTGGTTGCTGGCACCAACGCCAGTGCCGCCAAACAGCCACTCCGc                                              <  1:631253/143‑1 (MQ=255)                                                                                                  taGTTAATCGCTACTGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCAAAGTGGTTGCTGGCACCAACGCCAGTGCCGCCAAACAGCCACTCCGCACCGATTTTCgg                                  <  1:843710/143‑1 (MQ=255)                                                                                                     ttAATCGCTACTGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGAGAATAAATCCACATGCTTCAAAGTGGTTGCTGGCACCGACGCCAGTGCCGCCAGACAGCCACTCCGCACCGAGTTTCggggg                               >  1:146850/1‑140 (MQ=255)                                                                                                        aTCGCTACTGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCAAAGTGGTTGCTGGCACCAACGCCAGTGCCGCCAAACAGCCACTCCGCACCGATTTTCGGCtcgct                            >  2:414737/1‑139 (MQ=255)                                                                                                            cTACTTGCAGGAATTGGTACTGAATTTCTGTCCGCGCAGATTCCTCTCGGCTGCCAATAACTCCTCCGGCGTCAAAGCGGTTGCTGGCCCCAACGCCAGTGCCGCCAAACAGCCACTCCGCACCGATTTTCGGCTGGTTccac                        <  2:434942/143‑1 (MQ=255)                                                                                                                tGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCAAAGTGGTTGCTGGCACCAACGCCAGTGCCGCCAAACAGCCACTCCGCACCGATTTTCGGCTGGTTCCACCACt                    <  2:564564/143‑1 (MQ=255)                                                                                                                tGGCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCAAAGTGGTTGCTGGCACCAACGCCAGTGCCGCCAAACAGCCACTCCGCACCGATTTTCGGCTGGTTCCACCACt                    >  1:793956/1‑143 (MQ=255)                                                                                                                  gCAGGAAATGGTACTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCAAAGTGGTTGCTGGCACCAACGCCAGTGCCGCCAAACAGCCACTCCGCACCGATTTTCGGCTGGTTCCACCACTgc                  >  1:798718/1‑143 (MQ=255)                                                                                                                  gCAGGAAATGGTACTGAATATTCCTCCACGCAAATTCCTCACGACAGCGACAACTTCGACCTGCTTAAACGTGTAGACTTGCACCAACGCCAGTGCCGCCAAACAGCCACTCCGCACCGATTTTCGGCTGGTTCCACCACTgc                  >  1:905734/1‑143 (MQ=255)                                                                                                                    aGGAAATGGTACTGAATATTCGTCCACGCCAATGCCTAACGGCTGCGGATAAATCCACCTGCTACAACGTGGTTGGTGGCACCAAAGCCAGTGCGGCCAAACAACGACTCCGCACAGATTTTCGGCGGGTCCTACCCGTgcag                >  2:176394/1‑143 (MQ=255)                                                                                                                               cTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCAAAGTGGTTGCTGGCACCACCCAAAACGCTGGCTTACAGCCCCTCCGCACCGATTTTCGGCTGGTTCCACCACTGCAGGGCAGGGTagc     >  1:565691/1‑142 (MQ=255)                                                                                                                               cTGAATATTCGTCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCAAAGTGGTTGCTGGCACCAACGCCAGTGCCGCCAAACAGCCACTCCGCACCGATTTTCGGCTGGTTCCACCACTGCGGGGCAGGGTaga     >  1:223002/1‑143 (MQ=255)                                                                                                                                tGAATATTCGTCCACGGAAATTCCGCACGGCTGCGAATAAATCCACCTGCTTCAAAGTGGTTGCTGGCACCAACGCCGGTGCCGCCAAACAGCCACTCCGCACCGATTTTCGGCTGGTTCCACCACTGCAGGGCAGGGTagag    <  1:34588/143‑1 (MQ=255)                                                                                                                                ggactatTCGTCCACGCAAATTCCGCACGGCGGCGAATAAATCCACCTGCTTCAAAGTGGTTGCTGGCACCAACGCCAGTGCCGCCAAACAGCCACTCCGCACCGATTTTCGGCTGGTTCCACCACTGCAGGGCAGGGTagag    <  2:856073/139‑1 (MQ=255)                                                                                                                                  aaTATTCGTCCACGCAAATTCCTCCCGGCTGCGAATAAATCGACCTGCTTCAAAGTGGTTGCTGGCACAAACGACAGTGCCGCCAAACAGCCACTCCGCACCGATTTTCGGCTGGTTCCACCACTGCGGGGCAGGGTGGAggg  >  2:788561/1‑143 (MQ=255)                                                                                                                                  aaTATTCGTCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCAAAGTGGTTGCTGGCACCAACGCCAGTGCCGCCAAACAGCCACTCCGCACCGATTTTCGGCTGGTTCCACCACTGCAGGGCAGGGTAGAggg  >  1:670720/1‑143 (MQ=255)                                                                                                                                          |                                                                                                                                       GTCGCGGGATTTAATCCGCACATGCCCACGGCTTGGCGAGCGCATTGAGCCGACGTGGCACTGGAAACCGTGCTCTTTCACTGCATTCGAGCCGTTATAGTTAATCGCTACTGGCAGGAAATGGTACTGAATATTCGGCCACGCAAATTCCTCACGGCTGCGAATAAATCCACCTGCTTCAAAGTGGTTGCTGGCACCAACGCCAGTGCCGCCAAACAGCCACTCCGCACCGATTTTCGGCTGGTTCCACCACTGCAGGGCAGGGTAGAGGG  >  NC_000913/326039‑326310

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