Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 987,561 G→C 17.8% intergenic (‑579/+24) ompF ← / ← asnS outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,5610GC17.8% 417.2 / 81.5 164intergenic (‑579/+24)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (59/74);  new base C (16/13);  total (76/88)
Fisher's exact test for biased strand distribution p-value = 3.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

Only 100 of 176 total aligned reads displayed.

TATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCA  >  NC_000913/987462‑987660
                                                                                                   |                                                                                                   
  ttCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgcg                                                                                                   >  1:1647/1‑100 (MQ=255)
    ccGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCt                                                                                                 >  2:2038740/1‑100 (MQ=255)
    ccGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCt                                                                                                 >  2:2516762/1‑100 (MQ=255)
     cGTAACAGGATGATGGTTATGGCGGGGTTGCAAATAGATTGCTTAACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCt                                                                                                 <  1:3575749/99‑1 (MQ=255)
     cGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGAATTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTg                                                                                                >  1:2433281/1‑100 (MQ=255)
      gTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGTTTGAGTTAATAAAAAGATAACTTTATAGAGTAGGTAAAAAAAGAGCAAGCGGGCGCt                                                                                                 >  1:15702/1‑98 (MQ=255)
       tAACAGGATGATCGTTATGGCGGGGTTGCAAATAGACTGCTTGACTTAACAAACAGATAACTTGAAAGAAAAGATAAAAAAAGAGCCAGCGAgcgc                                                                                                  >  1:2002620/1‑96 (MQ=255)
         aCAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGctc                                                                                            <  2:2632226/100‑1 (MQ=255)
          cAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgcgagatcgga                                                                                           >  1:2311802/1‑92 (MQ=255)
          cAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGctct                                                                                           <  1:306788/100‑1 (MQ=255)
           aGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCtt                                                                                          <  1:1362266/100‑1 (MQ=255)
            ggATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTg                                                                                         >  2:3849189/1‑100 (MQ=255)
            ggATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTg                                                                                         >  2:2486466/1‑100 (MQ=255)
             gatgatCGTTATGGCGGGGTTGCAAATAGATTGCTTTACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGa                                                                                        <  1:1062795/100‑1 (MQ=255)
              atgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGaa                                                                                       <  2:297707/100‑1 (MQ=255)
               tgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAg                                                                                      <  2:1057606/100‑1 (MQ=255)
                gatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGa                                                                                     >  1:493315/1‑100 (MQ=255)
                gatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGa                                                                                     <  1:972306/100‑1 (MQ=255)
                gatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCttttttt                                                                                     <  2:2193246/100‑6 (MQ=38)
                 atCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGa                                                                                     >  1:261790/1‑99 (MQ=255)
                 atCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGAc                                                                                    <  2:555377/100‑1 (MQ=255)
                 atCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCttttttta                                                                                    >  2:3447680/1‑94 (MQ=38)
                 atCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCttttttat                                                                                    <  2:3320606/100‑7 (MQ=38)
                  tCGTTATGGCGGGGTTGCAAATAGATTGCTTTACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGa                                                                                        <  2:656638/95‑1 (MQ=255)
                  tCGGTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCtttttttat                                                                                   <  1:1890863/100‑8 (MQ=37)
                    gTTATGGCGGGGTTGCAAATAGATTGCTGGACTTAACAAACAGATAAATTGACAGAAAAGGTAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGAAGCt                                                                                 >  1:3444715/1‑100 (MQ=255)
                     ttATGGCGGGGTTGCAAATAGATTGCTTGACTAAACAAACAGATAACTTGACAGAAAATATAAAAAAAGAGCCAGCGCCCGCTGGCTCtttttttatc                                                                                  <  1:3804367/98‑9 (MQ=37)
                       aTGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTaaa                                                                              <  1:1186854/100‑1 (MQ=255)
                          gCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGCAGACGCTAAATTa                                                                           >  1:2180610/1‑100 (MQ=255)
                          gCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTa                                                                           <  2:4155453/100‑1 (MQ=255)
                            ggggTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGa                                                                         <  1:2219212/100‑1 (MQ=255)
                              gtttGTAAATAGATTGCTTGACTTAATAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAg                                                                       <  1:2101583/98‑1 (MQ=255)
                                   cAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCAAAATTAGAAGCTGGc                                                                  >  2:536366/1‑100 (MQ=255)
                                    aaaTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCg                                                                 <  1:3557320/100‑1 (MQ=255)
                                     aaTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGt                                                                <  1:4381820/100‑1 (MQ=255)
                                         gATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcg                                                            >  2:1860865/1‑100 (MQ=255)
                                         gATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcg                                                            >  1:2365883/1‑100 (MQ=255)
                                         gATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcg                                                            >  1:3151997/1‑100 (MQ=255)
                                          aTTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcgc                                                           >  1:2773466/1‑100 (MQ=255)
                                            tGCTTGACTCAACAAACAGATAACTTGACAGAAAAGCTAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcgcg                                                          >  2:1900209/1‑99 (MQ=255)
                                             gCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcgcg                                                          >  1:3946066/1‑98 (MQ=255)
                                              cTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGa                                                        >  1:212309/1‑99 (MQ=255)
                                              cTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAg                                                       <  1:271718/100‑1 (MQ=255)
                                              cTTGACTTAACAAACAGAGAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAg                                                       <  1:3849210/100‑1 (MQ=255)
                                                  aCTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACg                                                   <  1:4081365/100‑1 (MQ=255)
                                                  aCTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACg                                                   >  2:4064716/1‑100 (MQ=255)
                                                     tAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTgg                                                <  2:2227041/100‑1 (MQ=255)
                                                     tAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTgg                                                <  2:2227080/100‑1 (MQ=255)
                                                     tAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTgg                                                <  2:849240/100‑1 (MQ=255)
                                                     tAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTgg                                                <  1:1941223/100‑1 (MQ=255)
                                                      aacaaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTggg                                               <  1:2516003/100‑1 (MQ=255)
                                                         aaacaGATAACTTGACAGAAAGGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAAc                                            >  1:4217799/1‑100 (MQ=255)
                                                          aacaGATAACTTGACAGAAAAGATAAAAAAAGAGCAAACGAACGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAAc                                            >  1:2039302/1‑99 (MQ=255)
                                                           acaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACgg                                          >  2:3336656/1‑100 (MQ=255)
                                                             aGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGaa                                        <  1:829582/100‑1 (MQ=255)
                                                             aGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGaa                                        <  2:2078796/100‑1 (MQ=255)
                                                              gATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAat                                       >  2:3331616/1‑100 (MQ=255)
                                                              gATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAat                                       >  2:3878163/1‑100 (MQ=255)
                                                               aTAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatc                                      <  2:739979/100‑1 (MQ=255)
                                                               aTAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatc                                      <  2:266283/100‑1 (MQ=255)
                                                               aTAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatc                                      <  2:3566554/100‑1 (MQ=255)
                                                               aTAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatc                                      <  2:2718670/100‑1 (MQ=255)
                                                                 atcTTGACATACTAGATAAAAAAAGAGTCATCGGTCGCTTGCTCTTGAAGGCTCTAAATTAGCAGCTTGCGTTACGCGGATTACGTGTGAACGGAatca                                     <  1:3896006/97‑1 (MQ=255)
                                                                 aaCTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcac                                    <  2:4273735/100‑1 (MQ=255)
                                                                  aCTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcaca                                   <  1:2822763/100‑1 (MQ=255)
                                                                  aCTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcaca                                   >  2:1939467/1‑100 (MQ=255)
                                                                   cTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacat                                  >  1:3683398/1‑100 (MQ=255)
                                                                    ttGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatc                                 <  2:3110966/100‑1 (MQ=255)
                                                                     tGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatca                                <  2:1355060/100‑1 (MQ=255)
                                                                        cAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgt                             >  1:688194/1‑100 (MQ=255)
                                                                          gAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtac                           >  1:2864901/1‑100 (MQ=255)
                                                                          gAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtac                           >  2:487467/1‑100 (MQ=255)
                                                                           aaaaGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtacg                          >  2:1408041/1‑100 (MQ=255)
                                                                              aGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGttt                       <  2:1815763/100‑1 (MQ=255)
                                                                               gATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGtttt                      <  2:1700709/100‑1 (MQ=255)
                                                                                cgtcttcaaGAGCCAGCGCCCGGTGGCTCTTGAAGACGCTAACTTAGAAGCTGGTGTTACGCGGAGTAGGTGGGAAAGGAatcacatcac                               >  2:2176480/8‑90 (MQ=25)
                                                                                cgtcttcaaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTg                     >  2:4205613/8‑100 (MQ=38)
                                                                                aTAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTg                     <  1:1106661/100‑1 (MQ=255)
                                                                                aTAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTg                     <  2:2889545/100‑1 (MQ=255)
                                                                                aTAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTg                     <  2:2298569/100‑1 (MQ=255)
                                                                                 tAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                    <  2:3514242/100‑1 (MQ=255)
                                                                                  tcttcaaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCa                   >  2:3898748/6‑100 (MQ=38)
                                                                                    ttcaaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACg                 >  2:2445088/4‑100 (MQ=39)
                                                                                    aaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACg                 <  2:3238550/100‑1 (MQ=255)
                                                                                    aaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACg                 <  1:627813/100‑1 (MQ=255)
                                                                                     aaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGATGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGc                >  2:2400758/1‑100 (MQ=255)
                                                                                      caaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGcc               >  1:533063/2‑100 (MQ=255)
                                                                                      aaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGcc               <  2:2588902/100‑1 (MQ=255)
                                                                                       aaGAGCCAGCGCCCGCTGGCTCTTGAAGGCGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGCACGGAATCACATCACGTACGTTTTCCACGCCa              >  1:750140/1‑100 (MQ=255)
                                                                                       aaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCa              >  2:1348341/1‑100 (MQ=255)
                                                                                        agagCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCa              >  1:300616/1‑99 (MQ=255)
                                                                                        agagCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCTTTCCGCGGAGTACGTGGGAACGTAATCACATCACGTACGAATTACACGCCAg             >  2:38071/1‑100 (MQ=255)
                                                                                        agagCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAg             >  1:2441215/1‑100 (MQ=255)
                                                                                         gagCCAGCGCCCGCTGGCTCTTGAAAACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGt            >  2:3783940/1‑100 (MQ=255)
                                                                                           gCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTa          >  1:501968/1‑100 (MQ=255)
                                                                                              aGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGt       <  1:3480626/100‑1 (MQ=255)
                                                                                               gCGGGCGCTGGCTCTTGAAGACTCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACg        <  2:2326940/98‑1 (MQ=255)
                                                                                               gCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTa      <  1:3985826/100‑1 (MQ=255)
                                                                                                 gggCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAg    >  2:1291975/1‑100 (MQ=255)
                                                                                                 gggCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAg    >  2:2993051/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCA  >  NC_000913/987462‑987660

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: