Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,3510AT72.5% 363.6 / 152.8 178intergenic (‑46/+77)lacZ/lacIbeta‑galactosidase/DNA‑binding transcriptional repressor LacI
Reads supporting (aligned to +/- strand):  ref base A (22/27);  new base T (69/60);  total (91/87)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 191 total aligned reads displayed.

TCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTG  >  NC_000913/366302‑366400
                                                 |                                                 
 taTAGCTGTTTTCTGTGAGAAATTGTTATCCGCCCACAATTCCACacaa                                                   >  2:2424357/2‑49 (MQ=38)
 cATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaa                                                   >  1:2876219/1‑49 (MQ=255)
  aTAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaaca                                                 >  1:553861/1‑50 (MQ=255)
   tAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaaca                                                 >  1:3797315/1‑49 (MQ=255)
       tGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACATTATAc                                              >  1:4535229/1‑48 (MQ=255)
         tttCCTGTGTGAAATTGTTATCCGCTCACAATTCCACGCAACATACGAGc                                          <  1:432681/50‑1 (MQ=255)
         tttCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACATTATACGAGc                                          <  2:2590476/50‑1 (MQ=255)
         tttCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGc                                          <  2:115037/50‑1 (MQ=255)
          ttCCTGTGTGTAATTGTTATCCGCTCACAATTCCACACATTATACGCGcc                                         >  1:625691/1‑50 (MQ=255)
          ttCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACATTATACGAGcc                                         >  1:819888/1‑50 (MQ=255)
          ttCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACATTATACGAGcc                                         >  1:974003/1‑50 (MQ=255)
          ttCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACATTATACGAGcc                                         >  1:378494/1‑50 (MQ=255)
           tCCTGTGTGAAATTGTTATCCGCTTACAATTCCACACAACATACGAGcc                                         <  2:4665172/49‑1 (MQ=255)
           tCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACATTATACGAGcc                                         <  1:3341979/49‑1 (MQ=255)
           tCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACATTATACGAGCCg                                        >  1:1134350/1‑50 (MQ=255)
            ccTGTGTGAAATTGTTATCCGCTCACAATTCCACACATTATACGAGCCg                                        <  1:1491729/49‑1 (MQ=255)
             cTGTGTGAAATTGTTATCCGCTCACAATTCCACACCTTATACGAGCCGGA                                      >  1:2212312‑M1/1‑50 (MQ=255)
             cTGTGTGAAATTGTTATCCGCTCACAATTCCACACATTATACGAGCCGGA                                      >  1:1073061‑M1/1‑50 (MQ=255)
              tgtgtgAAATTGTTATCCGCTCACAATTCCACACATTATACGAGCCGGAt                                     >  1:1637149‑M1/1‑49 (MQ=255)
                tgtgAAATTGTTATCCGCTCACAATTCCACACATTATACGAGCCGGAtga                                   >  2:3803426‑M1/1‑47 (MQ=255)
                tgtgAAATTGTTATCCGCTCACAATTCCACACATTATACGAGCCGGAt                                     <  2:3550813‑M1/48‑2 (MQ=255)
                tgtgAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGc                                   <  2:3974656/50‑1 (MQ=255)
                 gtgAAATTGTTATCCGCTCACAATTCCACACATTATACGAGCCGGAtgat                                  >  2:52834‑M1/1‑46 (MQ=255)
                  tgAAATTGTTATCCGCTCACAATTCCACACATTATACGAGCCGGAtgat                                  <  2:4453613‑M1/49‑5 (MQ=255)
                   gAAATTGTTATCCGCTCACAATTCCACACATTATACGAGCCGGAtgatt                                 >  1:1867069‑M1/1‑44 (MQ=255)
                    aaaTTGTTATCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaa                               >  2:167381‑M1/1‑43 (MQ=255)
                    aaaTTGTTATCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaa                               <  2:2401860‑M1/50‑8 (MQ=255)
                    aaaTTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATaa                               <  1:4957213/50‑1 (MQ=255)
                    aaaTTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATaa                               <  1:2057560/50‑1 (MQ=255)
                     aaTTGTTATCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaat                              <  2:4803457‑M1/50‑9 (MQ=255)
                      aTTGTTATCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaatt                             <  2:1558310‑M1/50‑10 (MQ=255)
                      aTTGTTATCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaatt                             <  1:1458626‑M1/50‑10 (MQ=255)
                      aTTGTTATCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaatt                             >  1:3494314‑M1/1‑41 (MQ=255)
                      aTTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAg                             <  1:3404796/50‑1 (MQ=255)
                       ttGTTATCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaattg                            <  2:1651771‑M1/50‑11 (MQ=255)
                       ttGTTATCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaattg                            <  1:2478312‑M1/50‑11 (MQ=255)
                         gTTATCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaattgt                           <  2:3830564‑M1/49‑12 (MQ=255)
                         gTTATCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaattgt                           <  2:646381‑M1/49‑12 (MQ=255)
                          ttATCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaattgtca                         <  1:3762394‑M1/50‑14 (MQ=255)
                          ttATCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaattgtca                         <  1:3261398‑M1/50‑14 (MQ=255)
                          ttATCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaattgtca                         <  2:5383268‑M1/50‑14 (MQ=255)
                          ttATCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaattgtc                          >  1:5095634‑M1/1‑37 (MQ=255)
                          ttATCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaattgtc                          >  1:2607463‑M1/1‑37 (MQ=255)
                          ttATCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaattgtc                          >  2:898521‑M1/1‑37 (MQ=255)
                          ttATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgtgt                          <  1:3736392/49‑1 (MQ=255)
                           tATCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaattgtcaa                        >  1:1787913‑M1/1‑36 (MQ=255)
                           tATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTa                         >  1:2702939/1‑49 (MQ=255)
                            ctcCGCTCACAATTCCACACATTCTACGAGCCGGAtgattaattgtcaa                        <  2:2895154‑M1/48‑15 (MQ=255)
                            aTCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaattgtcaa                        >  2:4561917‑M1/1‑35 (MQ=255)
                            aTCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaattgtcaa                        >  2:2670356‑M1/1‑35 (MQ=255)
                            aTCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaattgtcaa                        >  2:3257270‑M1/1‑35 (MQ=255)
                            aTCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaattgtcaa                        >  1:4614740‑M1/1‑35 (MQ=255)
                             tCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaattgtcaaca                      <  2:541429‑M1/50‑17 (MQ=255)
                             tCCGCTCACAATTCCACACATTATACGAGCCGGAtgattaattgtcaac                       <  2:4974539‑M1/49‑16 (MQ=255)
                             tCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAg                      <  1:2972418/50‑1 (MQ=255)
                              ccGCTCACAATTTCACACAACATACGACCCGGAAGCATAAAGTGTAAAGc                     >  1:4848672/1‑50 (MQ=255)
                               cGCTCACAATTCCACACATTATACGAGCCGGAtgattaattgtcaacagc                    <  2:4146055‑M1/50‑19 (MQ=255)
                               cGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGcc                    <  1:1250603/50‑1 (MQ=255)
                                gCTCACAATTCCACACATTATACGAGCCGGAtgattaattgtcaacagct                   >  2:845616‑M1/1‑31 (MQ=255)
                                gCTCACAATTCCACACATTATACGAGCCGGAtgattaattgtcaacagct                   >  2:251959‑M1/1‑31 (MQ=255)
                                gCTCACAATTCCACACATTATACGAGCCGGAtgattaattgtcaacagct                   <  2:178512‑M1/50‑20 (MQ=255)
                                gCTCACAATTCCACACATTATACGAGCCGGAtgattaattgtcaacagct                   <  1:3686770‑M1/50‑20 (MQ=255)
                                 cTCACAATTCCACACATTATACGAGCCGGAtgattaattgtcaacagct                   >  1:3813261‑M1/1‑30 (MQ=255)
                                   cacaATTCCACACATTATACGAGCCGGAtgattaattgtcaacagctcat                >  2:2988714‑M1/1‑28 (MQ=255)
                                    acaATTCCACACATTATACGAGCCGGAtgattaattgtcaacagctcatt               <  1:3995780‑M1/50‑24 (MQ=255)
                                    acaATTCCACACATTATACGAGCCGGAtgattaattgtcaacagctcatt               >  2:3323130‑M1/1‑27 (MQ=255)
                                    acaATTCCACACATTATACGAGCCGGAtgattaattgtcaacagctcat                <  1:2831043‑M1/49‑23 (MQ=255)
                                    acaATTCCACACATTATACGAGCCGGAtgattaattgtcaacagctcat                <  2:1175129‑M1/49‑23 (MQ=255)
                                    acaATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTgggg               >  1:1678691/1‑50 (MQ=255)
                                    acaATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTgggg               >  2:4730929/1‑50 (MQ=255)
                                     caATTCCACGCAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGt              <  1:3528701/50‑1 (MQ=255)
                                     caATTCCACACATTATACGAGCCGGAtgattaattgtcaacagctcattt              >  2:683328‑M1/1‑26 (MQ=255)
                                     caATTCCACACATTATACGAGCCGGAtgattaattgtcaacagctcattt              <  2:4782933‑M1/50‑25 (MQ=255)
                                     caATTCCACACATTATACGAGCCGGAtgattaattgtcaacagctcattt              >  1:2305401‑M1/1‑26 (MQ=255)
                                     caATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGt              <  2:1065608/50‑1 (MQ=255)
                                     caATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGt              >  1:4705130/1‑50 (MQ=255)
                                      aaTTCCACACATTATACGAGCCGGAtgattaattgtcaacagctcatttc             <  1:1726308‑M1/50‑26 (MQ=255)
                                      aaTTCCACACATTATACGAGCCGGAtgattaattgtcaacagctcatttc             <  2:4238484‑M1/50‑26 (MQ=255)
                                      aaTTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTg             >  1:785433/1‑50 (MQ=255)
                                      aaTTCCACACAACATACCATCCGGCAGCATACAGTGTAAAGCCTggggag             >  2:232096/1‑48 (MQ=25)
                                       aTTCGACACATTATACGAGCCGGAtgattacttgtcaacagctcattt              >  2:3114441‑M1/1‑24 (MQ=255)
                                       aTTCCACACATTATACGAGCCGGAtgattaattgtcaacagctcatttca            <  2:275738‑M1/50‑27 (MQ=255)
                                       aTTCCACACATTATACGAGCCGGAtgattaattgtcaacagctcattt              >  2:1380188‑M1/1‑24 (MQ=255)
                                       aTTCCACACATTATACGAGCCGGAtgattaattgtcaacagctcattt              >  1:3204782‑M1/1‑24 (MQ=255)
                                        ttCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGcc           <  1:1932866/50‑1 (MQ=255)
                                         tCCACACATTATACGAGCCGGAtgattaattgtcaacagctcctttcag           >  2:4042721‑M1/1‑22 (MQ=255)
                                         tCCACACATTATACGAGCCGGAtgattaattgtcaacagctcatttcag           >  1:4162764‑M1/1‑22 (MQ=255)
                                         tCCACACATTATACGAGCCGGAtgattaattgtcaacagctcatttcag           >  2:3365461‑M1/1‑22 (MQ=255)
                                         tCCACACATTATACGAGCCGGAtgattaattgtcaacagctcatttcag           >  1:1297247‑M1/1‑22 (MQ=255)
                                          ccACACATTATACGAGCCGGAtgattaattgtcaacagctcatttcaga          <  1:3157526‑M1/49‑29 (MQ=255)
                                          ccACACATTATACGAGCCGGAtgattaattgtcaacagctcatttcaga          <  2:3201686‑M1/49‑29 (MQ=255)
                                           cacacaTTATACGAGCCGGAtgattaattgtcaacagctcatttcagaat        <  1:5117494‑M1/50‑31 (MQ=255)
                                           cacacaTTATACGAGCCGGAtgattaattgtcaacagctcatttcagaat        <  1:2918512‑M1/50‑31 (MQ=255)
                                           cacacaTTATACGAGCCGGAtgattaattgtcaacagctcatttcagaa         >  1:105612‑M1/1‑20 (MQ=255)
                                           cacacaTTATACGAGCCGGAtgattaattgtcaacagctcatttcagaa         >  2:1335664‑M1/1‑20 (MQ=255)
                                           cacacaACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTaa        >  2:2046311/1‑50 (MQ=255)
                                            acacaTTATACGAGCCGGAtgattaattgtcaacagctcatttcagaata       >  1:1556362‑M1/1‑19 (MQ=255)
                                              acaacaTACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAAt       >  1:5033802/1‑48 (MQ=255)
                                                aacaTACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAAtgagt   <  1:38434/50‑1 (MQ=255)
                                                 acaTGCGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAAtgagtg  <  2:1296052/50‑1 (MQ=255)
                                                 |                                                 
TCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTG  >  NC_000913/366302‑366400

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