Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	780,929	C→T	intergenic (+11/‑139)	lysW → / → valZ	tRNA‑Lys/tRNA‑Val

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	780,929	0	C	T	98.1%	186.6 / NA	53	intergenic (+11/‑139)	lysW/valZ	tRNA‑Lys/tRNA‑Val
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (27/25);  minor base G (1/0);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATT  >  NC_000913/780783‑781075                                                                                                                                                   |                                                                                                                                                   ggAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAg                                                                                                                                                <  2:761629/151‑1 (MQ=21) gcagagCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGa                                                                                                                                                  <  2:759881/147‑1 (MQ=18)               ttACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCg                                                                                                                                  >  2:227799/1‑151 (MQ=255)                 aCAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCCCCAGTTTTAACATTGAAGACAGATGTTAAGCgtg                                                                                                                                >  2:15289/1‑151 (MQ=255)                 aCAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCgtg                                                                                                                                <  2:738047/151‑1 (MQ=255)                               ggcggTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGGTGGTAGAGCAGTTTACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAATATTGGAAAGAGGTGTTAATCGGATAGGATAAcccag                                                                                                                   >  1:441892/1‑146 (MQ=11)                               gacggATCAATCCCATCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGc                                                                                                                  <  1:749318/149‑1 (MQ=255)                                        aTCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAAc                                                                                                         >  2:574479/1‑151 (MQ=255)                                          cccGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACgg                                                                                                       <  2:653076/151‑1 (MQ=255)                                             gTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGccc                                                                                                    <  2:776055/151‑1 (MQ=255)                                             gTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGAATAACGTTGCGTCAGCAACGGccc                                                                                                    >  2:774584/1‑151 (MQ=255)                                                 tcaCCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGt                                                                                                  >  2:737521/1‑149 (MQ=255)                                                   acccaccACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAggg                                                                                              <  2:661079/151‑1 (MQ=255)                                                            cGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGagcgaag                                                                                     <  2:689964/151‑1 (MQ=255)                                                               gTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGagcgaagcga                                                                                  <  1:445730/151‑1 (MQ=255)                                                               gTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGagcgaagcg                                                                                   >  2:520235/1‑150 (MQ=255)                                                                tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGagcgaagcg                                                                                   >  2:253122/1‑149 (MQ=255)                                                                 cGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGTCAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGGGAGGGAAGCGCGt                                                                                >  1:672058/1‑151 (MQ=255)                                                                  gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTc                                                                               >  2:418796/1‑151 (MQ=255)                                                                   ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTc                                                                               >  2:276191/1‑150 (MQ=255)                                                                        tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGGCGGATGTTAAGCGTGAAGGATAACGTTGCGTCAGCAACGGCTCGGAAGgggagcgggaagggaaatagtg                                                                         >  1:663303/1‑136 (MQ=11)                                                                         cAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTg                                                                         <  2:164481/150‑1 (MQ=255)                                                                         cAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGc                                                                        >  2:462970/1‑151 (MQ=255)                                                                          aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGc                                                                        >  2:629274/1‑150 (MQ=255)                                                                          aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGc                                                                        >  2:276765/1‑150 (MQ=255)                                                                                   agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc                                                              >  2:600179/1‑151 (MQ=255)                                                                                      aGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACt                                                           <  2:542499/151‑1 (MQ=255)                                                                                        ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTa                                                          >  2:561359/1‑150 (MQ=255)                                                                                        ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTa                                                          >  2:386377/1‑150 (MQ=255)                                                                                           aCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGa                                                      <  1:227766/151‑1 (MQ=255)                                                                                              tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACg                                                    <  1:554325/150‑1 (MQ=255)                                                                                                  aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTggg                                               <  2:72212/151‑1 (MQ=255)                                                                                                  aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTggg                                               >  1:748336/1‑151 (MQ=255)                                                                                                   tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACGTGGTACACAGAAGTTAAGCGTGTAGGAAAACGTTGCGTCAGAAACGGCCCGTAGGGCGAGCGCCGCGATTCATCCACCAGGACCTACCCCTAATGACg                                                    >  2:52949/1‑145 (MQ=14)                                                                                                      aTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCg                                           <  1:251941/151‑1 (MQ=255)                                                                                                      aTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGATc                                            >  1:308650/1‑150 (MQ=255)                                                                                                       ttGGTCGAAGGTTCGAATCCTGCTCGACCCAACAGTTTTAACATTGAAGACTGATGTTAAGCGTGTAGGATCACGTTGCGTCAGCAACGGCCCGTGGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCg                                           <  2:571131/150‑1 (MQ=17)                                                                                                            cGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGaa                                     <  1:241248/151‑1 (MQ=255)                                                                                                            cGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGa                                      >  2:159357/1‑150 (MQ=255)                                                                                                             gCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGa                                      >  2:545996/1‑149 (MQ=255)                                                                                                               aGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTCGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGCTCggag                                        >  1:696028/1‑143 (MQ=255)                                                                                                                       aTCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTgtag                          <  1:574424/151‑1 (MQ=255)                                                                                                                          cTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTgtagta                        >  2:260064/1‑150 (MQ=255)                                                                                                                             cACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCa                    <  1:161334/151‑1 (MQ=255)                                                                                                                             cACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCa                    <  2:350757/151‑1 (MQ=255)                                                                                                                                gACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAgcg                 >  2:356715/1‑151 (MQ=255)                                                                                                                                       caGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATc          <  2:369917/151‑1 (MQ=255)                                                                                                                                       caGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATc          >  2:381262/1‑151 (MQ=255)                                                                                                                                        aGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGATCTAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTGGTAGTATCCAGCGCAGTATc          <  1:561721/150‑1 (MQ=255)                                                                                                                                          ttttAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGATCGGTGGAAGAAGTTTGTAGTATCCAGCGCAGTATCggg       >  1:716763/1‑151 (MQ=255)                                                                                                                                              aaCATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGAt   <  2:135434/151‑1 (MQ=255)                                                                                                                                               aCATTGAAGGCAGATGTTAAGCGTGAAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGAtt  <  1:542838/151‑1 (MQ=255)                                                                                                                                               aCATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGAtt  <  2:314297/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATT  >  NC_000913/780783‑781075

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