Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,833,231	A→C	M436R (ATG→AGG)	yicI ←	putative alpha‑glucosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,833,231	0	A	C	98.0%	174.1 / NA	49	M436R (ATG→AGG)	yicI	putative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (29/19);  minor base G (0/1);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCATTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTAAA  >  NC_000913/3833088‑3833373                                                                                                                                                |                                                                                                                                               aaCAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTg                                                                                                                                         <  1:38357/151‑1 (MQ=255) aaCAATAGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTg                                                                                                                                         >  1:8847/1‑151 (MQ=255)    aaTCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCgg                                                                                                                                      >  2:560961/1‑151 (MQ=255)       cGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGAt                                                                                                                                    <  2:534030/150‑1 (MQ=255)       cGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATc                                                                                                                                   >  1:416940/1‑151 (MQ=255)       cGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATc                                                                                                                                   <  2:260605/151‑1 (MQ=255)         ccgccccAGTGGACCAGGAATTTCTGCGCAGCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCAGACCAGTTCGTTGTCGATGTACGCATAATGATTATGCGTTTTCTGCGGagggg                                                                                                                                 <  1:269056/148‑5 (MQ=255)                aGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTc                                                                                                                          <  2:418878/151‑1 (MQ=255)                 gtgtACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTc                                                                                                                          >  1:817217/1‑150 (MQ=255)                       cGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCAc                                                                                                                   >  2:569646/1‑151 (MQ=255)                                 tGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATc                                                                                                           >  2:677874/1‑149 (MQ=255)                                      aCCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTggg                                                                                                    >  2:729513/1‑151 (MQ=255)                                         gACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACAGTCAAACCACTGAACATCAGTTGGGAt                                                                                                  >  1:861127/1‑150 (MQ=255)                                         gACGAAGGCCGAGCGGGCAAACAACGCAGCGTTTTCCACACCAACGGTTTCCTTGAGCAGGTTCCACACCAGGTCGTAGTAGATGTACGCGTTATGGTTATGCCGTCTCTGCGTATCGGAACCGTCAAACCACTGCACATCAGTTGGtata                                                                                                 <  1:229324/151‑1 (MQ=255)                                           cGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACg                                                                                               >  1:482795/1‑151 (MQ=255)                                               ggCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCg                                                                                           <  1:634943/151‑1 (MQ=255)                                                 ccGAGCGGGCAAACAAGACAGCTTCTTCCTCTCCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTGTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGcc                                                                                         <  1:851951/151‑1 (MQ=255)                                                  cGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCa                                                                                        <  1:812390/151‑1 (MQ=255)                                                     gCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAg                                                                                     <  1:929848/151‑1 (MQ=255)                                                        ggCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCg                                                                                  <  1:425222/151‑1 (MQ=255)                                                        ggCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCg                                                                                  <  1:173967/151‑1 (MQ=255)                                                         gCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCgg                                                                                 >  1:146464/1‑151 (MQ=255)                                                          cAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGt                                                                                >  1:402276/1‑151 (MQ=255)                                                             aCAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCtt                                                                             >  1:27447/1‑151 (MQ=255)                                                               aaGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTa                                                                            >  2:441195/1‑150 (MQ=255)                                                                        tcttcCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCa                                                                   <  1:427448/150‑1 (MQ=255)                                                                               cACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCAt                                                           >  1:264659/1‑151 (MQ=255)                                                                                 ccAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATcg                                                         >  1:640437/1‑151 (MQ=255)                                                                                 ccAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATc                                                          >  2:349870/1‑150 (MQ=255)                                                                                    aCGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGa                                                      <  1:212788/151‑1 (MQ=255)                                                                                         gtCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGAcc                                                    >  2:926206/1‑148 (MQ=255)                                                                                         gtCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGAcc                                                    >  1:472214/1‑148 (MQ=255)                                                                                           ccTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGAcc                                               <  2:777405/151‑1 (MQ=255)                                                                                               gAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTc                                           >  1:692101/1‑151 (MQ=255)                                                                                                 gCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAg                                         <  1:333737/151‑1 (MQ=255)                                                                                                   aCGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGtt                                       >  2:447/1‑151 (MQ=255)                                                                                                       tCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGt                                    <  2:169804/150‑1 (MQ=255)                                                                                                       tCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTc                                   >  1:771378/1‑151 (MQ=255)                                                                                                              cAGTTCGTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTAc                            >  1:208576/1‑151 (MQ=255)                                                                                                                    tttGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTg                      >  1:169414/2‑151 (MQ=255)                                                                                                                    gTTGTAGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTg                      >  2:390456/1‑151 (MQ=255)                                                                                                                         aGATGTACGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGc                 <  2:416955/151‑1 (MQ=255)                                                                                                                                cGCATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCgg           <  1:811169/150‑1 (MQ=255)                                                                                                                                  cATAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGAtt        >  1:13354/1‑151 (MQ=255)                                                                                                                                   aTAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTg       >  2:789375/1‑151 (MQ=255)                                                                                                                                    tAATGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTgg      >  1:551027/1‑151 (MQ=255)                                                                                                                                      aTGGTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGt     >  2:568942/1‑150 (MQ=255)                                                                                                                                        ggTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTaaa  <  2:262389/151‑1 (MQ=255)                                                                                                                                        ggTTATGCCTTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTaa   >  1:884709/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                               AACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCATTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTAAA  >  NC_000913/3833088‑3833373

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