Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,337,723	0	A	C	40.8%	71.4 / 15.8	51	G566G (GGT→GGG)	gyrA	DNA gyrase (type II topoisomerase), subunit A
Reads supporting (aligned to +/- strand):  ref base A (27/2);  new base C (0/20);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 8.17e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACG  >  NC_000913/2337574‑2337869                                                                                                                                                      |                                                                                                                                                   ccGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTa                                                                                                                                                    <  2:510751/150‑1 (MQ=255)  cgcgcgcgCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAcc                                                                                                                                                  >  2:673037/1‑151 (MQ=255)           cACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGcc                                                                                                                                         >  2:420511/1‑151 (MQ=255)               aGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCacg                                                                                                                                      >  1:219920/1‑150 (MQ=255)                gTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCacgac                                                                                                                                    >  2:669916/1‑151 (MQ=255)                 tGGCTTCCGGCAACTGAGAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATCGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTACCGCCacgac                                                                                                                                    >  2:309161/1‑150 (MQ=255)                             aCTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAAAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATCCGGGCGGCAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTa                                                                                                                       <  1:659412/151‑1 (MQ=255)                              cTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTAt                                                                                                                      >  2:801342/1‑151 (MQ=255)                                 aaaaaCTTTCATCGAATAGCCGCGCCCACGGCTGAGGAAGCACAAATTATGGTCGGGAGTCTTCCCCCCCGGCAGTCGGTCGTAAAAGTCTTCTTCTTGAATCCGTGCGCCATATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCa                                                                                                                   <  2:816248/149‑1 (MQ=255)                                   aaaCTTTCATCGAATAGACGCGACCCCGGCTGGAGAAGCACAAAATAGGGTCGTGAGTGTTCGCCCCCAGCAGTCGGTCAAAAAAGTCTTCTTTTTAAATACGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAg                                                                                                                  <  1:751267/150‑1 (MQ=255)                                   aaaCTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTc                                                                                                                    >  1:213049/1‑148 (MQ=255)                                    aaCTTTCATTAAATGGACGCGCCCGCGGCTGGAGAACCACAGAAGATGTGCGGGAGTGTCCGCACCCAGCAGTCGGCCGATAAAGTCTTCTGTTTTAATAGGTGCGGCAGTTTTCCCTTACCCGCCACGACGCTGCGCTTCGTATTc                                                                                                                    <  2:491843/147‑1 (MQ=255)                                       tttCATCGAATAGACGCGACCACGGCTGGAAAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTTTTTAATACGTGGGGCAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGaaa                                                                                                               <  1:700135/149‑1 (MQ=255)                                             cGAATAGACGCGCCCACGCCGGGAGAAGCACAAAATGTGGCCGTAAGTGTTCCCCACCAGCAGTCGGTCAATAAAGTCTTCTTTTTTAATCCGTGCGGAAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTg                                                                                                        <  2:654769/150‑1 (MQ=255)                                             cGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCt                                                                                                         >  1:418697/1‑149 (MQ=255)                                              gAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGAt                                                                                                      >  1:6712/1‑151 (MQ=255)                                               aaTAGACGCGCCCGCGGCTGGAGAAGCACAAAATATGGTCGGGAGTTTTCGCCCCCAGCGGTGGGTCAATAAAGTCTTCTTCTTTAATACGTGCGGAAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATa                                                                                                     <  1:200719/151‑1 (MQ=255)                                                aTAGACGCACCCAGGGTGGGAGAAGCAAAAAATATGTTGGGGGGTGTTCCCCACAAGCAGTCGGTCGAAAAATTTTTTTTTTTTAATACGTGGGGAAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATAc                                                                                                    <  2:741155/151‑1 (MQ=255)                                                          ccACGGCTGGAGAAGCACAGAATATGGTCGTGATTGTTCGCCCCCAGCAGTCGGTCAATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACg                                                                                              <  2:252390/147‑1 (MQ=255)                                                           cACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGcc                                                                                         >  1:331740/1‑151 (MQ=255)                                                              ggCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCCCCAGCAGTCGGTCGATAAAGTCTTTTTCTTTAATACGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTg                                                                                      <  1:375493/151‑1 (MQ=255)                                                                 tGGAGAAGCAAAGAATGTTGTCGGGAGTTTCCGCCCCCAGGAGGTGGTCGATAAAGTTTTCTTCTTTAGTACGTGCGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCt                                                                                       <  2:28690/147‑1 (MQ=255)                                                                                tatGGGCGTGTGTGTTCGCCACCGGCAGCCGGTCGATAAACTTTTTTTTTTTAAAACGTGGGGCAAATTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCa                                                                        <  1:733032/147‑1 (MQ=255)                                                                                  tGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGAcc                                                                   >  2:167837/1‑150 (MQ=255)                                                                                        tGAGTGTTCCCCACCAGCAGTGGGTCAATAAAGTCTTCTTTTTTAATCCGGGCGGCAGATTCCCTTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACcac                                                                 <  1:402192/146‑1 (MQ=255)                                                                                         gAGTGTTCGCAACCAGCAGGCGGTCGATAAAGTTTTATTATTTAATACGTGCGGCAGATTTACCCTTCCCGCCACGACGCTGCGCCTCGTATTCCGAACGCGACTGGTACTTAACGTAGCCCTGGTGAGCGAGCGTCACGACCACATCtt                                                            >  2:206805/1‑150 (MQ=255)                                                                                             gtTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTgg                                                       >  2:233732/1‑151 (MQ=255)                                                                                              ttCCCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGAAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTg                                                        <  1:732364/149‑1 (MQ=255)                                                                                                 gCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTg                                                    >  2:279613/1‑150 (MQ=255)                                                                                                 ccccccaGCAGTCGGTCAAAAAAGTCTTCTTCTTTAATACGTGCGGCATTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTg                                                        <  2:844141/142‑1 (MQ=255)                                                                                                   caccaGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATc                                                 >  1:369034/1‑151 (MQ=255)                                                                                                     ccaGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATc                                                 >  2:858799/1‑149 (MQ=255)                                                                                                         cagTCGTTCGATAAAGTCTTCTTCTTAAATCCGGGCGGAAGATTTCCTTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATc                                            <  1:478812/150‑1 (MQ=255)                                                                                                           gtcggtcgATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTc                                         >  1:150287/1‑151 (MQ=255)                                                                                                           ggcgctggATAAAGTCTTTTTCTTTAATAGGGGGGCAAGTTTTTCCTTGCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCtt                                          <  1:168572/143‑1 (MQ=255)                                                                                                              gtcctaaaaaGTCTTTTTCTTTACTTGTTGCGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCa                                       <  2:279579/143‑1 (MQ=255)                                                                                                                   aTAAAGTCTTTTTTTTTAATGCGTGGGGCATTTTTGCCTTACCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTg                                  <  1:634935/150‑1 (MQ=255)                                                                                                                    tAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGAt                                >  1:789914/1‑151 (MQ=255)                                                                                                                     aaaGTCTTCTTCTTTAATCCGTGCGGCAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTg                                  <  1:311565/148‑1 (MQ=255)                                                                                                                     aaaGTCTTCTTCTTAAATACGGGCGAAAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATg                               <  1:816479/151‑1 (MQ=255)                                                                                                                      aaGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGt                              >  2:23803/1‑151 (MQ=255)                                                                                                                         tcttcttctTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTg                           >  1:218604/1‑151 (MQ=255)                                                                                                                          cttcttcttTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTgc                          >  1:112882/1‑151 (MQ=255)                                                                                                                           ttcttcttTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTgcg                         >  1:468365/1‑151 (MQ=255)                                                                                                                              ttcttAAATGCGTGCGGAAGTTTTACTTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGa                                 <  2:330037/140‑1 (MQ=255)                                                                                                                                      tACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGt              >  1:817325/1‑151 (MQ=255)                                                                                                                                        cGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTg             >  1:435961/1‑150 (MQ=255)                                                                                                                                             ggtataTTTACTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCa       <  1:359469/146‑1 (MQ=255)                                                                                                                                              gCACATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGAttt         <  1:824437/148‑1 (MQ=255)                                                                                                                                                  aTTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCCTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTacg  >  1:519182/1‑151 (MQ=255)                                                                                                                                                  aTTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTACCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTacg  >  2:342016/1‑151 (MQ=255)                                                                                                                                                   tttACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGCTTTCAGTacg  >  2:421334/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACG  >  NC_000913/2337574‑2337869

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