Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,133,302	0	T	G	45.5%	45.8 / 13.9	44	H219P (CAC→CCC)	yghS	putative ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (7/17);  new base G (20/0);  total (27/17)
Fisher's exact test for biased strand distribution p-value = 6.75e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGT‑CAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGC  >  NC_000913/3133152‑3133451                                                                                                                                                       |                                                                                                                                                       gAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgt                                                                                                                                                        <  1:256818/151‑1 (MQ=255)     gTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTc                                                                                                                                                     <  1:259230/150‑1 (MQ=255)              cTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTGGCCCATc                                                                                                                                              >  1:54681/1‑148 (MQ=255)                     gcgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAgg                                                                                                                                    <  1:461408/151‑1 (MQ=255)                     gcgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTTaa                                                                                                                                      >  2:509301/1‑149 (MQ=255)                          tGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTTGCCCCTCCAGGTCAAAGGTCtt                                                                                                                               >  1:112561/1‑151 (MQ=255)                              cTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGTGCCCATCCAGTTCAAGGATCttt                                                                                                                              >  1:156765/1‑148 (MQ=255)                                  tAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCa                                                                                                                        <  1:302951/150‑1 (MQ=255)                                  tAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCATCCAATTTGAAGATCTTTgcgc                                                                                                                          >  2:101482/1‑148 (MQ=255)                                       gtttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCATCCAGTTTCAGGATCTTTTCGCCCtt                                                                                                                      >  1:283446/1‑147 (MQ=255)                                          tgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCATCCAGTTCCAGGATCTTTTCGCCCtttgaatt                                                                                                                >  2:130493/1‑144 (MQ=255)                                            agtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCCCCAGTTCCAGGGTCTTTTCGCCCTTTGATGt                                                                                                                >  1:577560/1‑148 (MQ=255)                                                 gTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCATCCCCGTTAAGGGTCTTTGTGCCCTTTAATGT‑CAGtt                                                                                                          >  2:225313/1‑148 (MQ=255)                                                  tAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGTGCCCATTCAGTTTAAGGATCTTTGGGGCATTTGATGT‑GAGTTgggg                                                                                                      >  1:156126/1‑151 (MQ=255)                                                        aTCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATACAGTTCAAAGATCTTTGCGCCCTTTAAATG‑CAGTTGCGGGGTaa                                                                                                 >  2:760994/1‑150 (MQ=255)                                                          ccATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGT‑CAGTTGCGGCGTAACgg                                                                                              <  1:243506/151‑1 (MQ=255)                                                            aTTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAgg                                                                                                                                    <  1:665019/112‑1 (MQ=255)                                                            aTTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGGCGCCCATCCAGTTCAAGGATCTTTGCGCCCTTTTATTT‑CCGGTGCGGCGTTACgggg                                                                                            >  1:356875/1‑149 (MQ=255)                                                            aTTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAgg                                                                                                                                    >  2:665077/1‑112 (MQ=255)                                                                 ttAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCCTCCAACTCAAGGAAACTTGCGCCATTTTATTTGCAG‑TGCGGCGTaaacgggggttt                                                                                       >  2:548377/1‑141 (MQ=255)                                                                          ccccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGT‑CAGTTGCGGCGTAACGGCGATTTTTTCCTGTa                                                                               <  1:717291/150‑1 (MQ=255)                                                                                  cATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGT‑CAGTTGCGGCGTa                                                                                                  <  1:658430/123‑1 (MQ=255)                                                                                  cATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCCCCAGGTCTAGGATCTTTTTGCCCTTGAATGT‑CAGTTGCGGGGGGACGGGGAGTTTTTCCTGTAAAGCTGcc                                                                       >  1:514834/1‑150 (MQ=255)                                                                                   aTTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGT‑CAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAg                                                                     <  1:466110/151‑1 (MQ=255)                                                                                     tttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGT‑CAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGt                                                                    <  2:1331/150‑1 (MQ=255)                                                                                         gTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGGGGCCATCCAGTTCAAGGGACTTTTGGCCATTTAATGG‑CAGGTGGGGGGGGAAGGGGGGTTTTTTCTGTAAAGCTgtccggt                                                                   >  1:269070/1‑141 (MQ=255)                                                                                          tGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCACCAGTTTAAGGATCTTTTCGCCCTTGGAAGT‑CAGTTGGGGGGGAAAGGGGGATTTTTTCTTTAAACCTGCCCGGtggtgg                                                              >  2:294066/1‑151 (MQ=255)                                                                                           gATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTTAAAGATCTTTGCGCCCTTGAATGT‑CAGGTGCGGCGGAAACGCGGTTTTTTTCTGTTAAACCGCCCGTtggtgg                                                              >  1:766886/1‑150 (MQ=255)                                                                                              cAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAAGATCTTTTCGCCATTTAATGT‑CAGTTGCGGGGTAAAGGCGGTTTTTTTCTGTTAAGCTGCCCGtttggggggcg                                                          >  2:511682/1‑142 (MQ=255)                                                                                               aGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGGCGCCCATCCCGGTCAAAGATATTTTCGCCCTTGTATGT‑CAGTTGCGGCGGAACGGCGGTTTTTTTCTGTAAAGCTGCCAGGtggtggggcgg                                                         >  2:303649/1‑146 (MQ=255)                                                                                                    aGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAGGATCTTTTCGCCCTTGAATGT‑TAGTTGCGGCGGAAAGGCGGTTTTTTTCTTTTACGCCGCCCGTtggggggccgggtt                                                      >  1:370148/1‑141 (MQ=255)                                                                                                      cgcgGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCCGTTTAAGGATCTTTTCGCCATTGGATGG‑CAGTTGCGGCGGAACGGCGATTTTTTTCTGTAAAGCTGCCCGTTGGGGGGGAGGttt                                                      >  2:317751/1‑147 (MQ=255)                                                                                                               aaTTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCCTCCAGTTCAAGGATCTTTTCGCCATTGAATGT‑CAGTTGCGGCGGTACGGCGATTTTTTTCTGTAACGCTGCCAGTTGGTGGGCAGGTTTTCGCGCaa                                              >  2:205077/1‑146 (MQ=255)                                                                                                                   gcgcGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGT‑CAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGttt                                     <  1:225280/151‑1 (MQ=255)                                                                                                                            gACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCCGTTCAAGGATCTTTTCGCCATTTAAATT‑CAGTTGCGGCGGAACGGCGGTTTTTTTCTGTAACGCTGCCCGTTGGGGGGCAGGGTTACGCGCAAACGCGGGTTGTtcatca                             >  2:257079/1‑150 (MQ=255)                                                                                                                             aCGCTTGCAGAATTTCATCGGCGGGGGGTCGGCCATTCCGTTCAAGGAGCTTTTCCCCATTGTAAGT‑CTGTTGCGGCGGTACGGGGGATTTTTTCTTTAACGCTTCCAGGTGGGGGGGGGGGTTTTGGGCCAACGCTGGTTGTTCATTaa                            >  1:656870/1‑150 (MQ=255)                                                                                                                              cGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGT‑CAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATg                          <  1:205881/151‑1 (MQ=255)                                                                                                                                   gCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGT‑CAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAg                     <  1:596639/151‑1 (MQ=255)                                                                                                                                    cAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGT‑CAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGa                    <  1:717381/151‑1 (MQ=255)                                                                                                                                    cAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGGTCCAGGATCTTTTGGTCATTGAATGT‑CAGTTGCGGCGGAAAGGGGGTTTTTTTCTGTAAAAGTGCCCGCTGGGGGGCCGGGTTTCGCGCCAACGCGGGTTGTTCATTAATGTCaa                      >  1:408748/1‑149 (MQ=255)                                                                                                                                      gAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGT‑CAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACg                  <  1:511630/151‑1 (MQ=255)                                                                                                                                         tttCATCGGCGGGGTGTCGCCCATCCAGTTCAAAGATCTTTTTGCCATTTAATGT‑CAGTTGCGGCGTAAAGGCGCTTTTTTTCTGTTACGCTGCCAGTTGGGGGGCAGGGTTTCGCGCAAACGCGGGTTGTTTATTAATGCCCAGACgga                >  1:150995/1‑148 (MQ=255)                                                                                                                                               cggcgtGGTGTCGCGCATACAGTTCAAGGATCTTTGCGCCATTGAATGT‑CAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATcaaca          <  1:514523/150‑1 (MQ=255)                                                                                                                                               cggcggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGT‑CAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATcaacaa         <  2:151028/151‑1 (MQ=255)                                                                                                                                                     ggTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGT‑CAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACggg   <  1:747036/151‑1 (MQ=255)                                                                                                                                                       tgtCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGT‑CAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGc  >  1:513432/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                       GAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGT‑CAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGC  >  NC_000913/3133152‑3133451

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