Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,0140AC20.6% 332.3 / 11.2 97intergenic (‑32/+571)ompF/asnSouter membrane porin F/asparagine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (60/17);  new base C (0/20);  total (60/37)
Fisher's exact test for biased strand distribution p-value = 5.84e-11
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 111 total aligned reads displayed.

CATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATA  >  NC_000913/986980‑987046
                                  |                                
catTATTTATTACCCTCATGGTTTTTTTTATGaca                                  <  1:3413120/35‑1 (MQ=255)
catTATTTATTACCCTCATGGTTTTTTTTATGaca                                  <  1:1213283/35‑1 (MQ=255)
 attattTATTACCCTCATGGTTTTTTTTATGacac                                 >  1:800324/1‑35 (MQ=255)
 attattTATTACCCTCATGGTTTTTTTTATGacac                                 <  1:272791/35‑1 (MQ=255)
 attattTATTACCCTCATGGTTTTTTTTATGacac                                 >  1:7861240/1‑35 (MQ=255)
  ttatttatTACCCTCATGGTTTTTTTTATGACAcc                                <  1:1470759/35‑1 (MQ=255)
  ttatttatTACCCTCATGGTTTTTTTTATGACAcc                                <  1:4338455/35‑1 (MQ=255)
   tatttattACCCTCATGGTTTTTTTTATGACACCt                               <  1:7663194/35‑1 (MQ=255)
   tatttattACCCTCATGGTTTTTTTTATGACACCt                               <  1:6562171/35‑1 (MQ=255)
   tatttattACCCTCATGGTTTTTTTTATGACACCt                               <  1:4311284/35‑1 (MQ=255)
    atttattACCCTCATGGTTTTTTTTATGACACCTg                              <  1:4855551/35‑1 (MQ=255)
     tttattACCCTCATGGTTTTTTTTATGACACCTGc                             <  1:6898715/35‑1 (MQ=255)
     tttattACCCTCATGGTTTTTTTTATGACACCTGc                             >  1:6031187/1‑35 (MQ=255)
     tttattACCCTCATGGTTTTTTTTATGACACCTGc                             >  1:4836192/1‑35 (MQ=255)
     tttattACCCTCATGGTTTTTTTTATGACACCTGc                             >  1:290448/1‑35 (MQ=255)
     tttattACCCTCATGGTTTTTTTTATGACACCTGc                             >  1:1702600/1‑35 (MQ=255)
     tttattACCCTCATGGTTTTTTTTATGACACCTGc                             >  1:3711293/1‑35 (MQ=255)
     tttattACCCTCATGGTTTTTTTTATGACACCTGc                             <  1:7438115/35‑1 (MQ=255)
     tttattACCCTCATGGTTTTTTTTATGACACCTGc                             <  1:7742391/35‑1 (MQ=255)
     tttattACCCTCATGGTTTTTTTTATGACACCTGc                             >  1:7860279/1‑35 (MQ=255)
      ttattaCCCTCATGGTTTTTTTTATGACACCTGcc                            <  1:1735356/35‑1 (MQ=255)
      ttattaCCCTCATGGTTTTTTTTATGACACCTGcc                            <  1:4144631/35‑1 (MQ=255)
      ttattaCCCTCATGGTTTTTTTTATGACACCTGcc                            <  1:1995190/35‑1 (MQ=255)
      ttattaCCCTCATGGTTTTTTTTATGACACCTGcc                            <  1:6775220/35‑1 (MQ=255)
      ttattaCCCTCATGGTTTTTTTTATGACACCTGcc                            <  1:3950515/35‑1 (MQ=255)
       tattcCCCTCATGGTTTTTTTTATGACACCTGCCa                           <  1:6220087/35‑1 (MQ=255)
       tattcCCCTCATGGTTTTTTTTATGACACCTGCCa                           <  1:1424707/35‑1 (MQ=255)
       tattaCCCTCATGGTTTTTTTTATGACCCCTGCCa                           <  1:5122009/35‑1 (MQ=255)
       tattaCCCTCATGGTTTTTTTTATGACACCTGCCa                           >  1:7386146/1‑35 (MQ=255)
        attaCCCTCATGGTTTTTTTTATGACACCTGCCAc                          >  1:4433020/1‑35 (MQ=255)
        attaCCCTCATGGTTTTTTTTATGACACCTGCCAc                          >  1:6843758/1‑35 (MQ=255)
        attaCCCTCATGGTTTTTTTTATGACACCTGCCAc                          >  1:7337315/1‑35 (MQ=255)
        attaCCCTCATGGTTTTTTTTATGACACCTGCCAc                          >  1:339410/1‑35 (MQ=255)
         ttaCCCTCATGGTTTTTTTTATGACACCTGCCACt                         >  1:1230588/1‑35 (MQ=255)
         ttaCCCTCATGGTTTTTTTTATGACACCTGCCACt                         >  1:7079523/1‑35 (MQ=255)
           ccccTCATGTTTTTTTTTATGACACCTGCCACTGc                       <  1:6711255/34‑1 (MQ=38)
           aCCCTCATGGTTTTTTTTATGACACCTGCCACTGc                       >  1:5831899/1‑35 (MQ=255)
           aCCCTCATGGTTTTTTTTATGACACCTGCCACTGc                       >  1:2464193/1‑35 (MQ=255)
           aCCCTCATGGTTTTTTTTATGACACCTGCCACTGc                       >  1:3396459/1‑35 (MQ=255)
            cccTCATGGTTTTTTTTATGACCCCTGCCACTGcc                      <  1:3283618/35‑1 (MQ=255)
            cccTCATGGTTTTTTTTATGACACCTGCCACTGcc                      >  1:6851702/1‑35 (MQ=255)
            cccTCATGGTTTTTTTTATGACACCTGCCACTGcc                      >  1:916109/1‑35 (MQ=255)
            cccTCATGGTTTTTTTTATGACACCTGCCACTGcc                      >  1:2105598/1‑35 (MQ=255)
            cccTCATGGTTTTTTTTATGACACCTGCCACTGcc                      >  1:1367786/1‑35 (MQ=255)
            cccTCATGGTTTTTTTTATGACACCTGCCACTGcc                      >  1:1248360/1‑35 (MQ=255)
              cTCATGGTTTTTTTTATGACCCCTGCCACTGCCGt                    <  1:1763903/35‑1 (MQ=255)
              cTCATGGTTTTTTTTATGACACCTGCCACTGCCGt                    >  1:1044882/1‑35 (MQ=255)
               tCATGTTTTTTTTTATGACCCCTGCCACTGCCGTc                   <  1:4797209/35‑1 (MQ=38)
               tCATGGTTTTTTTTATGACCCCTGCCACTGCCGTc                   <  1:2491472/35‑1 (MQ=255)
               tCATGGTTTTTTTTATGACCCCTGCCACTGCCGTc                   <  1:1423703/35‑1 (MQ=255)
               tCATGGTTTTTTTTATGACACCTGCCACTGCCGTc                   >  1:7166897/1‑35 (MQ=255)
               tCATGGTTTTTTTTATGACACCTGCCACTGCCGTc                   >  1:2283831/1‑35 (MQ=255)
                 aTGGTTTTTTTTATGACCCCTGCCACTGCCGTCaa                 <  1:3179997/35‑1 (MQ=255)
                 aTGGTTTTTTTTATGACCCCTGCCACTGCCGTCaa                 <  1:2769457/35‑1 (MQ=255)
                 aTGGTTTTTTTTATGACACCTGCCACTGCCGTCaa                 >  1:4908037/1‑35 (MQ=255)
                 aTGGTTTTTTTTATGACACCTGCCACTGCCGTCaa                 >  1:2172819/1‑35 (MQ=255)
                 aTGGTTTTTTTTATGACACCTGCCACTGCCGTCaa                 >  1:2172861/1‑35 (MQ=255)
                 aTGGTTTTTTTTATGACACCTGCCACTGCCGTCaa                 >  1:7816568/1‑35 (MQ=255)
                  tGGTTTTTTTTATGACCCCTGCCACTGCCGTCAAt                <  1:1655435/35‑1 (MQ=255)
                  tGGTTTTTTTTATGACCCCTGCCACTGCCGTCAAt                <  1:6802440/35‑1 (MQ=255)
                    tttttttttATGACCCCTGCCACTGCCGTCAATaa              <  1:6427633/34‑1 (MQ=39)
                    gTTTTTTTTATGACACCTGCCACTGCCGTCAATaa              >  1:4671360/1‑35 (MQ=255)
                    gTTTTTTTTATGACACCTGCCACTGCCGTCAATaa              >  1:1792519/1‑35 (MQ=255)
                    gTTTTTTTTACGCCCCCTGCCACTGCCGTCAATaa              <  1:4074217/35‑1 (MQ=255)
                     ttttttttATGACACCTGCCACTGCCGTCAATAAg             >  1:1673600/1‑35 (MQ=255)
                      tttttttATGCCCCCTGCCACTGCCGTCAATAAGt            <  1:3229961/35‑1 (MQ=38)
                      tttttttATGACCCCTGCCACTGCCGTCAATAAGt            <  1:7129471/35‑1 (MQ=255)
                      tttttttATGACACCTGCCACTGCCGTCAATAAGt            >  1:4714777/1‑35 (MQ=255)
                       ttttttATGCCCCCTGCCACTGCCGTCAATAAGtt           <  1:6077142/35‑1 (MQ=38)
                       ttttttATGACACCTGCCACTGCCGTCAATAAGtt           >  1:1180322/1‑35 (MQ=255)
                       ttttttATGACACCTGCCACTGCCGTCAATAAGtt           >  1:4563898/1‑35 (MQ=255)
                       ttttttATGACACCTGCCACTGCCGTCAATAAGtt           >  1:2822582/1‑35 (MQ=255)
                        tttttATGACACCTGCCACTGCCGTCAATAAGTTc          >  1:1356792/1‑35 (MQ=255)
                        tttttATGACACCTGCCACTGCCGTCAATAAGTTc          >  1:6763773/1‑35 (MQ=255)
                        tttttATGACACCTGCCACTGCCGTCAATAAGTTc          >  1:6763728/1‑35 (MQ=255)
                         ttttATGACACCTGCCACTGCCGTCAATAAGTTCt         >  1:856689/1‑35 (MQ=255)
                          tttATGCCCCCTGCCACTGCCGTCAATAAGTTCTg        <  1:3605876/35‑1 (MQ=38)
                          tttATGACCCCTGCCACTGCCGTCAATAAGTTCTg        <  1:3605917/35‑1 (MQ=255)
                          tttATGACCCCTGCCACTGCCGTCAATAAGTTCTg        <  1:2401494/35‑1 (MQ=255)
                          tttATGACACCTGCCACTGCCGTCAATAAGTTCTg        >  1:7405940/1‑35 (MQ=255)
                           ttATGACCCCTGCCACTGCCGTCAATAAGTTCTGt       <  1:7069127/35‑1 (MQ=255)
                           ttATGACACCTGCCACTGCCGTCAATAAGTTCTGt       >  1:4150411/1‑35 (MQ=255)
                           ttATGACACCTGCCACTGCCGTCAATAAGTTCTGt       >  1:1716876/1‑35 (MQ=255)
                            tATGACCCCTGCCACTGCCGTCAATAAGTTCTGTc      <  1:1094371/35‑1 (MQ=255)
                            tATGACACCTGCCACTGCCGTCAATAAGTTCTGTc      >  1:3146536/1‑35 (MQ=255)
                            tATGACACCTGCCACTGCCGTCAATAAGTTCTGTc      >  1:2284331/1‑35 (MQ=255)
                            tATGACACCTGCCACTGCCGTCAATAAGTTCTGTc      >  1:2128523/1‑35 (MQ=255)
                             aTGACCCCTGCCACTGCCGTCAATAAGTTCTGTCa     <  1:1713988/35‑1 (MQ=255)
                              tGACCCCTGCCACTGCCGTCAATAAGTTCTGTCaa    <  1:5547377/35‑1 (MQ=255)
                              tGACCCCTGCCACTGCCGTCAATAAGTTCTGTCaa    <  1:6993880/35‑1 (MQ=255)
                              tGACCCCTGCCACTGCCGTCAATAAGTGCTGTCaa    <  1:8221/35‑1 (MQ=38)
                              tGACACCTGCCACTGCCGTCAATAAGTTCTGTCaa    >  1:6355808/1‑35 (MQ=255)
                              tGACACCTGCCACTGCCGTCAATAAGTTCTGTCaa    >  1:4944590/1‑35 (MQ=255)
                              tGACACCTGCCACTGCCGTCAATAAGTTCTGTCaa    >  1:7071648/1‑35 (MQ=255)
                              tGACACCTGCCACTGCCGTCAATAAGTTCTGTCaa    >  1:2019479/1‑35 (MQ=255)
                              tGACACCTGCCACTGCCGTCAATAAGTTCTGTCaa    >  1:2006973/1‑35 (MQ=255)
                               gACACCTGCCACTGCCGTCAATAAGTTCTGTCAAt   >  1:5063446/1‑35 (MQ=255)
                                acacCTGCCACTGCCGTCAATAAGTTCTGTCAATa  >  1:6274214/1‑35 (MQ=255)
                                acacCTGCCACTGCCGTCAATAAGTTCTGTCAATa  >  1:6483407/1‑35 (MQ=255)
                                acacCTGCCACTGCCGTCAATAAGTTCTGTCAATa  >  1:5224538/1‑35 (MQ=255)
                                  |                                
CATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATA  >  NC_000913/986980‑987046

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