Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,337,723	0	A	C	43.2%	61.0 / 12.9	44	G566G (GGT→GGG)	gyrA	DNA gyrase (type II topoisomerase), subunit A
Reads supporting (aligned to +/- strand):  ref base A (21/4);  new base C (0/19);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 6.43e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGA‑TTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAG  >  NC_000913/2337578‑2337865                                                                                                                                                    |                                                                                                                                               gcgcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGA‑TTTACCttt                                                                                                                                           <  2:624849/151‑1 (MQ=255) gcgcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGA‑TTTACCtt                                                                                                                                            <  2:350982/150‑1 (MQ=255)       cACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGA‑TTTACCTTTCCCGcc                                                                                                                                     >  1:309845/1‑151 (MQ=255)         cGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGA‑TTTACCTTTCCCGCCac                                                                                                                                   >  2:446374/1‑151 (MQ=255)               gCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGA‑TTTACCTTTCCCGCCACGACGCt                                                                                                                             >  1:749906/1‑151 (MQ=255)                 ttCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGA‑TTTACCTTTCCCGCCACGACGCTgc                                                                                                                           >  1:204291/1‑151 (MQ=255)                  tCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGA‑TTTACCTTTCCCGCCACGACGCTgcg                                                                                                                          >  2:671955/1‑151 (MQ=255)                       cAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGA‑TTTACCTTTCCCGCCACGACGCTGCGCTTCg                                                                                                                     >  2:126684/1‑151 (MQ=255)                        aaCTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGA‑TTTACCTTTCCCGCCACGACGCTGCGCTTCg                                                                                                                     >  2:918287/1‑150 (MQ=255)                                agcgtcaaTCGAAAAAATGCAACCACGTTTGGAGAAGTACAGAATATGTTCGTGATTATTTGCTAACAGCCGTCAGTCGATAAAGTCTTCTTCTTT‑AATACGGCCTGAAGA‑TTTACCTTTCTCACTACGACGCTGCGCTTAGTATTCAg                                                                                                              >  2:757018/8‑149 (MQ=255)                                    ttCATCTAATAGCCGCGACCACGGTTGGAAAAGCACAGAATAGGGTGGTAAGTGTTCGCCCCCAGCAGTCGGTCGATAAAGTCTTCTTTTTT‑AATACGTGCGGCAGA‑TTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGc                                                                                                         <  2:868874/150‑1 (MQ=255)                                      cATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCATCTTCTTT‑AATACGTGCGGCAGA‑TTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGc                                                                                                      >  2:577015/1‑151 (MQ=255)                                       agcGAATAGAGGCGCCCGCGGTGGGAGAAGCAAAGAATTTGGTCGTGAGTTTTCGCCACCATCAGTCGGTCAAAAAAGTCTTTTTTTTT‑AATGCGTGCGGAATT‑TTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCt                                                                                                     <  2:795817/149‑1 (MQ=255)                                           aaTGGACGCGCCCAGGGCTGGAGAAGCACAGATTTTGGTGGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAATTCTTCTTTTTT‑AATACGTGGGGAAGT‑TTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCt                                                                                                     <  1:573755/147‑1 (MQ=255)                                           aaTAGCCGCGACCCGGGCGGGAGAAGCAAAAAATAGGGTCATGGGTGTTCGCCCCCAGCAGTCGGTCGATAAAGTCTTCTTCTTA‑AATCCGTGGGGAAGT‑TTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGc                                                                                                      <  2:662781/146‑1 (MQ=255)                                                 cgcgCCCGCGGCGGAAGAACCACAAAATATGGTCGTGTGTGTCCGCCACCAGCAGTGGGTCGATAATTTTTTTTTCTTA‑AATGCGTGCGGCAGA‑TTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTaa                                                                                            <  2:826483/150‑1 (MQ=255)                                                      ccgcGGCTGAAGAAGCACAAAATATGGTCGGGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGA‑TTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAg                                                                                       <  1:201571/147‑1 (MQ=255)                                                      ccccGGCTGGAAAAGCAAAGAATATGGTCGGGAGGTTCCCCCCCCAGCAGTCGGTCAATAAAGTCTTTTTCTTT‑AATCCGTGCGGAAGA‑TTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACtt                                                                                              <  2:14170/140‑1 (MQ=255)                                                      ccACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGA‑TTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGc                                                                                      >  1:804425/1‑151 (MQ=255)                                                      ccACGGCTGGAAAAGCACAGAATATGGTCGTGAGTGTTCCCCCCCAGCGGTCGGTCGATAAAGTCTTCTTCTTT‑AATACGTGCGGAAGA‑TTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAg                                                                                       <  1:1061198/150‑1 (MQ=255)                                                           ggtGGGGAAACAAAAAATGTGGTCGGGTGGGTCCGCCACCGGCATGCGCGCGATAAAGTCTTTTTTTTT‑AACACGGGCGGAAGT‑TGTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCt                                                                                   <  2:55297/147‑1 (MQ=255)                                                                  agaccgaGAATAGGGCCGGGATTGTCCCCCACCAGCAGTGGGTC‑AAAAAATCCTCTTTTTTAAATAGGGGCGG‑AGATTTTCCCTTTCCCGCCACGACGCTGCGCTTCTTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTgag                                                                             <  2:739833/142‑1 (MQ=255)                                                                     cacaAAATATGGTCGTGATTTTTCCCCCCCAGCGGTCGTTCGATAAAGTTTTCTTTTTT‑AATCGGTGCGGCAGA‑TTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGc                                                                        <  1:335346/150‑1 (MQ=255)                                                                         gAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTT‑GATACGTGCGGCAGA‑TTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCAc                                                                   >  1:192228/1‑151 (MQ=255)                                                                                     gggtgtTCGCCCCCAGCAGTCGGTCGATAAAGTCTTCTTCTTT‑AATCCGGGCGGCAGT‑TTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCtt                                                        <  1:513167/148‑1 (MQ=255)                                                                                     gAGTGTCCGCCACCAGCAGTCGGTCGATAAATTCTTCTTCTTT‑AATCCGTGCGGCAGT‑TTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCtt                                                        <  1:904393/150‑1 (MQ=255)                                                                                             gCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGA‑TTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTg                                                >  2:655654/1‑150 (MQ=255)                                                                                             gCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGA‑TTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTg                                                >  1:1034040/1‑150 (MQ=255)                                                                                              cccccaGCAGTCGGTCGAAAAAGTCTTCTTCTTT‑AATACGTGCGGAAGA‑TTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGa                                               <  1:848072/147‑1 (MQ=255)                                                                                                    gcagTCGGTCGATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGA‑TTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGAt                                         <  1:786338/150‑1 (MQ=255)                                                                                                          ggtcgATAAATTCTTCTTTTTT‑AAGACGGGCGGAAGT‑TTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCtt                                      <  2:277692/147‑1 (MQ=255)                                                                                                          ggtcgATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGT‑TTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTcc                                    <  1:114442/149‑1 (MQ=255)                                                                                                            tcgATAAAGTCTTCTTTTTT‑AATCCTTGCGGCAGT‑TTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAgg                                 <  1:546553/150‑1 (MQ=255)                                                                                                            tcgATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGA‑TTTACCTTTCCCGCCACGACTCTGCTATTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGt                                >  2:226069/1‑151 (MQ=255)                                                                                                             cgATAAGGTCTTTTTCTTA‑AATCCGTGCGGGAGA‑TTTCCTTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAgg                                 <  1:571603/149‑1 (MQ=255)                                                                                                               aaaaaGTTTTCTTCTTT‑AAGACGTGCGGCAGA‑TTTACCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTg                              <  2:749612/148‑1 (MQ=255)                                                                                                                 aaaTTCTTCTTCTTT‑AATCCGTGCGGCAGA‑TTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAgg                                 <  1:4185/145‑1 (MQ=255)                                                                                                                   aGTCTTCTTCTTT‑AATACGTGCGGCAGA‑TTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTc                         >  2:740775/1‑151 (MQ=255)                                                                                                                      cttcttcttT‑AATACGTGCGGCAGA‑TTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTgc                      >  2:1066701/1‑151 (MQ=255)                                                                                                                         cttcttT‑AATACGTGCGGCAGA‑TTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCt                   >  1:129651/1‑151 (MQ=255)                                                                                                                             ttt‑AATACGTGCGGCAGA‑TTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTg               >  1:405542/1‑151 (MQ=255)                                                                                                                              tt‑AATACGTGCGGCAGA‑TTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTgg              >  2:107544/1‑151 (MQ=255)                                                                                                                                      gTGCGGCAGA‑TTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGAt       >  2:13510/1‑151 (MQ=255)                                                                                                                                            cAGA‑TTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAg  >  2:619036/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                               GCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGA‑TTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAG  >  NC_000913/2337578‑2337865

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