Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,876,744	T→C	N70S (AAC→AGC)	yidA ←	sugar phosphatase YidA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,876,744	0	T	C	97.8%	158.1 / NA	46	N70S (AAC→AGC)	yidA	sugar phosphatase YidA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (21/24);  minor base G (1/0);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGTTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAGTTGTGCACACCTGCATACGGGCGACCCGTCGTTAGCACGACATTCACGCCACGGGCGCGAGCTGCGGCAATCGCATTTTTA  >  NC_000913/3876607‑3876878                                                                                                                                          |                                                                                                                                       aCAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTATGCGCCACAGTGTTACAATCAGCGGCCTTCTGTACCAGCGCGCCGCta                                                                                                                                      >  2:334362/1‑140 (MQ=255)     cGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAg                                                                                                                                  >  2:94054/1‑140 (MQ=255)     cGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAg                                                                                                                                  >  1:87187/1‑140 (MQ=255)     cGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAg                                                                                                                                  >  1:699093/1‑140 (MQ=255)          tgCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATg                                                                                                                             <  1:794491/140‑1 (MQ=255)                   gggCGTGGAAATGAGAACCGACCTCCCGATAGAGTTTCTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCg                                                                                                                    <  1:253184/140‑1 (MQ=255)                    ggCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGc                                                                                                                   <  2:724635/140‑1 (MQ=255)                     gCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGcc                                                                                                                  <  2:676729/140‑1 (MQ=255)                         ggAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGc                                                                                                              >  2:499179/1‑140 (MQ=255)                         ggAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGc                                                                                                              >  1:835943/1‑140 (MQ=255)                         ggAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGACGCTATAAGTAATTCAGTAGTCGCCCGGc                                                                                                              >  1:412054/1‑140 (MQ=255)                          gAAATGAGCACCGACTTCGCGAAAGAGTTTTTCGATGAATCGATAGTCCTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCt                                                                                                             <  2:412054/140‑1 (MQ=255)                          gAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCt                                                                                                             <  2:835943/140‑1 (MQ=255)                          gAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCt                                                                                                             <  2:164622/140‑1 (MQ=255)                            aaTGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGt                                                                                                           <  1:696638/140‑1 (MQ=255)                             aTGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGtt                                                                                                          <  1:94054/140‑1 (MQ=255)                               gagaACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCACCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTcc                                                                                                        <  2:13626/140‑1 (MQ=255)                                  aaCCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGAACCAGCGCACCGGTATTAGTAATGCAGTAGACGCCCGGCTGTTCCata                                                                                                     >  2:301132/1‑140 (MQ=255)                                       aCTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAg                                                                                                <  2:92600/140‑1 (MQ=255)                                       aCTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAg                                                                                                <  2:710920/140‑1 (MQ=255)                                          tCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGctc                                                                                             <  1:499179/140‑1 (MQ=255)                                           cgcgAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGctct                                                                                            >  1:172532/1‑140 (MQ=255)                                                    gTTTTTCCAGGAAACTATAGTCGTCTTAGTTGAGAGCTGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAg                                                                                   <  1:301132/140‑1 (MQ=255)                                                    gTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAg                                                                                   <  1:53074/140‑1 (MQ=255)                                                    gTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAg                                                                                   <  1:497708/140‑1 (MQ=255)                                                    gTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAg                                                                                   <  1:190750/140‑1 (MQ=255)                                                          ccAGGAAACGATAGTCGGCATAACTGAGAGCTGTTGGCGTCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAGTTGTGc                                                                             <  2:294066/140‑1 (MQ=255)                                                           cAGGAAACGATAGTCGTTATAGCTGAGAGCATTTTTCGCCACGGTGCTTCCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGTAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAGTTGTGca                                                                            <  2:346628/140‑1 (MQ=255)                                                             ggAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAGTTGTGcaca                                                                          >  1:179851/1‑140 (MQ=255)                                                              gAAACGATAGTCGTCCTAGCTGAGACCAGTTTGCGCCGCGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAGTTGTGCacac                                                                         <  1:311363/140‑1 (MQ=255)                                                              gAAACGATAGTCGTCATCGCTGGGAGCAGTTTGCGGCACAGTGCTACCCTCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAGTTGTGCacac                                                                         <  2:276825/140‑1 (MQ=255)                                                              gAAACGATAGTCGTCATAGCTGTGAGCAGCTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAGTTGTGCacac                                                                         <  1:334362/140‑1 (MQ=255)                                                              gAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAGTTGTGCacac                                                                         <  2:50981/140‑1 (MQ=255)                                                               aaaCGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGTAGTAGTCGCCCGGCAGTTCCATCTGCTGCACTTTCAAGTAGTTGTGCACAcc                                                                        >  1:346628/1‑140 (MQ=255)                                                               aaaCGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAGTTGTGCACAcc                                                                        >  2:151028/1‑140 (MQ=255)                                                                                      agCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAGTTGTGCACACCTGCATACGGGCGACCCGTCGTTa                                                 >  2:49975/1‑140 (MQ=255)                                                                                      agCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGATAGTTGTGCACGCCAGCATACAGGCGACCCGTCGTTa                                                 >  1:401723/1‑140 (MQ=255)                                                                                           tttGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAGTTGTGCACACCTGCATACGGGCGACCCGTCGTTAGCACg                                            <  2:699093/140‑1 (MQ=255)                                                                                                                 aGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAGTTGTGCACACCTGCATACGGGCGACCCGTCGTTAGCACGACATTCACGCCACGGGCGCGAg                      >  2:591784/1‑140 (MQ=255)                                                                                                                 aGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAGTTGTGCACACCTGCATACGGGCGACCCGTCGTTAGCACGACATTCACGCCACGGGCGCGAg                      >  1:202491/1‑140 (MQ=255)                                                                                                                 aGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAGTTGTGCACACCTGCATACGGGCGACCCGTCGTTAGCACGACATTCACGCCACGGGCGCGAg                      >  1:689713/1‑140 (MQ=255)                                                                                                                 aGCGGCCTTCTGTACCAGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAGTTGTACACACCTGCATACGGGCGACCCGTCGTTAGCACGACATTCACGCCACGGGCGCGAg                      >  2:425386/1‑140 (MQ=255)                                                                                                                                 aGCGCGCCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAGTTGTGCACACCTGCATACGGGCGACCCGTCGTTAGCACGCCATTCACGCCACGGGCGCGAGCTGCGGCAATCGCAtt      >  1:385331/1‑140 (MQ=255)                                                                                                                                  gcgcgcCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAGTTGTGCACACCTGCATACGGGCGACCCGTCGTTAGCACGACATTCACGCCACGGGCGCGAGCTGCGGCAATCGCAttt     >  1:623339/1‑140 (MQ=255)                                                                                                                                  gcgcgcCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAGTTGTGCACACCTGCATACGGGCGACCCGTCGTTAGCACGACATTCACGCCACGGGCGCGAGCTGCGGCAATCGCAttt     >  1:752799/1‑140 (MQ=255)                                                                                                                                   cgcgcCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAGTTGTGCACACCTGCATACGGGCGACCCGTCGTTAGCACGACATTCACGCCACGGGCGCGAGCTGCGGCAATCGCAtttt    >  1:108999/1‑140 (MQ=255)                                                                                                                                   cgcgcCGCTATAAGTAATGCAGTAGTCGCCCGGCTGTTCAATATGCAGCTCTTTCAGGTAGTTGTGCACACCTGCATACGGGCGACCCGTCGTTAGCACGACATTCACGCCACGGGCGCGAGCTGCGGCAATCGCCtttt    >  2:150438/1‑140 (MQ=255)                                                                                                                                     cgccgcTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAGTTGTGCACACCTGCATACGGGCGACCCGTCGTTAGCACGACATTCACGCCACGGGCGCGAGCTGCGGCAATCGCATTTTTa  >  2:728023/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                       ACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTGTACCAGCGCGCCGTTATAAGTAATGCAGTAGTCGCCCGGCTGTTCCATATGCAGCTCTTTCAGGTAGTTGTGCACACCTGCATACGGGCGACCCGTCGTTAGCACGACATTCACGCCACGGGCGCGAGCTGCGGCAATCGCATTTTTA  >  NC_000913/3876607‑3876878

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: ‑