Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 987,560 G→C 16.1% 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,5600GC16.1% 346.8 / 55.9 130intergenic (‑578/+25)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (49/60);  new base C (7/14);  total (56/74)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

Only 100 of 136 total aligned reads displayed.

ATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGC  >  NC_000913/987461‑987659
                                                                                                   |                                                                                                   
atatTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATCGATTGCTTGCCTGTACAAACAGATAACTTGACAGAAAAGATAGAAAAAGAGCCAGCgg                                                                                                     >  2:1481550/1‑100 (MQ=255)
 tatTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCggg                                                                                                    >  1:745675/1‑100 (MQ=255)
 tatTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCggg                                                                                                    <  2:1926495/100‑1 (MQ=255)
   ttCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgcg                                                                                                  <  2:1732218/100‑1 (MQ=255)
        tAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGc                                                                                             <  1:1941839/100‑1 (MQ=255)
        tAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGc                                                                                             >  1:1376194/1‑100 (MQ=255)
         aaCAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGct                                                                                            <  2:3578890/100‑1 (MQ=255)
          aCAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGctc                                                                                           <  1:3000246/100‑1 (MQ=255)
          aCAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGctc                                                                                           <  1:497239/100‑1 (MQ=255)
           cAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGctct                                                                                          <  1:777139/100‑1 (MQ=255)
           cAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGctct                                                                                          <  2:1429701/100‑1 (MQ=255)
            aGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACTGATAACTTGACAGAAAAGATAATAAAAGAGCCAGCGGGCGCTGGCTCtt                                                                                         <  1:3789684/100‑1 (MQ=255)
            aGGATGACCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGctct                                                                                          <  1:687037/99‑1 (MQ=255)
             ggATGATCGTTATGGCGGGGTTGCAAATAGATTGGTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTg                                                                                        <  2:1173247/100‑1 (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)
                  atCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCttttttta                                                                                   <  1:838547/100‑7 (MQ=38)
                  atCGTTATGGCGGGGTTGCAAAAAGATTGCATGACTTAACAAGCAGATAACATGAAACAAAAGATATAAACGGAGCCCACGGGCGCTGTctctctacgac                                                                                   >  2:3398273/1‑93 (MQ=255)
                     gTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCt                                                                                <  2:3206132/100‑1 (MQ=255)
                     gTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCt                                                                                <  1:3264352/100‑1 (MQ=255)
                     gTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGACAGCGGGCGCTGGCTCTTGAAGACGCt                                                                                >  2:2998852/1‑100 (MQ=255)
                         tGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAAt                                                                            <  1:233387/100‑1 (MQ=255)
                         tGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAAt                                                                            <  1:2561257/100‑1 (MQ=255)
                          ggCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAAtt                                                                           >  2:1422109/1‑100 (MQ=255)
                           gCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTa                                                                          >  1:86716/1‑100 (MQ=255)
                           gCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTa                                                                          >  1:2512344/1‑100 (MQ=255)
                            cGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAg                                                                         <  2:1654916/100‑1 (MQ=255)
                              gggTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGaa                                                                       <  1:3292878/100‑1 (MQ=255)
                               ggTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAg                                                                      <  2:2015341/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)
                                             tGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                                        <  1:2124746/100‑1 (MQ=255)
                                             tGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                                        <  1:453970/100‑1 (MQ=255)
                                              gCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGa                                                       <  1:164457/100‑1 (MQ=255)
                                               cTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAg                                                      >  1:1416295/1‑100 (MQ=255)
                                               cTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAg                                                      >  1:2103100/1‑100 (MQ=255)
                                               cTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAg                                                      <  1:2502047/100‑1 (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                                              <  2:2184419/100‑1 (MQ=255)
                                                       aacaaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTggg                                              >  1:1173404/1‑100 (MQ=255)
                                                       aacaaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTggg                                              <  1:1831478/100‑1 (MQ=255)
                                                       aacaaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTggg                                              >  1:2550699/1‑100 (MQ=255)
                                                       aacaaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTggg                                              <  1:3689946/100‑1 (MQ=255)
                                                       aacaaacaGAGAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTgg                                               <  2:3663760/99‑1 (MQ=255)
                                                          aaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAAc                                           >  1:3441842/1‑100 (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)
                                                                 tAACTTGCCAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatca                                    >  1:2003146/1‑100 (MQ=255)
                                                                 tAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatca                                    <  2:1679221/100‑1 (MQ=255)
                                                                   aCTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcaca                                  >  2:2187914/1‑100 (MQ=255)
                                                                      tGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatca                               <  2:2885351/100‑1 (MQ=255)
                                                                         cAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgt                            >  1:1027307/1‑100 (MQ=255)
                                                                          aGAAAAGATAAAAAAAGAGCCAGCGGTCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCTGAGTACGTGGGAACGGAATCACATCacgta                           >  2:148527/1‑100 (MQ=255)
                                                                           gAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtac                          >  2:3503696/1‑100 (MQ=255)
                                                                            aaaaGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtac                          <  1:2708036/99‑1 (MQ=255)
                                                                               aGATAAAAAAAGAGCCAGCGGTCGCTGTCTCTTGAATACGCTAAATTAGAATCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGttt                      <  2:1882823/100‑1 (MQ=255)
                                                                                gATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGtttt                     <  1:3010354/100‑1 (MQ=255)
                                                                                 aTAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTg                    <  2:1416285/100‑1 (MQ=255)
                                                                                 aTAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTg                    >  2:2454302/1‑100 (MQ=255)
                                                                                  tAAACAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTACATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                   <  1:3580224/100‑1 (MQ=255)
                                                                                  tAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                   <  2:2750882/100‑1 (MQ=255)
                                                                                   tcttcaaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCa                  >  1:3034573/6‑100 (MQ=38)
                                                                                   aaaaaaaGGGGGAGCGGGCGCTGGCTCTTGAAGACGCTAAATGAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGCAATCACATCACGTACGTTTTGc                   <  2:2209399/99‑1 (MQ=255)
                                                                                    cttcaaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAc                 >  2:1473367/5‑100 (MQ=39)
                                                                                    aaaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAc                 >  1:1843703/1‑100 (MQ=255)
                                                                                     ttcaaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCa                  >  2:861778/4‑98 (MQ=255)
                                                                                     ttcaaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACg                >  1:489954/4‑100 (MQ=255)
                                                                                     aaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACg                >  1:1220425/1‑100 (MQ=255)
                                                                                       caaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGcc              >  1:918340/2‑100 (MQ=255)
                                                                                       aaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGcc              <  1:631760/100‑1 (MQ=255)
                                                                                        aaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCa             >  2:291538/1‑100 (MQ=255)
                                                                                         agagCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAg            >  1:1520735/1‑100 (MQ=255)
                                                                                         agagCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAg            >  2:3385449/1‑100 (MQ=255)
                                                                                           agCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGtt          >  2:1522081/1‑100 (MQ=255)
                                                                                              cAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACg       <  2:1173395/100‑1 (MQ=255)
                                                                                                 cGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGGACGGAATCACATCACGTACGTTTTGGACGCCAGTTACGt      >  1:1511999/1‑98 (MQ=255)
                                                                                                   ggCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCCCGTACGTTTTGCACGCCAGTTACGTAAGc  >  2:2012853/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGC  >  NC_000913/987461‑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: