Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 987,560 G→C 17.4% intergenic (‑578/+25) 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,5600GC17.4% 434.2 / 82.3 170intergenic (‑578/+25)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (62/77);  new base C (16/13);  total (80/90)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-01

Only 100 of 175 total aligned reads displayed.

TATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGC  >  NC_000913/987462‑987659
                                                                                                  |                                                                                                   
  ttCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgcg                                                                                                  >  1:1647/1‑100 (MQ=255)
  ttCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgcg                                                                                                  >  2:3889357/1‑100 (MQ=255)
  tccgatctcAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCTGCGGgcg                                                                                                  <  2:2311812/92‑1 (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)
         aCAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGctc                                                                                           <  1:498122/100‑1 (MQ=255)
          cAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgcgagatcgga                                                                                          >  1:2311802/1‑92 (MQ=255)
           aGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCtt                                                                                         <  1:1362266/100‑1 (MQ=255)
            ggATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTg                                                                                        >  2:2486466/1‑100 (MQ=255)
             gatgatCGTTATGGCGGGGTTGCAAATAGATTGCTTTACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGa                                                                                       <  1:1062795/100‑1 (MQ=255)
              atgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCttttt                                                                                      <  2:2489128/100‑4 (MQ=39)
               tgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAg                                                                                     <  2:1057606/100‑1 (MQ=255)
                gatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGa                                                                                    <  1:972306/100‑1 (MQ=255)
                gatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGa                                                                                    >  1:493315/1‑100 (MQ=255)
                 atCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGa                                                                                    >  1:261790/1‑99 (MQ=255)
                 atCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCttttttta                                                                                   <  2:3805829/100‑7 (MQ=38)
                  tCGTTATGGCGGGGTTGCAAATAGATTGCTTTACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGa                                                                                       <  2:656638/95‑1 (MQ=255)
                  tCGGTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCtttttttat                                                                                  <  1:1890863/100‑8 (MQ=37)
                   cGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGc                                                                                 <  2:2526840/100‑1 (MQ=255)
                   cGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCtttttttatc                                                                                 <  2:1738395/100‑9 (MQ=38)
                    gTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCt                                                                                <  2:4058111/100‑1 (MQ=255)
                    gTTATGGCGGGGTTGCAAATAGATTGCTGGACTTAACAAACAGATAAATTGACAGAAAAGGTAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGAAGCt                                                                                >  1:3444715/1‑100 (MQ=255)
                     ttATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTa                                                                               <  2:1924993/100‑1 (MQ=255)
                     ttATGGCGGGGTTGCAAATAGATTGCTTGACTAAACAAACAGATAACTTGACAGAAAATATAAAAAAAGAGCCAGCGCCCGCTGGCTCtttttttatc                                                                                 <  1:3804367/98‑9 (MQ=37)
                       aTGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTaaa                                                                             <  1:1186854/100‑1 (MQ=255)
                       aTGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTaaa                                                                             <  2:1258235/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)
                              ggTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAg                                                                      >  2:2690249/1‑100 (MQ=255)
                               gTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGc                                                                     >  2:344980/1‑100 (MQ=255)
                                   cAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGc                                                                 >  1:3048825/1‑100 (MQ=255)
                                   cAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCAAAATTAGAAGCTGGc                                                                 >  2:536366/1‑100 (MQ=255)
                                     aaTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGt                                                               <  1:4381820/100‑1 (MQ=255)
                                         gATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcg                                                           >  1:3151997/1‑100 (MQ=255)
                                         gATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcg                                                           >  2:1860865/1‑100 (MQ=255)
                                         gATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcg                                                           >  1:2365883/1‑100 (MQ=255)
                                          aTTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcgc                                                          >  1:2773466/1‑100 (MQ=255)
                                            tGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                                        <  2:3244959/100‑1 (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)
                                               ttGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGt                                                     >  2:1294442/1‑100 (MQ=255)
                                                  aCTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACg                                                  >  2:2327244/1‑100 (MQ=255)
                                                     tAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTgg                                               <  2:849240/100‑1 (MQ=255)
                                                     tAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTgg                                               <  2:2227080/100‑1 (MQ=255)
                                                     tAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTgg                                               <  1:1941223/100‑1 (MQ=255)
                                                      aacaaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTggg                                              <  1:2516003/100‑1 (MQ=255)
                                                      aacaaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTggg                                              >  2:1973494/1‑100 (MQ=255)
                                                        caaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGaa                                            >  2:3063028/1‑100 (MQ=255)
                                                         aaacaGATAACTTGACAGAAAGGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAAc                                           >  1:4217799/1‑100 (MQ=255)
                                                         aaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAAc                                           <  2:4270089/100‑1 (MQ=255)
                                                         aaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAAc                                           >  2:3415254/1‑100 (MQ=255)
                                                          aacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACg                                          <  2:3334502/100‑1 (MQ=255)
                                                          aacaGATAACTTGACAGAAAAGATAAAAAAAGAGCAAACGAACGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAAc                                           >  1:2039302/1‑99 (MQ=255)
                                                             aGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGaa                                       <  2:2078796/100‑1 (MQ=255)
                                                             aGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGaa                                       <  1:829582/100‑1 (MQ=255)
                                                             aGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGaa                                       >  2:3640943/1‑100 (MQ=255)
                                                              gATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAat                                      <  2:3124603/100‑1 (MQ=255)
                                                              gATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAat                                      >  2:3331616/1‑100 (MQ=255)
                                                              gATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAat                                      >  2:3878163/1‑100 (MQ=255)
                                                               aTAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatc                                     <  2:3814922/100‑1 (MQ=255)
                                                               aTAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatc                                     <  2:266283/100‑1 (MQ=255)
                                                               aTAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatc                                     <  2:3566554/100‑1 (MQ=255)
                                                               aTAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatc                                     <  2:739979/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)
                                                                   cTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacat                                 >  1:3683398/1‑100 (MQ=255)
                                                                   cTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacat                                 >  1:4050505/1‑100 (MQ=255)
                                                                     tGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatca                               >  1:72966/1‑100 (MQ=255)
                                                                     tGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatca                               <  2:1355060/100‑1 (MQ=255)
                                                                        cAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgt                            >  1:688194/1‑100 (MQ=255)
                                                                          gAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtac                          >  2:487467/1‑100 (MQ=255)
                                                                          gAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtac                          >  1:2864901/1‑100 (MQ=255)
                                                                           aaaaGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtacg                         >  2:1408041/1‑100 (MQ=255)
                                                                              aGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGttt                      <  2:2864917/100‑1 (MQ=255)
                                                                                cgtcttcaaGAGCCAGCGCCCGGTGGCTCTTGAAGACGCTAACTTAGAAGCTGGTGTTACGCGGAGTAGGTGGGAAAGGAatcacatcac                              >  2:2176480/8‑90 (MQ=25)
                                                                                cgtcttcaaGAGCCAGCGCCCGCTGGCTCTTTAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTg                    >  2:643757/8‑100 (MQ=37)
                                                                                aTAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTg                    <  1:1106661/100‑1 (MQ=255)
                                                                                aTAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTg                    <  2:2298569/100‑1 (MQ=255)
                                                                                 tgaaaaaaGAGCCAGCGGGCGCTGGCTCTTGAAAACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAAGCACATCACGTACGTTTTGc                   >  2:2976765/3‑100 (MQ=255)
                                                                                 tAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                   <  1:3499692/100‑1 (MQ=255)
                                                                                 tAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                   <  2:2187655/100‑1 (MQ=255)
                                                                                    ttcaaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACg                >  2:2445088/4‑100 (MQ=39)
                                                                                    aaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACg                <  1:627813/100‑1 (MQ=255)
                                                                                      caaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGcc              >  1:533063/2‑100 (MQ=255)
                                                                                      aaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGcc              <  2:2226350/100‑1 (MQ=255)
                                                                                        agagCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCa             >  1:300616/1‑99 (MQ=255)
                                                                                        agagCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCCTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAg            >  1:3827881/1‑100 (MQ=255)
                                                                                        agagCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAg            >  1:2441215/1‑100 (MQ=255)
                                                                                        agagCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAg            >  2:162523/1‑100 (MQ=255)
                                                                                              aGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGt      <  1:3480626/100‑1 (MQ=255)
                                                                                               gCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTa     <  1:3985826/100‑1 (MQ=255)
                                                                                                 gggCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAg   >  2:2993051/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
TATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGC  >  NC_000913/987462‑987659

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: