Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	366,515	0	G	T	40.9%	49.4 / 14.8	45	P332P (CCC→CCA)	lacI	lactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (9/17);  new base T (18/0);  total (28/17)
Fisher's exact test for biased strand distribution p-value = 4.94e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

ATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATA  >  NC_000913/366368‑366665                                                                                                                                                    |                                                                                                                                                        aTAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCGCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGAGGGGAGAGGCGGTTTGCGTATTGGgc                                                                                                                                                      <  2:729393/151‑1 (MQ=255) aTAAAGGGTAAAGCCGGGGGGGCATAAGGAGTGAGCTAACTCACAGTAAAAGAGAGGGGCCGACGGGCCGCTAGAAAGCAGGGAAAACTGTCGCGCCAGATGCATTAATGAAGCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTGTTtgg                                                                                                                                                       >  2:731939/1‑150 (MQ=255)    aaGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTTTTGGGGGcc                                                                                                                                                   >  1:318525/1‑151 (MQ=255)      gtgtAAAGCCTGGGGGGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTACTGCATCGGCCAACGCGCGGGGAGAGGCGGGCTGCGGATTAGgcgaccg                                                                                                                                                 >  2:711978/1‑147 (MQ=255)         tAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTTGGGGCCCggg                                                                                                                                               >  1:57909/1‑150 (MQ=255)               cTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGGGGgttt                                                                                                                                        >  1:56571/1‑151 (MQ=255)                tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCCGGGGGGgtttt                                                                                                                                       >  1:506631/1‑151 (MQ=255)                      gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttcttttc                                                                                                                                 <  1:529045/151‑1 (MQ=255)                         tAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGTGTATTTGGGGCCACGGGGGTTTTTTTTTTTAcc                                                                                                                              >  1:443558/1‑151 (MQ=255)                         tAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGGATTGGGGGGCAGGGGGGGttttttttttaa                                                                                                                               >  1:400239/1‑147 (MQ=255)                         tAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGCCCGGGGGGGTTTTCTTTTTAcc                                                                                                                              >  1:491660/1‑151 (MQ=255)                             gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGGGCCCGGGGGGGTTTTTTTTTTAcc                                                                                                                              >  1:766595/1‑147 (MQ=255)                                 gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGTGGATTTGGGGCCAGGGGGGGTTTTTTTTTTACCAGTg                                                                                                                          >  1:267949/1‑147 (MQ=255)                                     tAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCAGGGGGGGTTTTTTTTTTACCAGTGAG‑ACggg                                                                                                                  >  1:189177/1‑150 (MQ=255)                                      aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACggg                                                                                                                  <  1:208855/149‑1 (MQ=255)                                       aCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGTGGGTTTGGGGCCCGGGGGGGttttttttttccccgggag                                                                                                                        >  2:443237/1‑133 (MQ=255)                                          cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTTGGCGCCCGGGGGGGTTTTTTTTTTTCCCGTGAG‑ACggggca                                                                                                               >  1:334474/1‑145 (MQ=255)                                          cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGATTTCGGATTGGGGGCCCGGGTGGTTTTTTTTTTTCCCAGCTGG‑ACGGGCaa                                                                                                               >  2:678854/1‑148 (MQ=255)                                            caTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTGTTTGGGGCCAGGGGGGGTTTTCTTTTCACCAGCGAG‑CAGGGGAACa                                                                                                             >  2:255599/1‑148 (MQ=255)                                                aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCAGGGGGGTTTTTTTTTTTACCCGTGAG‑ACGGGGAACAGCTGAtt                                                                                                      >  1:724306/1‑151 (MQ=255)                                                      gttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGTATTTGGGGCCCGGGGGGGTTTTTTTTTTACCCGTGGG‑AAGGGCAACAACCGATTGCCCtt                                                                                                >  1:653216/1‑151 (MQ=255)                                                        tgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTc                                                                                               <  2:734822/150‑1 (MQ=255)                                                         gcgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTTATTGCCCTTCa                                                                                              <  2:592428/150‑1 (MQ=255)                                                           gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGTGAATTGGGGGGCCGGGGGGGTTTTCTTTTTACCCAGGAGAAGGGGCAACAG‑TGATTGCCCTTCTcc                                                                                            >  2:83930/1‑150 (MQ=255)                                                                cTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTATTTGGCGCCAGGGGGGGTTTTTTTTTCACCAGTGAG‑ACGGGCAAAAGCTGATTGCCCTTCACCGcc                                                                                         >  1:455231/1‑148 (MQ=255)                                                                  gCCCGCTTTCCAGTCGGGAAACCTGTCTTGCCAGCTGCATTAATGAAGCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCGCCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTgg                                                                                      <  2:443299/149‑1 (MQ=255)                                                                        tttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCAGGGTGGGTTTTTTTTTTACCAGTGAG‑AAGGGCAACAGCTGATTGCCCCTCACCCCCCGGccccgcg                                                                               >  1:695792/1‑146 (MQ=255)                                                                             aGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGTGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGtt                                                                         <  1:217685/151‑1 (MQ=255)                                                                                cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGGGGTTTTTTTTTTTACCAGTGGG‑AAGGGGCACAACTGATTGCCCTTCACCCCCCGGGCCTgggaggggt                                                                         >  1:213024/1‑144 (MQ=255)                                                                                  ggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGCTTTGGGGCCCGGGGGGGTTTTCCTTTCACCCGGGAG‑ACGGGCAACAGCTGATTGCCCCTCACCGCCTGGCCCCGagagaggt                                                                         >  2:748793/1‑144 (MQ=255)                                                                                     aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAg                                                                 <  1:165748/151‑1 (MQ=255)                                                                                          gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTc                                                            <  2:489998/151‑1 (MQ=255)                                                                                           tCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGGGGGTTTTTTTTTTACCCGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTcc                                                           >  1:443407/1‑151 (MQ=255)                                                                                                     tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGTGTATTGGGGGCCAGGGGGGGGTTTCTTTTTACCACGGAG‑AAGGGCAAGCGCTGCTTGCCCCTCCCCGCCCGGCCCCGGGGGGGGTGGAGGAAACGGTccgcact                                                      >  1:220347/1‑141 (MQ=255)                                                                                                         ttAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGGGGTTTTTCTTTTTACCCGTGAG‑AAGGGCAACAGCTGATTTCCCTTCACCGCCCGGCCCTGAGAGAGTTGCAGCAAACGGGCCACGCTGGGTTGccc                                             >  1:454989/1‑151 (MQ=255)                                                                                                           aaTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCca                                           <  1:127380/151‑1 (MQ=255)                                                                                                                    gCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaa                                  <  2:113913/151‑1 (MQ=255)                                                                                                                    gCCAACGCGCGGGGAGAGGCGGCTTGCGGCTTGGGGGCCAGGGGGGGTTTTTTTTTCTCCAGGGAG‑ACGGGCAACAGCTGGTTGGCCCTCTCCGCCTGGCCCTGAGAGAGTTGCCGCAAGCGGGCCCCGCCGGGTTTCCCCCGCCGGCg                                    >  2:136496/1‑149 (MQ=255)                                                                                                                               gggAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTg                       <  1:184699/151‑1 (MQ=255)                                                                                                                                     ggCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtg                 >  2:655541/1‑151 (MQ=255)                                                                                                                                      gCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtgg                <  2:340068/151‑1 (MQ=255)                                                                                                                                       cGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATggtggt               >  2:554954/1‑151 (MQ=255)                                                                                                                                            tGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAc           <  2:481894/150‑1 (MQ=255)                                                                                                                                            tGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAc           <  1:651555/150‑1 (MQ=255)                                                                                                                                              cGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc        <  1:316244/151‑1 (MQ=255)                                                                                                                                               gTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc        <  2:664595/150‑1 (MQ=255)                                                                                                                                                    gggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTGCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGata  <  2:178620/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                        ATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATA  >  NC_000913/366368‑366665

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