Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	356,058	A→G	L379L (TTA→TTG)	codB →	cytosine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	356,058	0	A	G	98.0%	172.2 / NA	50	L379L (TTA→TTG)	codB	cytosine transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (18/31);  minor base T (1/0);  total (19/31)
Fisher's exact test for biased strand distribution p-value = 3.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACA  >  NC_000913/355921‑356190                                                                                                                                          |                                                                                                                                     ccAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGcc                                                                                                                                    <  2:98118/140‑1 (MQ=255)   aGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCgg                                                                                                                                  >  1:417608/1‑140 (MQ=255)   aGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGCCGGACGCATTGGGGATGGCTGTGGGCCACTGGTTGCCgg                                                                                                                                  >  2:325763/1‑140 (MQ=255)       ggtggCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAAt                                                                                                                              >  2:361783/1‑140 (MQ=255)       ggtggCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAAt                                                                                                                              >  1:113640/1‑140 (MQ=255)            cGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTc                                                                                                                         >  1:596215/1‑140 (MQ=255)              tGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCg                                                                                                                       >  2:41867/1‑140 (MQ=255)                   atcGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCaa                                                                                                                  >  2:142369/1‑140 (MQ=255)                   atcGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGAGTGCAGCACGCCACTGGGTGCCGGGAATTGTTCCGGTCaa                                                                                                                  >  2:263153/1‑140 (MQ=255)                    tcGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAAc                                                                                                                 <  2:411347/140‑1 (MQ=255)                     cGCCGACTATCTGATGAACCGTTGCAGCTCTGAGCACATTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAAcg                                                                                                                <  1:263153/140‑1 (MQ=255)                     cGCCGACTATCTGATGAACCGTGGCGGCTATGAGCACCTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAAcg                                                                                                                <  1:239420/140‑1 (MQ=255)                     cGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAAcg                                                                                                                <  1:114314/140‑1 (MQ=255)                     cGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAAcg                                                                                                                <  1:73896/140‑1 (MQ=255)                     cGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAAcg                                                                                                                <  1:713223/140‑1 (MQ=255)                     cGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAAcg                                                                                                                <  1:448672/140‑1 (MQ=255)                     cGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGATCAAcg                                                                                                                <  2:596215/140‑1 (MQ=255)                     cGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCATGGGGATTGCTGTAGGACACTGGTTTCCGGAAATTGTTCCGGTCAAcg                                                                                                                >  2:235291/1‑140 (MQ=255)                     cGCCGACTATCTGATGAACCGTCGCCGCTATCAGCACTTGGCGACCACACGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAAcg                                                                                                                <  2:277009/140‑1 (MQ=255)                        cGACTATCTTAAGAACCGGCGCCGCTATTAGCACCTTGCGGCCACGCGTATGATGAGTGGCAATTGGGTGGCGATTCTGGCGGGCGGCTTGGGGATTGCTGCAGGCCGCTGGTTGCCGGGAATTGTTCCGGTCAAccggg                                                                                                             <  2:571771/140‑5 (MQ=255)                          aCTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCaa                                                                                                                  >  1:425338/1‑133 (MQ=255)                          aCTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCaa                                                                                                                  <  2:425338/133‑1 (MQ=255)                                  atgaACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAAcg                                                                                                                >  1:496806/1‑127 (MQ=255)                                  atgaACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAAcg                                                                                                                <  2:496806/127‑1 (MQ=255)                                     aaCCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGgcgc                                                                                                <  1:41867/140‑1 (MQ=255)                                        cGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCTCTGGTTGCAGGGGATTTTTCCGGTCAACGCCGTATTAGGTGGCGCGCt                                                                                             >  2:252734/1‑140 (MQ=255)                                           cgccgcTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAg                                                                                          <  1:53696/140‑1 (MQ=255)                                                  aTGAGCACTTTGCGACCACGCGCATGATGCGTGTCAAGTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTg                                                                                   <  1:252734/140‑1 (MQ=255)                                                       cACTTTGCGACCACGCGTATGACGACTGTCAATTGGGAGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCt                                                                              <  1:189415/140‑1 (MQ=255)                                                           ttGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAAc                                                                          >  2:119661/1‑140 (MQ=255)                                                            tGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAAcc                                                                         <  2:597999/140‑1 (MQ=255)                                                            tGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAAcc                                                                         <  1:35150/140‑1 (MQ=255)                                                             gCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAACGCGGTATTAGGGGGCACGATGAGCTGTCTGATCCTTAAccc                                                                        >  1:218216/1‑140 (MQ=255)                                                                               aGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCGGAGTTATATGATCCTTGACCCGATTTTGAATCGTAAAAc                                                      >  2:256939/1‑140 (MQ=255)                                                                                  gtCAATTGGGTGGCGATTCTGGGGGGCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGAc                                                   <  1:330264/140‑1 (MQ=255)                                                                                   tCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGAca                                                  >  2:613040/1‑140 (MQ=255)                                                                                   tCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTCCGAATCGTAAAACGAca                                                  >  2:330264/1‑140 (MQ=255)                                                                                    cAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGAcag                                                 <  2:413886/140‑1 (MQ=255)                                                                                     aactGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGAcagc                                                <  2:491162/137‑1 (MQ=255)                                                                                      aTTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCAGATCCTTCACCCGAATTTGAAGCGTAAAACGAcagca                                               >  1:300738/1‑140 (MQ=255)                                                                                         gggTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAgcagca                                            <  1:534278/140‑1 (MQ=255)                                                                                                           tCGCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCt                          <  2:417608/140‑1 (MQ=255)                                                                                                             gCCTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTaa                        <  1:6792/140‑1 (MQ=255)                                                                                                               cTTGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACa                      <  1:142369/140‑1 (MQ=255)                                                                                                                 tGGGGATTGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAgt                    >  2:484628/1‑140 (MQ=255)                                                                                                                       ttGCTGCAGGCCACTGGTTGCCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGaa              <  2:440596/140‑1 (MQ=255)                                                                                                                       ttGCTGCAGGCCACTGGTTGCCCGGAATTGTTCCGCTCATCGCGGTATGAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGaa              <  2:225502/140‑1 (MQ=255)                                                                                                                        tgctgcAGGCCACTGGTTGCCGCGAGTTGATCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAAt             <  1:256939/140‑1 (MQ=255)                                                                                                                                   aCTGGTTGCCGGGAATTTTTCCCGTCAACGCAGTATTAGTTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTaca  <  2:218216/140‑1 (MQ=255)                                                                                                                                   aCTGGTTGCCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTaca  <  1:416872/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                     CCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACA  >  NC_000913/355921‑356190

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑