Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	184,723	G→A	A116V (GCC→GTC)	yaeI ←	phosphodiesterase YaeI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	184,723	0	G	A	98.4%	206.7 / ‑4.2	61	A116V (GCC→GTC)	yaeI	phosphodiesterase YaeI
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (23/37);  total (24/37)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

TTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGGCAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCC  >  NC_000913/184590‑184847                                                                                                                                      |                                                                                                                             ttCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACg                                                                                                                        >  2:458574/1‑140 (MQ=255)  tCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGt                                                                                                                       <  1:401819/140‑1 (MQ=255)   cGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTc                                                                                                                      >  1:194651/1‑140 (MQ=255)      aTTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGgc                                                                                                                   >  1:12586/1‑140 (MQ=255)       ttGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAGACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGgcg                                                                                                                  >  2:882980/1‑140 (MQ=255)          ccTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGcac                                                                                                               >  2:836406/1‑140 (MQ=255)              ttCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTc                                                                                                           >  1:730346/1‑140 (MQ=255)              ttCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCAGAAGCAGACAAACGTCGGCGCACATTc                                                                                                           <  2:526021/140‑1 (MQ=255)                cGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCgg                                                                                                         <  2:359379/140‑1 (MQ=255)                 ggCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGc                                                                                                        <  2:730346/140‑1 (MQ=255)                  gCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCa                                                                                                       <  1:458574/140‑1 (MQ=255)                   cGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCaa                                                                                                      <  2:770838/140‑1 (MQ=255)                     tGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAgg                                                                                                    <  2:414039/140‑1 (MQ=255)                       gCGATCACCGTAGCCTGGTTAGACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAgggg                                                                                                  <  1:882980/140‑1 (MQ=255)                        cGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAggggg                                                                                                 <  2:417094/140‑1 (MQ=255)                          aTCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGGCAAACGTCGGCGCACATTCGGCAAGGGGGga                                                                                               >  1:653769/1‑140 (MQ=255)                               cGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTa                                                                                          <  2:760269/140‑1 (MQ=255)                                 tAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACg                                                                                        >  1:253745/1‑140 (MQ=255)                                 tAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACg                                                                                        >  1:540151/1‑140 (MQ=255)                                      tGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACt                                                                                   >  2:459756/1‑140 (MQ=255)                                         ttAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTaaa                                                                                >  1:621830/1‑140 (MQ=255)                                         ttAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTaaa                                                                                >  2:128779/1‑140 (MQ=255)                                           aaaCAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGga                                                                                               <  2:670694/123‑1 (MQ=255)                                           aaaCAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGga                                                                                               >  1:670694/1‑123 (MQ=255)                                              caacaCCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCg                                                                           >  2:465069/1‑140 (MQ=255)                                                acacCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGaa                                                                         <  1:542542/140‑1 (MQ=255)                                                  acCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGaaaa                                                                       >  1:775510/1‑140 (MQ=255)                                                     gcgATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAAtt                                                                    <  2:621830/138‑1 (MQ=255)                                                          gCCCGCTGATTTCAGCGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCg                                                               <  2:189716/140‑1 (MQ=255)                                                          gCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCg                                                               <  2:194651/140‑1 (MQ=255)                                                          gCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCg                                                               <  1:465069/140‑1 (MQ=255)                                                          gCCCGCTGATTTCAACGTCTCGCCAATTAAGTCATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGCGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCg                                                               <  2:810192/140‑1 (MQ=255)                                                           cccGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGa                                                              >  1:565374/1‑140 (MQ=255)                                                            ccGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGAc                                                             <  1:832140/140‑1 (MQ=255)                                                              gCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACat                                                           >  1:858077/1‑140 (MQ=255)                                                                tGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACatat                                                         >  2:651596/1‑140 (MQ=255)                                                                   tttCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCaa                                                      >  1:117482/1‑140 (MQ=255)                                                                    ttCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCaaa                                                     <  2:297897/140‑1 (MQ=255)                                                                    ttCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCaaa                                                     <  1:619711/140‑1 (MQ=255)                                                                    ttCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCaaa                                                     <  1:836406/140‑1 (MQ=255)                                                                     tCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAaat                                                    >  1:419837/1‑140 (MQ=255)                                                                     tCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAaat                                                    <  2:599913/140‑1 (MQ=255)                                                                       aaCGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAaataa                                                  <  2:221591/140‑1 (MQ=255)                                                                           tctcGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAAATAATACg                                              <  2:12586/140‑1 (MQ=255)                                                                            ctcGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAAATAATACGt                                             >  1:623482/1‑140 (MQ=255)                                                                              cGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAAATAATACGTaa                                           >  1:200508/1‑140 (MQ=255)                                                                                 cAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAAATAATACGTAATcg                                        <  2:419837/140‑1 (MQ=255)                                                                                    ttAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAAATAATACGTAATcgccg                                     <  2:540151/140‑1 (MQ=255)                                                                                     tAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAAATAATACGTAATcgccgc                                    <  2:117482/140‑1 (MQ=255)                                                                                     tAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAAATAATACGTAATcgccgc                                    <  1:651596/140‑1 (MQ=255)                                                                                     tAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAAATAATACGTAATcgccgc                                    <  2:309708/140‑1 (MQ=255)                                                                                         tGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAAATAATACGTAATCgccgcc                                   >  1:157782/1‑137 (MQ=255)                                                                                         tGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAAATAATACGTAATCgccgcc                                   <  2:157782/137‑1 (MQ=255)                                                                                         agATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAg                                <  2:215906/139‑1 (MQ=255)                                                                                          gATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGt                               <  1:675864/140‑1 (MQ=255)                                                                                           aTTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTa                              >  1:168050/1‑140 (MQ=255)                                                                                              tttttCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAata                           <  1:427121/140‑1 (MQ=255)                                                                                               ttttCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAatat                          <  2:168050/140‑1 (MQ=255)                                                                                                     gTACCAACAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATc                    <  1:646614/140‑1 (MQ=255)                                                                                                           aCAGGGCGATCGTGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCtt              <  1:262366/140‑1 (MQ=255)                                                                                                                       tGGTTGCCGAAGCAGACAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGcc  <  2:764262/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                             TTCGAATTGCCTGTTCGGCGTGGCGATCACCGTAGCCTGGTTAAACAACACCGTGATGCCCGCTGATTTCAACGTCTCGCCAATTAAGTGATTTTTTTCTGTACCAACAGGGCGATCGTGGTTGCCGAAGCAGGCAAACGTCGGCGCACATTCGGCAAGGGGGGAGAGTACGTCACTAAACGCCGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCC  >  NC_000913/184590‑184847

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