Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,5260TG46.7% 61.3 / 13.9 45T329P (ACC→CCC) lacIlactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (5/19);  new base G (21/0);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 3.06e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCT  >  NC_000913/366382‑366676
                                                                                                                                                |                                                                                                                                                       
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)
                                                                                                                                      ggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat            >  1:565139/1‑151 (MQ=255)
                                                                                                                                         gcCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGAtata           <  1:304594/149‑1 (MQ=255)
                                                                                                                                           cAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAt       <  2:719490/151‑1 (MQ=255)
                                                                                                                                                tggtTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCt  >  1:439190/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                       
CTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCT  >  NC_000913/366382‑366676

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