Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913312,7670GA46.0% 34.2 / inf 388intergenic (‑766/‑814)ykgM/eaeHputative ribosomal protein/attaching and effacing protein
Reads supporting (aligned to +/- strand):  ref base G (104/105);  new base A (79/99);  total (184/204)
Fisher's exact test for biased strand distribution p-value = 3.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 398 total aligned reads displayed.

GCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAAC  >  NC_000913/312692‑312842
                                                                           |                                                                           
 cATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGa                                                                            >  1:2960631/1‑76 (MQ=255)
  atatCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGaa                                                                           >  2:4632994/1‑76 (MQ=255)
   tatCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAAt                                                                          <  2:4445381/76‑1 (MQ=255)
   tatCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAAt                                                                          >  1:5141492/1‑76 (MQ=255)
   tatCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAAt                                                                          <  2:3505961/76‑1 (MQ=255)
    atCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAAtt                                                                         <  2:641954/76‑1 (MQ=255)
    atCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAAtt                                                                         <  2:3238627/76‑1 (MQ=255)
    agcGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAAtt                                                                         <  2:328340/74‑1 (MQ=255)
      cGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGa                                                                       >  2:153033/1‑76 (MQ=255)
      cGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGa                                                                       <  1:1637772/76‑1 (MQ=255)
      cGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGa                                                                       >  1:4636852/1‑76 (MQ=255)
       gCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAg                                                                      >  2:5534731/1‑76 (MQ=255)
       gCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAg                                                                      <  1:3416067/76‑1 (MQ=255)
       gCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAg                                                                      >  1:1191856/1‑76 (MQ=255)
         aTCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGtt                                                                    >  1:3860835/1‑76 (MQ=255)
         aTCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGtt                                                                    >  1:3551941/1‑76 (MQ=255)
          tctcATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAAtt                                                                         <  1:423009/70‑1 (MQ=255)
          tctcATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTc                                                                   <  2:1199893/76‑1 (MQ=255)
          tctcATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTc                                                                   <  2:4179990/76‑1 (MQ=255)
          tctcATCATAATGTCGCATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTc                                                                   >  2:2506476/1‑76 (MQ=255)
            tcatcaTAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTa                                                                 <  1:1925325/76‑1 (MQ=255)
            tcatcaTAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTa                                                                 >  1:1409460/1‑76 (MQ=255)
              atcatAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTAtt                                                               >  1:2867256/1‑76 (MQ=255)
              atcatAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTAtt                                                               >  1:2813994/1‑76 (MQ=255)
              atcatAATGTCGTATCTCCTGGGGTGTTATACAAGACATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTAtt                                                               >  1:2742796/1‑76 (MQ=255)
               tcatAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTa                                                              >  1:269887/1‑76 (MQ=255)
                catAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTaa                                                             <  1:863277/76‑1 (MQ=255)
                 atAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTaaa                                                            <  2:358914/76‑1 (MQ=255)
                 atAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTaaa                                                            >  1:785677/1‑76 (MQ=255)
                 atAATGCCCTATCTACTGGGGGGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATGGAGTTCTATTaca                                                            >  2:3642165/1‑74 (MQ=255)
                  tAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAAc                                                           >  1:5370072/1‑76 (MQ=255)
                   aaTGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAAcc                                                          >  2:5684170/1‑76 (MQ=255)
                   aaTGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAAcc                                                          <  2:5518333/76‑1 (MQ=255)
                    aTGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCg                                                         <  1:5786498/76‑1 (MQ=255)
                    aTGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCg                                                         <  2:601869/76‑1 (MQ=255)
                    aTGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCg                                                         <  2:4572911/76‑1 (MQ=255)
                     tGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTTTATTTAAACGt                                                        >  1:2009659/1‑76 (MQ=255)
                     tGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGt                                                        >  1:3287355/1‑76 (MQ=255)
                      gTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTaaa                                                            <  2:3026106/71‑1 (MQ=255)
                      gTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTc                                                       >  1:247887/1‑76 (MQ=255)
                      gTCGTAGCTCATGGGGGGTTATACAAGATATCGTTGGTGGGGACGTGGGAGAGGAATTGAGTTCTATTAAACCGTc                                                       >  1:414989/1‑76 (MQ=255)
                       tCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCa                                                      >  1:4275277/1‑76 (MQ=255)
                       tCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCa                                                      >  2:1055647/1‑76 (MQ=255)
                       tCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCa                                                      <  2:5147670/76‑1 (MQ=255)
                         gTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAAc                                                    >  1:4076339/1‑76 (MQ=255)
                          tATCTCCTGGGGTGTTATACACGATCTCGTTGTTGTTGACCTGGGCTAGGAATTGATTTCTATTAAACCGTCAACt                                                   <  2:2840913/76‑1 (MQ=255)
                          tATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACt                                                   <  2:541183/76‑1 (MQ=255)
                            tctcCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTAt                                                 >  1:3743093/1‑76 (MQ=255)
                            tctcCTGGGGTGTTATACAAGATATCGTTGTTGGAGAAATGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTAt                                                 >  2:1844585/1‑76 (MQ=255)
                              tcCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGc                                               <  2:5341285/76‑1 (MQ=255)
                               ccTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGcc                                              >  2:557524/1‑76 (MQ=255)
                                cTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCg                                             <  2:3145793/76‑1 (MQ=255)
                                    ggTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGata                                         >  1:438336/1‑76 (MQ=255)
                                    ggTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGata                                         >  2:5797111/1‑76 (MQ=255)
                                    ggTGTTAGAAAAGTTAGCGTTGGGGGTGACCTGGGAGAGAAATTGAGTCATATGAAACAGGCAACGATGCCGGata                                         >  2:1757071/1‑76 (MQ=255)
                                     gtgtTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTGTGCCGGatac                                        <  2:853747/76‑1 (MQ=255)
                                      tgtTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGataca                                       >  1:4559431/1‑76 (MQ=255)
                                      tgtTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGataca                                       <  2:4330327/76‑1 (MQ=255)
                                      tgtTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGataca                                       <  2:2277019/76‑1 (MQ=255)
                                       gtTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGatacat                                      <  1:316393/76‑1 (MQ=255)
                                       gtTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGatacat                                      >  1:2027003/1‑76 (MQ=255)
                                         tataCAAGATCTCTTTGTTGGTGCCCTGCGAGAGAAATTGAGTTCTCTTAAACCGTCAACTATGCCGGatacctac                                    <  1:3898360/76‑1 (MQ=255)
                                         tataCAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGatacatac                                    <  1:5053186/76‑1 (MQ=255)
                                          ataCAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACt                                   <  1:5414486/76‑1 (MQ=255)
                                            aCAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTgg                                 >  1:4844877/1‑76 (MQ=255)
                                             cAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGa                                >  1:4691385/1‑76 (MQ=255)
                                             ataGATATCGTGGGTAGTGACCTGGGAGAGGAATTGAGTGCTATTGAACAGCATACTATTCCCGATACAGACTGGa                                >  2:3765402/3‑76 (MQ=255)
                                               aGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGAtt                              >  1:4945489/1‑76 (MQ=255)
                                                gATATCGTTGTTGGTGACCTGGGAGAGGACTTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGAttc                             <  2:2304400/76‑2 (MQ=255)
                                                 atatCGTTGTTGGTGTCCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTac                            <  2:208966/76‑1 (MQ=255)
                                                 atatCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCTTCAACTATGCCGGATACATACTGGATTac                            >  1:3619767/1‑76 (MQ=255)
                                                 atatCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTac                            <  2:748733/76‑1 (MQ=255)
                                                  tatCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTaca                           <  1:1571713/76‑1 (MQ=255)
                                                    tCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACt                         <  1:2336618/76‑1 (MQ=255)
                                                     cGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTg                        <  1:4747960/76‑1 (MQ=255)
                                                      gttgttGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGc                       <  2:5033864/76‑1 (MQ=255)
                                                       ttgttgGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCa                      <  2:786778/76‑1 (MQ=255)
                                                        tgttgGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAg                     <  2:1984296/76‑1 (MQ=255)
                                                        tgttgGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAg                     <  1:5729414/76‑1 (MQ=255)
                                                        tgttgGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAg                     >  2:4918262/1‑76 (MQ=255)
                                                        agttgGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAg                     >  2:1627307/2‑76 (MQ=255)
                                                         gttgGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAgg                    <  1:5554456/76‑1 (MQ=255)
                                                         gttgGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAgg                    <  2:3960019/76‑1 (MQ=255)
                                                          ttgGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGc                   <  2:1256823/76‑1 (MQ=255)
                                                            ggTTACCTGGGAGAGAAAGTGAGTTCTAGTAAACCGTCAACTATGCCGGATACATACTGGAGTACACTGCAGGCAc                 <  1:3079652/76‑1 (MQ=255)
                                                            ggTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCAc                 <  2:1384939/76‑1 (MQ=255)
                                                                aCCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCt             <  1:600711/76‑1 (MQ=255)
                                                                aCCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCt             <  1:497314/76‑1 (MQ=255)
                                                                 ccTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCtt            <  1:2218802/76‑1 (MQ=255)
                                                                 ccTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCtt            <  1:1030479/76‑1 (MQ=255)
                                                                  cTGGGAGCGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTa           <  1:2113513/76‑1 (MQ=255)
                                                                  cTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTa           <  2:2960631/76‑1 (MQ=255)
                                                                  cTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTa           >  1:1296750/1‑76 (MQ=255)
                                                                  cTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTa           <  2:2653860/76‑1 (MQ=255)
                                                                      gagagaAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATgag       >  1:3025413/1‑76 (MQ=255)
                                                                      gagagaAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATgag       >  1:659586/1‑76 (MQ=255)
                                                                       agagGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATgaga      <  2:290148/76‑1 (MQ=255)
                                                                       agagAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATgaga      <  2:4847004/76‑1 (MQ=255)
                                                                        gagGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATgagag     >  2:877930/1‑76 (MQ=255)
                                                                         agGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATgagaga    >  1:1448655/1‑76 (MQ=255)
                                                                           |                                                                           
GCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAAC  >  NC_000913/312692‑312842

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