Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,193,178	A→C	V105V (GTA→GTC)	yqiK →	protein YqiK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,193,178	0	A	C	96.2%	176.2 / NA	53	V105V (GTA→GTC)	yqiK	protein YqiK
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base C (31/19);  minor base G (2/0);  total (34/19)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCATGAAATAATCCCCATCAATATGAATACTCTGAAGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTAAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGACTTACGTATGTTGGTTGAAGATAAATTTGTCGATGCCCTCCGTGCAACAGCTGCGCAA  >  NC_000913/3193053‑3193304                                                                                                                              |                                                                                                                               tttCATGAAATAATCCCCATCAATATGAATACTCTGAAGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCtt                                                                                                                         <  1:550689/133‑1 (MQ=255)                 ccATCAATCTGAATATTCGGAACCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGc                                                                                                         <  1:195479/133‑1 (MQ=255)                 ccATCAATATGAATACTCTGAAGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGc                                                                                                         <  2:406170/133‑1 (MQ=255)                      aaTATGAATACTCTGAAGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCg                                                                                                    <  1:336503/133‑1 (MQ=255)                        tatGAATACTCTGAAGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCt                                                                                                  >  2:35971/1‑133 (MQ=255)                               aCTCTGAAGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGa                                                                                           >  2:128612/1‑133 (MQ=255)                                ctctGAAGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGAc                                                                                          >  2:399858/1‑133 (MQ=255)                                  ctGAAGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGc                                                                                        >  2:350717/1‑133 (MQ=255)                                  ctGAAGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGc                                                                                        >  2:31656/1‑133 (MQ=255)                                  ctGAAGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGc                                                                                        >  1:342057/1‑133 (MQ=255)                                   tGAAGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTGAAACCTTCGGTAGAAGGGATTGCCAACGCTTCCCAGACGCt                                                                                       >  2:196365/1‑133 (MQ=255)                                   tGAAGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTGAAACCTTCAGTAGAAGGGATTGCCACCGCTACCCAGACGCt                                                                                       >  2:207560/1‑133 (MQ=255)                                   tGAAGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTAAAACCTTCGGTAAAAGCGATTGCCAACGCTGCCCAGACGCt                                                                                       >  2:191010/1‑133 (MQ=255)                                    gAAGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTg                                                                                      <  1:31656/133‑1 (MQ=255)                                     aaGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTgg                                                                                     >  2:483227/1‑133 (MQ=255)                                     aaGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTgg                                                                                     >  2:404922/1‑133 (MQ=255)                                     aaGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGAAGAGGGGGTTGCCACCTCTGCCCAGACGCTgg                                                                                     >  2:235605/1‑133 (MQ=255)                                     aaGCTGGAAGTCAGCCGCTCAACCATTAATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTgg                                                                                     >  2:543803/1‑133 (MQ=255)                                      aGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGGAGGGATTTCCAGCGGTGCCCAGACGCTggg                                                                                    >  1:180753/1‑133 (MQ=255)                                      aGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTggg                                                                                    >  1:441525/1‑133 (MQ=255)                                      aGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTggg                                                                                    <  1:35971/133‑1 (MQ=255)                                           gAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTTAGTAGAAGGGATGGACTCCGCTGCCTAGACGTTGGGGCAAc                                                                               >  2:271021/1‑133 (MQ=255)                                           gAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAAc                                                                               >  2:67676/1‑133 (MQ=255)                                           gAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAAc                                                                               >  2:43252/1‑133 (MQ=255)                                                    cGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGt                                                                      <  1:67676/133‑1 (MQ=255)                                                          aCCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTg                                                                <  1:483227/133‑1 (MQ=255)                                                              ttGATAGCCTGATTATGACAGAACGTATGCACGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGa                                                            <  1:196365/133‑1 (MQ=255)                                                              ttGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGa                                                            <  1:350717/133‑1 (MQ=255)                                                                   aGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGACTTAc                                                       >  2:521462/1‑133 (MQ=255)                                                                     ccTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGTACCCTGTCCCCTGTAGACTTACGt                                                     >  2:254900/1‑133 (MQ=255)                                                                     ccTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGaa                                                              <  1:524291/124‑1 (MQ=255)                                                                     ccTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGaa                                                              >  2:524291/1‑124 (MQ=255)                                                                     ccTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGACTTACGt                                                     >  2:389343/1‑133 (MQ=255)                                                                     ccTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGACTTACGt                                                     >  1:17251/1‑133 (MQ=255)                                                                     ccTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGACTTACGt                                                     >  1:114751/1‑133 (MQ=255)                                                                         aTTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGACTTACGTATgt                                                 >  1:386443/1‑133 (MQ=255)                                                                         aTTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGACTTACGTATgt                                                 >  2:602173/1‑133 (MQ=255)                                                                            aCGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGACTTACGTATgttgg                                              <  1:43252/133‑1 (MQ=255)                                                                             cGAAAGATCGCATGCCCGTCGATGCAGTCGTCACTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGACTTACGTATgttggt                                             <  2:166523/133‑1 (MQ=255)                                                                              gAAAGATCTTATGCGTGTGGATTTAGTAGTTGCTTTCTTTGGGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGACTTACGTATgttggtt                                            <  1:235605/133‑1 (MQ=255)                                                                              gAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGACTTACGTATgttggtt                                            >  1:461007/1‑133 (MQ=255)                                                                              gAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGACTTACGTATgttggtt                                            <  1:404922/133‑1 (MQ=255)                                                                                   aTCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGACTTACGTATGTTGGTTGAAGa                                       >  1:432941/1‑133 (MQ=255)                                                                                   aTCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGACTAACGTATGTTGGTTGAAGa                                       >  1:166523/1‑133 (MQ=255)                                                                                   aTCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTAGCCTGAATACTCACCTATGTTGGTTGAAGa                                       >  1:194575/1‑133 (MQ=255)                                                                                           cgcgTCGATGTAGTAATCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGACTTACGTATGTTGGTTGAAGATAAATTTg                               <  1:254900/133‑1 (MQ=255)                                                                                            gcgTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGACTTACGTATGTTGGTTGAAGATAAATTTGt                              >  2:145425/1‑133 (MQ=255)                                                                                               tCGATGTAGTAGTCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGACTTACGTATGTTGGTTGAAGATAAATTTGTCGa                           >  2:436837/1‑133 (MQ=255)                                                                                                      agtagtCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGACTTACGTATGTTGGTTGAAGATAAATTTGTCGATGCCCTc                    <  2:505290/133‑1 (MQ=255)                                                                                                       gtagtCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGGCTTACGTATGTTGGTTAAAGATAAAGTTGTAGATGCCCTcc                   >  2:197919/1‑133 (MQ=255)                                                                                                        tagtCGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGACTTACGTATGTTGGTTGAAGATAAATTTGTCGATGCCCTCCg                  <  2:17251/133‑1 (MQ=255)                                                                                                            cGCTTTCTTTGTGCGGGTCAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGACTTAGGTATGTTGTTTGGAGATAAATTTGTCGATGTCCTCCGTGCa              >  2:233381/1‑133 (MQ=255)                                                                                                                        gCGGGTCACACCTTCCGTCGAATGGATTGCTACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGACTTACGTATGTTGGTTGAAGATAAATTTGTCGATGCCCTCCGTGCAACAGCTGCGCaa  <  1:191010/133‑1 (MQ=255)                                                                                                                              |                                                                                                                               TTTCATGAAATAATCCCCATCAATATGAATACTCTGAAGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTAAAACCTTCAGTAGAAGGGATTGCCACCGCTGCCCAGACGCTGGGGCAACGCACCCTGTCGCCTGAAGACTTACGTATGTTGGTTGAAGATAAATTTGTCGATGCCCTCCGTGCAACAGCTGCGCAA  >  NC_000913/3193053‑3193304

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