Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,177,748	0	C	A	43.4%	49.3 / 31.2	55	E55* (GAG→TAG)	nudF	ADP‑ribose pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (19/11);  new base A (0/23);  total (20/35)
Fisher's exact test for biased strand distribution p-value = 8.36e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGT  >  NC_000913/3177598‑3177894                                                                                                                                                       |                                                                                                                                                   cacTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTc                                                                                                                                                    >  1:83260/1‑151 (MQ=255)   cTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCa                                                                                                                                                   <  1:739534/150‑1 (MQ=255)     ttCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCaaa                                                                                                                                                 <  2:69320/150‑1 (MQ=255)         cccTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAAttt                                                                                                                                            >  1:597846/1‑151 (MQ=255)          ccTCGTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTc                                                                                                                                           <  1:354017/151‑1 (MQ=255)                cAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgc                                                                                                                                       >  2:675043/1‑149 (MQ=255)                 aaTCACCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTAAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgccg                                                                                                                                     <  2:86696/150‑1 (MQ=255)                                   ctcCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGc                                                                                                                  >  1:429167/1‑151 (MQ=255)                                       agaaGCCAGGGGGTTGCGCGGGGGCCGTACGCGCCAATCCGAATCTGCTCAGTCAGCCAACCTCCATCACGCTCTGGGTCAAAGGGTAGCAAGCCCGCGGCGTGCCCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCg                                                                                                               <  1:839579/147‑1 (MQ=255)                                         tAGCCAGGGGTTTTCGTTGGTGTCGTACGCGGCAACCCGAATCTGCTAAATCAGCACACCTTCATCCCGCACGGGGTAAAAGGGTAGCAAGACGGCGGCGTGACCGCGCTCAAAATTTTCCCGCCGCACCTCATGACTCATTTg                                                                                                                   <  2:562254/144‑1 (MQ=255)                                            ccAGGGGGTTCCGCTGGTGTCGTCCGCGGCAATCCGAATCTGCTCAATCAGCACAACTCCACCACGCACTGGGTCAAAGGGTAGAAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAAt                                                                                                         <  2:262923/151‑1 (MQ=255)                                              agggggTGTCGCGGGGGCCGTACGCGGCATCCCAAATCTCCTCAATCAGCCAAACTGCACCACGCATTGGGACAAAGGGTACCAAGACGGCGGCGTGACCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAAt                                                                                                         <  2:597912/147‑1 (MQ=255)                                              aaGGGTTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCCCGCACTGGGTAAAAGGGTAGCAAGACGGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAg                                                                                                       <  2:156874/151‑1 (MQ=255)                                                        gCTGGTGTCCTACGCGGCATCGAAAATCTCACCAAACAGCACAATTACACCCCGCTCGGGGTGAAGGGGTAGCATTACTCTGGCGTGCCCGCCCTCGAAAATTTACCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACTGTGACg                                                                                              <  1:276212/150‑1 (MQ=255)                                                          tGGTGTCGTCCGCGGCAATCCGAATCTGCCAAACCAGAACAACTTCACCACGCACGGGGTCAAAGGGTAGAAAAACGGCGGCGTGACCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGa                                                                                             <  1:779196/149‑1 (MQ=255)                                                          tGGTGTCGTACGCGGCAATCCGAATCGGCTCAATCAGCACAACTTCATCACGCACTGGGACAAAGGGGAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaa                                                                                            >  1:656464/1‑150 (MQ=255)                                                           ggtgcctaaCGCGCCACTCCGAATCTGCTCAATCACCAAACCTTCACCCCGTCCGGGGTAAAGGGGTAGAAAAACGGGGGGCTGACCGCGATAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGa                                                                                             <  1:658187/140‑1 (MQ=255)                                                           ggTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTCCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa                                                                                           <  2:788314/150‑1 (MQ=255)                                                           ggTGTCGTACGCGGAAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAGGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa                                                                                           <  2:737283/150‑1 (MQ=255)                                                             tgtCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCt                                                                                        >  1:369105/1‑151 (MQ=255)                                                             tgtCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACATCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCt                                                                                        >  1:489216/1‑151 (MQ=255)                                                               tCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtat                                                                                      >  1:698683/1‑151 (MQ=255)                                                                cGTACGCGGCAATCCGAATCTGCTCAATCACCACACCTTCACCCCGCACGGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCta                                                                                       <  2:892734/149‑1 (MQ=255)                                                                cGTACGCGGCAACCCAAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGAAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCta                                                                                       <  2:826243/149‑1 (MQ=255)                                                                cGTACGCGCCAATCCAATTCTCCTCAATCACCACAACTCCACCCCCCCCGGGGACAAAGGTTAGCAAGACTGCGGCGTGCCCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtata                                                                                     <  1:309229/151‑1 (MQ=255)                                                                 gTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtata                                                                                     >  1:474791/1‑150 (MQ=255)                                                                      cgGCAATCCGAATCTGCTAAATCACCACAATTTCAGCACGCACTGGGTCAAAGGGTAGCAAGCCTGCGGCGTGCCCGCGCTCAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATc                                                                               <  2:826856/151‑1 (MQ=255)                                                                       ggaaaTCCAAATCTGCTCAATCAGCACACCTTCATCACGCACGGGGTCAAAGGTTAGCAAGACTGCGGGGTGACCGCGCTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGAt                                                                                <  1:183442/146‑1 (MQ=255)                                                                          attCCGATTCTGCCCAACCACCACACCTTCACCACGCACTGGGTAAAAGGGTACCAAGCCGGCGCCGTGCCCGCGCTAAAAAATTGCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATaaa                                                                                   <  1:823921/141‑1 (MQ=255)                                                                           atcgaaaTCCGCTCAATCAGCACAACTCCATCCCGCGTGGGGTCAAAGGGTACCAAGATTGCGGGTTGACCGGTCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACg                                                                                                    <  2:207635/120‑1 (MQ=255)                                                                           aTCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACg                                                                                                    >  1:207595/1‑125 (MQ=255)                                                                           aTCCGAATCTGCTCAATCAGCAAACCTTCATCCCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATg                                                                          <  2:886532/151‑1 (MQ=255)                                                                               gAATCTGCGCAATCAGCACAACTTCATCACGCACTAGGTCAAAGGGTAGCAAGACTGCGGCGTGACCCCGCTTAAAAATTTCCCGCCGCACCTCATGAGTCATTTGCCCGTTGAATAGACGATGACGAAATCTAGAACGATCTAATGaaa                                                                       >  1:67314/1‑150 (MQ=255)                                                                                  tCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaaa                                                                    >  1:751734/1‑150 (MQ=255)                                                                                     gCTCAATCAGCACACCTACATCCCGCACTGGTTCAAAGGGTAGCAACACTGCGGCGTGCCCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAg                                                                  <  2:791368/149‑1 (MQ=255)                                                                                                  aaCTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGtt                                                   >  2:858256/1‑151 (MQ=255)                                                                                                   aCTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGttt                                                  >  1:388445/1‑151 (MQ=255)                                                                                                   aCTTCACCCCGCACGGGGTCAAAGGGAAGAAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGtt                                                   <  1:280996/150‑1 (MQ=255)                                                                                                     ttCATCGCGCCTTGGCTGAAAGGGCAGCAAGCCTGCGGCGTGTCCGCGCTGAAAAATGTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGttt                                                  <  2:774448/149‑1 (MQ=255)                                                                                                     ttCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtct                                                >  2:952511/1‑151 (MQ=255)                                                                                                       caccaCGCACTGGGTCAAAGGTTAGCAAGACGGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtc                                                 <  2:138618/145‑1 (MQ=255)                                                                                                        atcaCACACTGGGTCAAAGGGTAGAAAGACTGCGGCGTGCCCGCGTTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtc                                                 <  1:770884/147‑1 (MQ=255)                                                                                                         tcaCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTg                                            >  1:448291/1‑151 (MQ=255)                                                                                                          caCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGc                                           >  1:831211/1‑151 (MQ=255)                                                                                                          caCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGc                                           >  1:243742/1‑151 (MQ=255)                                                                                                            cGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCa                                          >  2:827680/1‑150 (MQ=255)                                                                                                                       aaaGGGTAGCAAGACTGCGGCGTGCCCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACa                              <  1:386018/151‑1 (MQ=255)                                                                                                                       aaaGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTa                                <  2:966807/149‑1 (MQ=255)                                                                                                                       aaaGGGTAACAAGCCTGCGGCGTGCCCGCACTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTa                                <  1:298133/149‑1 (MQ=255)                                                                                                                           ggTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGt                          >  2:214846/1‑151 (MQ=255)                                                                                                                             tacataGACTGCGGCGTGACCGCCATAAAAAATTTCCCGCCGCACCTCATGTCTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATc                            <  1:569855/142‑1 (MQ=255)                                                                                                                                 aaGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTg                     >  2:319071/1‑150 (MQ=255)                                                                                                                                       gcggcgTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATg              <  2:874749/151‑1 (MQ=255)                                                                                                                                         ggcgTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTa            <  2:646165/151‑1 (MQ=255)                                                                                                                                                ccGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACggg       <  2:656348/149‑1 (MQ=255)                                                                                                                                                   cgcTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAgg   <  2:442302/150‑1 (MQ=255)                                                                                                                                                    gcacAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGt  <  1:638917/147‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                   CACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGT  >  NC_000913/3177598‑3177894

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