Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,026,941	0	A	C	40.0%	89.8 / 10.6	50	E392D (GAA→GAC)	guaD	guanine deaminase
Reads supporting (aligned to +/- strand):  ref base A (25/5);  new base C (0/20);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 2.25e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAAT  >  NC_000913/3026791‑3027089                                                                                                                                                       |                                                                                                                                                     aGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGaa                                                                                                                                                      <  1:386657/151‑1 (MQ=255)      aCAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCAc                                                                                                                                                 >  1:550135/1‑151 (MQ=255)        aaaGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCg                                                                                                                                               <  1:662316/151‑1 (MQ=255)         aaGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGc                                                                                                                                              >  2:554026/1‑151 (MQ=255)          aGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGcc                                                                                                                                             >  1:680788/1‑151 (MQ=255)            tATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCAc                                                                                                                                           >  1:244350/1‑151 (MQ=255)              ttGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTc                                                                                                                                         >  2:81871/1‑151 (MQ=255)              ttGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTc                                                                                                                                         >  2:464293/1‑151 (MQ=255)              ttGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTc                                                                                                                                         >  1:526083/1‑151 (MQ=255)                  aaTTACAAGCCTATCGCCTCCCGGCTTATACGTCGTTTTACCGGCCCCCGCTCGGCGGAGCGAAATCTCTGGCCCTTGTCGTTTTTCTTGGCCATTTTTTCCCGTGAAAGGAGGCTGTTTCCGTGGTGATGGACCCCACCGCCACTCCGCt                                                                                                                                     <  2:779637/151‑1 (MQ=255)                   aTTACAGGGCTACCGCCTCTCGCAATAGGAATCTTTTTCCCTGGCCCCGCTCGGCGGAGCGAAATCTCGGGGCTTTGACGTTTTGTTTGGCACTTTTTTCCCGGGCAAAGAGGTTGATTTCGTGGTGAGGAAACCCACCGCCACTcc                                                                                                                                        <  2:57294/147‑1 (MQ=255)                    ttACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTAc                                                                                                                                   >  2:865809/1‑151 (MQ=255)                       cAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTAcagc                                                                                                                                >  1:729719/1‑151 (MQ=255)                        aaGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTAcagc                                                                                                                                >  1:704153/1‑150 (MQ=255)                        aaGGCTACCGCCTCTCGGCATATGAAGCTTTTTCCCTGGCCCCGCCCGGGGGGGCGAAATCTCTGGGCCTTGACGATTGGATTGGCAACTTTTTACCGGGAAAAGAGGCTGATTTCGGGGTGATGGAACCCACCGCCACTCCGCTAcag                                                                                                                                 <  2:35648/149‑1 (MQ=255)                         aGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTAcagc                                                                                                                                >  1:347956/1‑149 (MQ=255)                           gCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCt                                                                                                                            >  1:665421/1‑151 (MQ=255)                             tATCGCCTCTCGGCATATGAAGCGTTTTCCCTGGCCACGCTCGGGGGAGCGAAATCTCGGGGCCTTGACGATTTGATTGGCAACTTTTTCCCTGGCAAAGAGGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCt                                                                                                                            <  1:480386/149‑1 (MQ=255)                                  ccTCGCGGCATATGAAGCTTTTTACCGGGCCACGCTCGGGGGAGCAAAATCTCTGGGCCTTGACGTTTTGATTGGAAACTTTTTCCCGGGAAAAGAGGCTGTTTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTa                                                                                                                      <  1:929419/150‑1 (MQ=255)                                        ggCATAGGAGGCGTTTTACCTGGCCCCGCTCGGGGGGGCGAAATCTCGGGGCCTTGACGATTTGATTGGCACCTTTTTCCCTGGAAAAGAGGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACa                                                                                                                 <  2:44143/149‑1 (MQ=255)                                            tatGAAGCTTTTTACCGGGCCCCGCTCGGGGGAGGAAAATCTCGGGGCCTTGCCGTTTTGATTGGCAACTTTTTCCCTGGAAAAGAGGCTGTTTTCGGGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAActct                                                                                                            <  2:859216/150‑1 (MQ=255)                                                    gTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCttt                                                                                                    >  2:730214/1‑150 (MQ=255)                                                    gTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCttt                                                                                                    >  1:706100/1‑150 (MQ=255)                                                     ttttcccgtgCCCCGCTGTGCGGGGCAAAATCTCGGGGCTTTGACTTTTTGATGGCCACCTTTTTCCCGGGCAAGGGGGCTTGTTTGGTGGTGAGGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCt                                                                                                      <  2:1483/138‑1 (MQ=255)                                                       ttcccTGGCCCCGCCCGGCGAGGCGAACTCTCTGGCCCTTTCCTTTTTTATGGACACCTTTTTCCCTGGAAAAGAGGCTGATTTCGTGGTGAGGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGt                                                                                                 <  2:845273/147‑1 (MQ=255)                                                          ccTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGAc                                                                                             >  1:778799/1‑151 (MQ=255)                                                               ccgcGCTCGGGGGGGCGAAATCCCGGGGCTTTGACGATTGGCTCGGCAATTTTTTCCCGGGAAAAGAGCCTGATTTCGTGGTGATGGAACCCTCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAAt                                                                                         <  1:98252/147‑1 (MQ=255)                                                               ccACTCTCGGCGGAGCAAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTTCCGGACAAAGAGGCTGTTTTCGTGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAAtt                                                                                        <  2:940801/151‑1 (MQ=255)                                                               ccACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAAtt                                                                                        >  2:844440/1‑151 (MQ=255)                                                                cACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGNNNNGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACNGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTg                                                                                       >  1:517521/1‑151 (MQ=255)                                                                    cgggggggggcaaaaTCTCGGGCCCTTGACGTTTTGATTGGCAACTTTTTCCCGGAAAAAGAGGCTGATTTCGTGGTGAGGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGt                                                                                      <  2:75032/136‑1 (MQ=255)                                                                     tCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGCCAAATTGTTCGt                                                                                  >  2:641417/1‑151 (MQ=255)                                                                       ggcggAGCGAAATCTCTGGCCCTTGACGTTTTGATTGGCAACTTTTCCCCGGCCAAAGGGCCTGATTTCTTGTTGATGGACCCCACCTCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGt                                                                                  <  2:35505/149‑1 (MQ=255)                                                                            aGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatgatg                                                                           >  2:373565/1‑151 (MQ=255)                                                                            aGCAAAATCTCTGGGCCTTGACGATTTGATTGGCATCTTTTTCCCTGGCAAAGAGGTTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatgat                                                                            <  1:693127/150‑1 (MQ=255)                                                                                aaaTCTCGGGGCCTTGACGTTTTGATTGGCAATTTTTTCCCGGGAAAAGAGGTTGATTTCGTGGGGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGt                                                                       <  1:754252/151‑1 (MQ=255)                                                                                      cggggCCTTGACGATTTGATTGGAAACTTTTTACCTGAAAAAGAGGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTggg                                                                   <  2:556724/147‑1 (MQ=255)                                                                                             ttGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCg                                                          >  1:471339/1‑151 (MQ=255)                                                                                                   aTTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGAt                                                    >  1:174487/1‑151 (MQ=255)                                                                                                      tGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCt                                                  >  1:652585/1‑150 (MQ=255)                                                                                                       gATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTAc                                                >  2:656579/1‑151 (MQ=255)                                                                                                         ttGGCATTTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTa                                                 <  1:9314/148‑1 (MQ=255)                                                                                                             cAATTTTTTACCTGGAAAAGAGGGTGATTTCGTGGTGATGAACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCg                                              <  2:558982/147‑1 (MQ=255)                                                                                                                       ccgggCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATg                                <  1:722288/148‑1 (MQ=255)                                                                                                                                        tttCGTGGTGATGGACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGgtgt                    <  1:396983/146‑1 (MQ=255)                                                                                                                                            gtggtgATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGc            >  1:929476/1‑150 (MQ=255)                                                                                                                                              ggtgATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAAc              <  1:464468/146‑1 (MQ=255)                                                                                                                                              ggtgATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAAc         >  1:460265/1‑151 (MQ=255)                                                                                                                                                tgtgggACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACg             <  2:549940/141‑1 (MQ=255)                                                                                                                                                    ggACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACtaat     <  2:665196/149‑1 (MQ=255)                                                                                                                                                    ggACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACtaa      <  1:247302/148‑1 (MQ=255)                                                                                                                                                      aaCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACtaataat  >  2:763544/1‑150 (MQ=255)                                                                                                                                                       aCCCCCCGCCAATCCGCTACACCAGCTGCGCTATGACAACTCTGTTTCCTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGACCTACCGCCCCTACGTTGATGGTCGTCTAGTGTACGAACGCAACtaata    <  2:815859/147‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                     AGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAAT  >  NC_000913/3026791‑3027089

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