Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,245,234	0	G	T	41.0%	50.7 / 15.8	39	L25F (TTG→TTT)	ycgY	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (3/20);  new base T (16/0);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 2.57e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTG  >  NC_000913/1245097‑1245381                                                                                                                                          |                                                                                                                                                    tAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTAttt                                                                                                                                        >  2:572942/1‑150 (MQ=255)   aaGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGGTGCGCTTAtt                                                                                                                                         >  2:347810/1‑148 (MQ=255)      tCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGTTGCGCTTATTTaaaa                                                                                                                                    >  1:23460/1‑148 (MQ=255)      tCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACa                                                                                                                                    <  1:3033/150‑1 (MQ=255)       cGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACAtt                                                                                                                                  <  1:575826/151‑1 (MQ=255)                gTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGa                                                                                                                         <  2:329980/151‑1 (MQ=255)                    caacaacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTTTTGAACATTTTGAATTGATCa                                                                                                                      >  1:67635/1‑150 (MQ=255)                     aacaacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTATTGAAAATTTAGAATTTATCaa                                                                                                                     >  2:119638/1‑150 (MQ=255)                           aGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGTTGCGCTTTTTTAAAAATTTGAATTgtttaaatctttt                                                                                                              >  1:444792/1‑139 (MQ=255)                           aGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCt                                                                                                              <  2:474615/151‑1 (MQ=255)                            gAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTATTGTACATTTAGAAATGATTAATTCtt                                                                                                                >  2:340038/1‑148 (MQ=255)                                ccTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGGGCTTATTGAACATTTTGAAATGATTAATCCCTCTgg                                                                                                            >  2:483661/1‑148 (MQ=255)                                  ttttttAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTATTTAACATTTAGAATTGATCAATCCTTCt                                                                                                              >  1:340482/1‑144 (MQ=255)                                    ttttAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGTGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGAtt                                                                                                      <  1:379225/150‑1 (MQ=255)                                       tAATTATTGTTTCGTTGTTTTTGGCGTGATTCAGAAAATAAATATTTTACACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCCTATTTAACATTTAGAATTGATCAATCCTTTTGGCCGAtt                                                                                                      >  1:525248/1‑147 (MQ=255)                                              tGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAAt                                                                                            <  2:434117/150‑1 (MQ=255)                                              tGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAAACCTTCTTGCAGATTAAGCCaaaa                                                                                             >  1:306854/1‑149 (MQ=255)                                                   cGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGt                                                                                      <  2:316809/151‑1 (MQ=255)                                                      tgttTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTTAACATTTAGAATTGATTAATACTTCTTGCAGATTAAGACAAAATACCCGtt                                                                                     >  1:553633/1‑149 (MQ=255)                                                           ttGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCCTATTTAACATTTTGAATTTATCCATCCTTCTTGGCGATTAAGCCAAAATAACCCTTCGGCaa                                                                               >  1:71379/1‑150 (MQ=255)                                                                        gAAAATAAATATTTTAGTCAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGtt                                                                     <  1:440380/147‑1 (MQ=255)                                                                            ataaataTTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTTTTTAAAATTTAGAATTGATCAATCCTTTTGGCAGATTAAGCCAAAAAACCCGTTTGGCAATGCTTCAGTTAcgcg                                                                >  2:247464/1‑148 (MQ=255)                                                                                  atTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGaa                                                        <  2:409634/150‑1 (MQ=255)                                                                                   ttttAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTATTGAAAATTTAAAATTTATCAATTCTTCTGGCAGATTAAGGCAAAATACCCGTTTGGGAATGCGTCAGTTTCGCGGAGAGGaa                                                        >  1:263717/1‑149 (MQ=255)                                                                                         aCAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAACCAAAATAACCGTTCGGCAATGCTTCAGGTACGCGAAGAgggaagg                                                     >  2:466422/1‑141 (MQ=255)                                                                                           aGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTATTTAAAATTTTGAAATTATCAATACTTCTTGGAGAATAAACCAAAAAAACCGTTTGGCAATGCCTCAGTTACGCGAAGAGGGAAGGTATGaa                                               >  1:529938/1‑150 (MQ=255)                                                                                             tAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGATTTTGCTGCGCTTATTTAAAATTTTGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCCTTCGGCAATGCTTCAGTTACACGAAGAGGGATGGTCTGAAtt                                             >  2:337066/1‑150 (MQ=255)                                                                                               gaaaGTAAGGGGTTAAAAATGGGGCGGATAGTGAGGAGGTTTTTGGTGGGGTTTTTGAAAAATTAGAATTTATTAATTATTTTTGTAGATTAAGAAAAAAAAAACGTTTGGGAATAGTTTAGTTTTGGGGAGAGGGATGGttttaaatttt                                          >  1:46963/4‑141 (MQ=255)                                                                                               gCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCtt                                          <  2:197674/151‑1 (MQ=255)                                                                                                   gTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGTTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCATATTAATCCAAAATACCCGTTCCGCAATGCTTCCCTTCCTCCACGACGAATTGTCTGCATTCtttttt                                      <  1:507902/151‑1 (MQ=255)                                                                                                     cAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGaaggaggggttt                                     <  2:498573/150‑11 (MQ=255)                                                                                                            tAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCt                              <  1:347747/150‑1 (MQ=255)                                                                                                                aaTGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAAtt                         <  1:30825/151‑1 (MQ=255)                                                                                                                        ggATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAg                 <  2:420349/151‑1 (MQ=255)                                                                                                                        ggATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAg                 <  1:247423/151‑1 (MQ=255)                                                                                                                         gATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAgt                <  2:426860/151‑1 (MQ=255)                                                                                                                          aTAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAgtg               <  1:325726/151‑1 (MQ=255)                                                                                                                                  gCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGAtt       <  1:490793/151‑1 (MQ=255)                                                                                                                                        ttGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTAtgtg  <  2:451885/150‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                    TAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTG  >  NC_000913/1245097‑1245381

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