Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	366,526	0	T	G	48.8%	52.0 / 15.0	44	T329P (ACC→CCC)	lacI	lactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (5/17);  new base G (21/0);  total (27/17)
Fisher's exact test for biased strand distribution p-value = 7.67e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGC  >  NC_000913/366393‑366675                                                                                                                                       |                                                                                                                                                       aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCCAGGGTGGTTTTTCTTTTCACCa                                                                                                                                       <  2:278594/151‑1 (MQ=255)    gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑TGGCGCCAGGGTGGGTTTTTTTTTTACCAGt                                                                                                                                     >  1:754258/1‑150 (MQ=255)    gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCCAGGGTGGTTTTTCTTTTCACCAGTg                                                                                                                                    <  2:375424/151‑1 (MQ=255)         agCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTATT‑TGGCGCCCGGGGGGGTTTTTTTTTTACCAGGGAGACg                                                                                                                               >  2:601639/1‑151 (MQ=255)           cTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGGGGCCGGGGGGTTTTTTTTTTCAcccaggagaaggg                                                                                                                             >  1:269162/1‑140 (MQ=255)             aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCGCATT‑GGGGGGCAGGGGGGGTTTTCTTTTTACCAGGGAGAgggg                                                                                                                             >  2:100939/1‑149 (MQ=255)                 cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑TGGGGCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAAc                                                                                                                         >  1:362402/1‑149 (MQ=255)                   caTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCt                                                                                                                     <  2:27608/151‑1 (MQ=255)                    aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATT‑TGGGGCCCGGGGGGGTTTTCTTTTCACCAGTGAGACGGGCaa                                                                                                                          >  2:651529/1‑145 (MQ=255)                     ttAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATT‑TGGGGCCCGGGGGGTTTTTCTTTTTACCAGTGAGACGGGCAACAGCTg                                                                                                                    >  1:178466/1‑150 (MQ=255)                      tAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGCGTTAGGGGGCCGGGGGGGTTTTT‑TTTTCGCCAAAGAGAAGGGCAACAg                                                                                                                       >  1:735633/1‑146 (MQ=255)                           gcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATT‑TGGGGCCAGGGGGGGTTTTCTTTTTACCAGTGAGACGGGGCACAGCGG‑ATTGgcc                                                                                                            >  2:689371/1‑151 (MQ=255)                            cgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTATT‑TGGCGCCAGGGGGGGTTTTCTTTTTACCAGGGGGACGGGGAACAGCTG‑GTTGgcc                                                                                                            >  2:544822/1‑150 (MQ=255)                             gttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATT‑TGGCGCCCGGGGGGTTTTTCTTTTTACCAGTGAGACGGGCAACAGCTG‑GTTGccc                                                                                                            >  2:115482/1‑149 (MQ=255)                                         gCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCGCCTGGcc                                                                                              <  2:315827/151‑1 (MQ=255)                                           ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGGCGCC‑GGGGGGGTTTTCTCTTCTCCGCGGAGAAGGGGAAAAGCGGAATTGCCCtccccgccctgcccc                                                                                             >  2:450183/1‑137 (MQ=255)                                            cGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCGCCTGGCCCTg                                                                                           <  1:322128/151‑1 (MQ=255)                                             gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCGGATC‑GGGGGCCCGGGTGGGTTTTTTTTTTACCAGTGGGACGGGGAACAGCTG‑ATTGCCCCTCACCCCCTGGccc                                                                                             >  2:494666/1‑148 (MQ=255)                                                 tCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCGCCTGGCCCTgagag                                                                                       <  1:607010/150‑1 (MQ=255)                                                        gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑TGGCGCCAGGGGGGGTTTTTTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCCCCCGGCCCTgggag                                                                                       >  2:351741/1‑143 (MQ=255)                                                        gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgcag                                                                               <  2:215411/151‑1 (MQ=255)                                                         ggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgcagc                                                                              <  2:123545/151‑1 (MQ=255)                                                         ggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgcagc                                                                              <  2:422219/151‑1 (MQ=255)                                                          gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgcagca                                                                             <  2:775/151‑1 (MQ=255)                                                                   cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATT‑TGGCGCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGGAAAAGCTG‑ATTGCCCTTCACCGCCTGGCCCTGGGAGAGGTGCAGCAAACGGGcc                                                                     >  1:710979/1‑150 (MQ=255)                                                                    gTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATT‑GGGCGCCCGGGGGGTTTTTTTTTTTACCAGTGAGACGGGCAACAGCTG‑GATTCCCTTCACCCCCCGGCCCCGGGAGAGTTGCAGCAAACGGGCCAc                                                                   >  1:242975/1‑151 (MQ=255)                                                                    gTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑TGGCGCCAGGGGGGGTTTTCTTTTCACCAGTGAGACGGGGAACAGCTG‑ATTGCCCTTCACCCCCCGGCCCTGAGAGAGGTGCAGCAAGCGGTCCAc                                                                   >  2:204724/1‑151 (MQ=255)                                                                       ccAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATT‑GGGGGCCCGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTTACCGCCCGGCCCTGAGAGGGGTTCAGCAAGCGGGCCACGc                                                                 >  2:798178/1‑150 (MQ=255)                                                                         aGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTgg                                                              <  2:666883/151‑1 (MQ=255)                                                                          gCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGTTT‑TGGTGGCAGGGAGGGTTTTTTTTTCTCCCATTCGACCGGAAACACCTG‑AATGCCCTTGACCTCCCGGCCCTGAGAGAGCTGAAGCAAGCGGTCCACGCTgg                                                              >  1:385919/1‑150 (MQ=255)                                                                            tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGCATT‑TGGGGCCCGGGGGGGTTTTTTTTTTACCAGCGAGACGGGCAACAGCCG‑GTTGCCCTTTACCCCCCGGCCCTTGGGGAGTTTCAGCAAGCGGTCCACGCCGGttt                                                           >  1:717506/1‑151 (MQ=255)                                                                                 tAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGccc                                                       <  1:587482/150‑1 (MQ=255)                                                                                  aaTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCca                                                     <  2:218294/151‑1 (MQ=255)                                                                                      aaTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagcag                                                 <  2:813828/151‑1 (MQ=255)                                                                                        tCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATT‑TGGCGCCCGGGGGGGTTTTTTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCCTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGGCCACCCTGGGTTGCCCCAGCAgg                                                >  1:61491/1‑150 (MQ=255)                                                                                          ggCCAACGCGCGGGGAGAGGCGGTTTTCGGATT‑GGGCGCCAGGGGGGTTTTTCTTTTCACCAAGGAGAAGGGCAACAGGTG‑ATTGCCCTTCCCCGCCTCGCCCCGAGAGAGGGGTAGCAAGCGGGCCACGCCGGTTTGcccc                                                      >  2:425492/1‑142 (MQ=255)                                                                                          ggCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑TGGCGCCAGGGGGGGTTTTCTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCGCCTGGCCCTGGGAGAGTTGCAGCAAGCGGGCCCCGCTGGGTTGCCCCAGCAGGCg                                              >  2:14652/1‑150 (MQ=255)                                                                                              aaCGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAAt                                         <  2:216459/151‑1 (MQ=255)                                                                                                 gcgcgGGGAGAGGCGGTTTGCGTATT‑TGGCGCCAGGGGGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGGCCACGCTGGTTTGCCCCAGCAGGCGAAAATcc                                       >  1:117546/1‑150 (MQ=255)                                                                                                      gggAGAGGCGGTTTGCGTATT‑GGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTg                                 <  1:293303/151‑1 (MQ=255)                                                                                                        gagaGGCGGTTTGCGTATT‑GGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAt                               <  2:851779/151‑1 (MQ=255)                                                                                                                          t‑GGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg               >  1:608310/1‑149 (MQ=255)                                                                                                                                gcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa         >  1:393589/1‑150 (MQ=255)                                                                                                                                 ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAAc        >  1:267626/1‑150 (MQ=255)                                                                                                                                       tggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGc  >  1:747464/1‑150 (MQ=255)                                                                                                                                       |                                                                                                                                                       AATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATT‑GGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTG‑ATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGC  >  NC_000913/366393‑366675

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