Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 987,561 G→C 17.3% 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.3% 333.3 / 56.7 128intergenic (‑579/+24)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (46/58);  new base C (8/14);  total (55/73)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 135 total aligned reads displayed.

TATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGC  >  NC_000913/987462‑987659
                                                                                                   |                                                                                                  
tatTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCggg                                                                                                    >  1:745675/1‑100 (MQ=255)
 atTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgc                                                                                                   >  1:1861060/1‑100 (MQ=255)
      gTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAGGTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTgg                                                                                              <  1:486438/100‑1 (MQ=255)
       tAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGc                                                                                             >  1:1376194/1‑100 (MQ=255)
        aaCAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGct                                                                                            <  2:3578890/100‑1 (MQ=255)
         aCAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTgg                                                                                              >  2:1498026/1‑97 (MQ=255)
         aCAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGctc                                                                                           <  1:497239/100‑1 (MQ=255)
          cAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGctct                                                                                          <  1:3760366/100‑1 (MQ=255)
          cAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGctct                                                                                          <  2:2153977/100‑1 (MQ=255)
          cAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGctct                                                                                          <  1:777139/100‑1 (MQ=255)
          cAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGctct                                                                                          <  2:1429701/100‑1 (MQ=255)
           aGGATGATCGTGATGGCGGGGTTGCAAATAGATTGCGTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGctct                                                                                          <  1:3211209/99‑1 (MQ=255)
            ngatgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTg                                                                                        >  1:870087/2‑100 (MQ=255)
            ggATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCttn                                                                                        <  2:1797720/100‑2 (MQ=255)
            ggATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACGTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTg                                                                                        <  2:3738610/100‑1 (MQ=255)
            ggATGATCGTTATGGCGGGGTTGCAAATAGATGGCTTGACTTAACAAACAGGTAGCTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTg                                                                                        <  2:3059787/100‑1 (MQ=255)
             gatgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGctct                                                                                          <  2:2852592/97‑1 (MQ=255)
             gaggatCGTTATGGCGGGGTTGTAAATAGATTGCGTGACCTAACAAGCAGATAGCTTGACAGAAAAAATAAAAAAAGAGCCAGCGCCCGCTGGCTCtt                                                                                         <  2:310182/95‑1 (MQ=255)
              atgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGctctattt                                                                                      <  1:1293961/100‑5 (MQ=255)
              atgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCttttt                                                                                      <  1:3581162/100‑4 (MQ=39)
               tgatGGTTATGGCGGGGTTGCAAATAGATTGGTTGACTTAACAAACAGATAACTGGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTg                                                                                        <  2:3738558/97‑1 (MQ=255)
               tgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCtttttt                                                                                     <  2:572341/100‑5 (MQ=39)
                 atCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGAc                                                                                   >  2:32595/1‑100 (MQ=255)
                   cGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGc                                                                                 >  1:801001/1‑100 (MQ=255)
                      taagacGGGATTACAAATAGATACCTTCACTTAACACACAGATAACTTGACAGCACAGATAAAAAAAGAGCCAGCGGGAGCTGGCTCTTGAAGACCCt                                                                                <  2:3241349/93‑1 (MQ=255)
                        tGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAAt                                                                            <  1:2561257/100‑1 (MQ=255)
                        tGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAAt                                                                            <  1:233387/100‑1 (MQ=255)
                          gCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTa                                                                          >  1:2512344/1‑100 (MQ=255)
                          gCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTa                                                                          >  1:86716/1‑100 (MQ=255)
                             gggTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGaa                                                                       <  1:3292878/100‑1 (MQ=255)
                               gTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGc                                                                     >  1:3795762/1‑100 (MQ=255)
                               gTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTATTGAAGACGCTAAATTAGAAg                                                                      >  2:3474215/1‑99 (MQ=255)
                                ttGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCt                                                                    >  1:2762419/1‑100 (MQ=255)
                                 tGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTg                                                                   >  1:3350007/1‑100 (MQ=255)
                                  gCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTgg                                                                  >  1:1970141/1‑100 (MQ=255)
                                    aaaTAGATAGCTTTACGTAACAAACAGATAACATGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGc                                                                 <  2:3615716/99‑1 (MQ=255)
                                      tcttccgatctGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGtt                                                              <  1:1976234/90‑1 (MQ=255)
                                           ttGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcgcg                                                         <  1:677321/100‑1 (MQ=255)
                                            tGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                                        <  1:453970/100‑1 (MQ=255)
                                            tGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                                        <  1:2124746/100‑1 (MQ=255)
                                             gCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGa                                                       <  1:164457/100‑1 (MQ=255)
                                              cTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAg                                                      >  1:2103100/1‑100 (MQ=255)
                                              cTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAg                                                      >  1:1416295/1‑100 (MQ=255)
                                              cTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAg                                                      <  1:2502047/100‑1 (MQ=255)
                                                tGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTa                                                    >  2:3473861/1‑100 (MQ=255)
                                                 gACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTagatcggaaga                                                   >  2:1976223/1‑90 (MQ=255)
                                                 gACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGACAGCGGGCGCTGGCTATTGAAGACGCTAGAGTAGAAGGTGGCGGTACGCGGAGTa                                                    >  2:656831/1‑99 (MQ=255)
                                                  aCTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACg                                                  >  2:3580197/1‑100 (MQ=255)
                                                    ttAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTg                                                <  1:1554279/100‑1 (MQ=255)
                                                    ttAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTg                                                <  2:2474920/100‑1 (MQ=255)
                                                      aacaaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTggn                                              <  2:870078/100‑2 (MQ=255)
                                                      aacaaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTggg                                              >  1:2550699/1‑100 (MQ=255)
                                                      aacaaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTggg                                              >  1:1173404/1‑100 (MQ=255)
                                                      aacaaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTggg                                              <  2:2184419/100‑1 (MQ=255)
                                                      aacaaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTggg                                              <  1:1831478/100‑1 (MQ=255)
                                                      aacaaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTggg                                              <  2:1783540/100‑1 (MQ=255)
                                                      aacaaacaGAGAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTgg                                               <  2:3663760/99‑1 (MQ=255)
                                                        caaacaGATAACTTGACAGAAACGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGaa                                            >  2:2215780/1‑100 (MQ=255)
                                                          aacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACg                                          <  1:2872224/100‑1 (MQ=255)
                                                          aacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACg                                          <  2:1460103/100‑1 (MQ=255)
                                                          aaaaGATAACTTCACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCAGGCTCTCGAATACGCAAAATTATAAGCTGGCGTTACTCGGAGTACGTTGGAACg                                          >  2:1698780/4‑100 (MQ=255)
                                                           atctaTAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACgg                                         <  2:2758284/96‑1 (MQ=255)
                                                            caGATAACTTGACAGTAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGa                                        <  1:3755420/100‑1 (MQ=255)
                                                              gATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAat                                      >  1:148548/1‑100 (MQ=255)
                                                              gATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAat                                      >  1:1679231/1‑100 (MQ=255)
                                                              gATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAat                                      >  2:631753/1‑100 (MQ=255)
                                                               aTAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGagat                                     >  1:2758304/1‑97 (MQ=255)
                                                                tAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatca                                    <  1:2674272/100‑1 (MQ=255)
                                                                tAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatca                                    <  2:3033226/100‑1 (MQ=255)
                                                                tAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatca                                    <  2:1679221/100‑1 (MQ=255)
                                                                  aCTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcaca                                  >  2:2187914/1‑100 (MQ=255)
                                                                  aCTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcaca                                  <  2:1679631/100‑1 (MQ=255)
                                                                     tGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatca                               <  2:2885351/100‑1 (MQ=255)
                                                                        cAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgt                            >  1:1027307/1‑100 (MQ=255)
                                                                          gAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtac                          >  2:3503696/1‑100 (MQ=255)
                                                                              aGATAAAAAAAGAGCCAGCGGTCGCTGTCTCTTGAATACGCTAAATTAGAATCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGttt                      <  2:1882823/100‑1 (MQ=255)
                                                                               gATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGtttt                     <  1:3010354/100‑1 (MQ=255)
                                                                                aTAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTg                    >  2:2454302/1‑100 (MQ=255)
                                                                                aTAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTg                    <  2:1416285/100‑1 (MQ=255)
                                                                                 tAAACAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTACATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                   <  1:3580224/100‑1 (MQ=255)
                                                                                 tAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                   <  2:2750882/100‑1 (MQ=255)
                                                                                  tcttcaaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCa                  >  1:3034573/6‑100 (MQ=38)
                                                                                   cttcaaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAc                 >  2:1473367/5‑100 (MQ=39)
                                                                                   aaaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAc                 >  1:1843703/1‑100 (MQ=255)
                                                                                    ttcaaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCa                  >  2:861778/4‑98 (MQ=255)
                                                                                    aaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACg                >  1:1220425/1‑100 (MQ=255)
                                                                                     tcaaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTCCGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAc                 >  1:1992073/3‑98 (MQ=39)
                                                                                      caaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGcc              >  1:918340/2‑100 (MQ=255)
                                                                                      aaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGcc              <  1:3404642/100‑1 (MQ=255)
                                                                                      aaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGcc              <  1:631760/100‑1 (MQ=255)
                                                                                       aaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCa             >  2:291538/1‑100 (MQ=255)
                                                                                        agagCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAg            >  2:3385449/1‑100 (MQ=255)
                                                                                        agagCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAg            >  1:1520735/1‑100 (MQ=255)
                                                                                        agagCCAGCGGGCGCTGGCTCTTGAAGACGAGAAATCAGAAGCTGGCGTTAGGCGGAGCACGAGGGAACGGAATCACATCACGTACGTTGTGc                   >  2:1951948/1‑93 (MQ=255)
                                                                                          agCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGtt          >  2:1522081/1‑100 (MQ=255)
                                                                                             cAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGGTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACACCAGTTACg       >  2:1399198/1‑100 (MQ=255)
                                                                                             cAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACg       <  2:1173395/100‑1 (MQ=255)
                                                                                               gCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTa     <  1:3650955/100‑1 (MQ=255)
                                                                                                cGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGGACGGAATCACATCACGTACGTTTTGGACGCCAGTTACGt      >  1:1511999/1‑98 (MQ=255)
                                                                                                  ggCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCCCGTACGTTTTGCACGCCAGTTACGTAAGc  >  2:2012853/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: