Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,950,201	0	T	G	46.8%	65.0 / 13.0	47	G320G (GGT→GGG)	argA	amino acid N‑acetyltransferase and inactive acetylglutamate kinase
Reads supporting (aligned to +/- strand):  ref base T (11/14);  new base G (22/0);  total (33/14)
Fisher's exact test for biased strand distribution p-value = 1.40e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGC  >  NC_000913/2950056‑2950341                                                                                                                                                   |                                                                                                                                             aaGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTGTTCTGGGG‑TTGATTAGCCCCCCCGAGCAGCAAGGGATTCt                                                                                                                                         >  2:119440/1‑151 (MQ=255) aaGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGGATTTTGGGG‑GTGGTTTGCCCACCGGGGCAGCAAGGGAtttt                                                                                                                                         >  2:733377/1‑149 (MQ=255)    aTGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTTTGGGG‑TTGATTTGCCCACTGGAGCAGCAAGGGATTCTg                                                                                                                                        >  2:84967/1‑149 (MQ=255)    aTGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGGG‑TTTATTTGCCCACTGGAGCAGCAAGGGATTCTgg                                                                                                                                       >  1:490262/1‑150 (MQ=255)      ggCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTa                                                                                                                                     <  1:339931/150‑1 (MQ=255)          cgcTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTTTTCTGGAG‑TTTTGTTGGCCCCCGGAGCAACAAGGGGTTCTGGTGCGCCg                                                                                                                                >  1:412096/1‑151 (MQ=255)                tGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGGATTTTGGTACGCCGTtctc                                                                                                                           >  1:589339/1‑150 (MQ=255)                 gCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTcgc                                                                                                                         <  2:76250/151‑1 (MQ=255)                   aaGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTTTGGAG‑GTGATTCGCCCACTGGAGCAGCAAGGGATTCTGGTACGCCGTTCTcg                                                                                                                          >  2:250422/1‑148 (MQ=255)                        ttgttCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGGATTCTGGTACGCCGTTCTCGCGagcag                                                                                                                   >  2:608182/1‑150 (MQ=255)                          gttCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTg                                                                                                                <  1:289109/151‑1 (MQ=255)                          gttCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAAGGATTCTGGGACGCCGTTCTCGCGagc                                                                                                                     >  1:384127/1‑146 (MQ=255)                           ttCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTTTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGGATTCTGGGACGCCGTTCTCGCGAGCAGCTgg                                                                                                               >  1:414939/1‑151 (MQ=255)                              tcACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGGG‑TTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGGACGCCGTTCTCGCGAGCAGCTgg                                                                                                               >  2:126884/1‑148 (MQ=255)                              tcACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGGATTCTGGTACGCCGTTCTCGCGAGCAGCTgg                                                                                                               >  2:530966/1‑148 (MQ=255)                                 cgcgACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATgg                                                                                                         <  2:483413/151‑1 (MQ=255)                                  gcgACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑GTGATTCGCCCACCGGGGCAGCAAGGGATTTTGGTACGCCCTTCTCGCGAGCAGCTGGAGATgg                                                                                                         >  2:636793/1‑150 (MQ=255)                                    gACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGGG‑TTGATTTGCCCACTGGAGCAGCAAGGGATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATgg                                                                                                         >  1:572719/1‑148 (MQ=255)                                         tATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑GTGATTCGCCCCCTGGAGCAGCAAGGGATTCTGGTACGCCGGTCTCGCGCGCAGCTGGAGATGGaaa                                                                                                      >  2:623593/1‑146 (MQ=255)                                             ggTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGGATTCTGGTACGCCGTTCTCGCGAGCAGCTGGGGATGGAAATCGACaaa                                                                                              >  1:176811/1‑150 (MQ=255)                                                aCGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGGATTCTGGTACGCCGTTCTCGCGAGCAGCTGGGGATGGAAATCGACAAAtt                                                                                            >  1:132729/1‑149 (MQ=255)                                                  gCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCAcc                                                                                        <  1:240358/151‑1 (MQ=255)                                                    aGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCat                                                                                      <  1:414524/151‑1 (MQ=255)                                                    aGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCat                                                                                      <  2:264812/151‑1 (MQ=255)                                                        tgtgATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTTTGGGG‑TTGATTTGCCCACTGGAGCAGCAAGGGATTCTGGGACGCCGTTCTCGCGGGCCGCTGGAGATGGGAATCGACAAATTCACCattatt                                                                                  >  1:713673/1‑151 (MQ=255)                                                           gATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCa                                                                                <  2:5144/150‑1 (MQ=255)                                                                 aaGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGAt                                                                         <  1:563706/151‑1 (MQ=255)                                                                       cGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGGATTCTGGTACGCCGTTCTCGCGAGCAGCTGGGGATGGGAATCGACAAATTCACCATTATTCAGCGCGATAacac                                                                    >  1:50354/1‑150 (MQ=255)                                                                        gAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGa                                                                  <  2:494526/151‑1 (MQ=255)                                                                         aGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGAc                                                                 <  1:258919/151‑1 (MQ=255)                                                                         aGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGGATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGGAATCGACAAATTCACCATTATTCAGCGCGATAACACGAc                                                                 >  2:466459/1‑151 (MQ=255)                                                                          gCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGCG‑GTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGAc                                                                 >  1:586168/1‑150 (MQ=255)                                                                            aGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑GTTATTTGCCCACTGGAGCAGCAAGGGATTCTGGGACGCCGTTCTCGCCAGCAGCTGGAGGTGGAAATCGACAAATTCACCATTATTTAGCGCGATAACACGACTa                                                               >  2:521189/1‑150 (MQ=255)                                                                                      cgcAACAATCAACGATATTGGCGGTATTCTGGGG‑TTGATTTGCCCCCTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGTGCGATAACACGACTATTGCCTGcgcc                                                    >  2:290042/1‑151 (MQ=255)                                                                                         aacaaTCAACGATATTGGCGGTATTCTGGGG‑TTTATTCGCCCACTGGAGCAGCAAGGGATTCTGGGACGCCGTTCTCGCGAGCAGCTGGGGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACCACTATTGCCTGcgccgc                                                  >  1:651462/1‑150 (MQ=255)                                                                                            aaTCAACGATATTGGCGGTCTTCTGGAGATTGTTTC‑CCCACTAGAACAACAAAGGAGTTTTGGAAACCCTTCTCGCGAGGAGCTGGAGATGGAAATCGACAAATTTACCAATATTTAGCGCGAAAAAACGACTAGTGACCGCGCCGCGctc                                              >  2:358679/1‑151 (MQ=255)                                                                                                  cGATATTGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCg                                        <  1:254784/151‑1 (MQ=255)                                                                                                       ttGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTccc                                   >  2:674470/1‑151 (MQ=255)                                                                                                                 tCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGgaagaga                           >  1:210769/1‑149 (MQ=255)                                                                                                                     gAG‑TTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTg                     >  1:52339/1‑151 (MQ=255)                                                                                                                               cGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATg             <  2:697403/150‑1 (MQ=255)                                                                                                                                gCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGc           >  1:410402/1‑151 (MQ=255)                                                                                                                                      tGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCtgtgt     >  2:707187/1‑151 (MQ=255)                                                                                                                                      tGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCtgtgt     >  1:473284/1‑151 (MQ=255)                                                                                                                                      tGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCtgtgt     >  1:632773/1‑151 (MQ=255)                                                                                                                                        gAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCtgtgt     <  2:281918/149‑1 (MQ=255)                                                                                                                                           cagcaAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGc  >  2:6831/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                             AAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAG‑TTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTACGCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGC  >  NC_000913/2950056‑2950341

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: ‑