Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,533,898 A→G intergenic (‑16/+63) yddE ← / ← narV hypothetical protein/cryptic nitrate reductase 2 gamma subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,533,8980AG87.3% 592.9 / 71.0 237intergenic (‑16/+63)yddE/narVhypothetical protein/cryptic nitrate reductase 2 gamma subunit
Reads supporting (aligned to +/- strand):  ref base A (14/16);  new base G (93/114);  total (107/130)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

Only 100 of 238 total aligned reads displayed.

GGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATC  >  NC_000913/1533806‑1533989
                                                                                            |                                                                                           
   aaaaCAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTg                                                                                          >  1:366023/1‑93 (MQ=255)
         acacCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAAc                                                                                    >  2:2110880/1‑93 (MQ=255)
         acacCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAAc                                                                                    <  1:288830/93‑1 (MQ=255)
           acCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTa                                                                                  <  1:141972/93‑1 (MQ=255)
           acCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTa                                                                                  <  1:1939582/93‑1 (MQ=255)
            ccGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTAt                                                                                 <  2:7056871/93‑1 (MQ=255)
             cGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATc                                                                                <  1:582431/93‑1 (MQ=255)
             cGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATc                                                                                <  2:1424517/93‑1 (MQ=255)
              ggCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCa                                                                               <  1:818121/93‑1 (MQ=255)
              ggCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCa                                                                               <  2:2420030/93‑1 (MQ=255)
               gCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCAt                                                                              <  1:102481/93‑1 (MQ=255)
                cAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATa                                                                             <  1:5154682/93‑1 (MQ=255)
                cAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATa                                                                             <  2:38994/93‑1 (MQ=255)
                cAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATa                                                                             <  1:3997719/93‑1 (MQ=255)
                 aGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAg                                                                            >  2:3882164/1‑93 (MQ=255)
                  gAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGc                                                                           >  1:5673091/1‑93 (MQ=255)
                   aaTTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGcc                                                                          >  1:5367503/1‑93 (MQ=255)
                     ttGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTg                                                                        >  2:6765352/1‑93 (MQ=255)
                         cgcgAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAg                                                                    >  1:6743680/1‑93 (MQ=255)
                          gcgAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGt                                                                   >  2:3129511/1‑93 (MQ=255)
                          gcgAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGt                                                                   >  1:543310/1‑93 (MQ=255)
                           cgAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTg                                                                  <  1:5552136/93‑1 (MQ=255)
                           cgAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTg                                                                  >  2:6071544/1‑93 (MQ=255)
                           cgAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTg                                                                  >  1:6893363/1‑93 (MQ=255)
                              aaCGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGcc                                                               <  2:4376995/93‑1 (MQ=255)
                              aaCGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGcc                                                               <  2:5764961/93‑1 (MQ=255)
                                 ggTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGga                                                            <  1:1008129/93‑1 (MQ=255)
                                  gTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGgag                                                           <  2:6483455/93‑1 (MQ=255)
                                  gTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGgag                                                           >  1:6102760/1‑93 (MQ=255)
                                  gTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGgag                                                           >  2:4039063/1‑93 (MQ=255)
                                   ttGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGgaga                                                          >  1:5235878/1‑93 (MQ=255)
                                   ttGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCACAGCCTGCAAGTGGCCGgaga                                                          <  2:6712971/93‑1 (MQ=255)
                                       gAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGa                                                      <  2:4582405/93‑1 (MQ=255)
                                         ggTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAg                                                    <  1:3233449/93‑1 (MQ=255)
                                           tAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAggg                                                  >  2:2027449/1‑93 (MQ=255)
                                           tAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAggg                                                  <  1:1107536/93‑1 (MQ=255)
                                            aaaaGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGc                                                 <  2:277218/93‑1 (MQ=255)
                                            aaaaGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGc                                                 <  1:957452/93‑1 (MQ=255)
                                            aaaaGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGc                                                 <  1:4135189/93‑1 (MQ=255)
                                            aaaaGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGc                                                 <  1:4097337/93‑1 (MQ=255)
                                            aaaaGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGc                                                 <  1:3602538/93‑1 (MQ=255)
                                            aaaaGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGc                                                 <  1:6490921/93‑1 (MQ=255)
                                              aaGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTa                                               >  1:5193611/1‑93 (MQ=255)
                                              aaGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTa                                               >  1:6254372/1‑93 (MQ=255)
                                              aaGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTa                                               >  2:6861727/1‑93 (MQ=255)
                                                gCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATc                                             <  2:1285103/93‑1 (MQ=255)
                                                   tCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCgg                                          >  1:2281333/1‑93 (MQ=255)
                                                   tCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCgg                                          >  2:4691518/1‑93 (MQ=255)
                                                   tCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAAATATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCgg                                          >  1:4623141/1‑93 (MQ=255)
                                                    cGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGc                                         <  1:7410393/93‑1 (MQ=255)
                                                     caCGTGGTACCCCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGGGGCCGGAGAGCGAAGGGCTATCCGGcc                                        <  1:2990933/92‑1 (MQ=255)
                                                      aCGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGGGGCCGGAGAGCGAAGGGCTATCCGGCCa                                       <  1:1807376/93‑1 (MQ=255)
                                                       cGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAg                                      <  2:5338430/93‑1 (MQ=255)
                                                        gTGGTCCACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGGGGCCGGAGAGCGAAGGGCTATCCGGCCAgg                                     <  2:6425954/93‑1 (MQ=255)
                                                         tGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCCggg                                    >  1:3498442/1‑93 (MQ=255)
                                                         tGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAggg                                    >  1:6997469/1‑93 (MQ=255)
                                                          ggTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGt                                   <  2:2928029/93‑1 (MQ=255)
                                                          ggTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGt                                   >  2:6111231/1‑93 (MQ=255)
                                                          ggTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGt                                   <  2:879736/93‑1 (MQ=255)
                                                          ggTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGt                                   <  1:2478647/93‑1 (MQ=255)
                                                          ggTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGt                                   <  1:5344790/93‑1 (MQ=255)
                                                             acacCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGaa                                >  2:1571532/1‑93 (MQ=255)
                                                              cacCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGaaa                               >  2:3419311/1‑93 (MQ=255)
                                                              cacCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGaaa                               <  2:7292834/93‑1 (MQ=255)
                                                              cacCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGaaa                               <  2:852337/93‑1 (MQ=255)
                                                                ccTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAAtt                             >  2:3805145/1‑93 (MQ=255)
                                                                 cTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTa                            <  1:7265575/93‑1 (MQ=255)
                                                                 cTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTa                            <  2:6743084/93‑1 (MQ=255)
                                                                  tGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTAt                           >  1:7324199/1‑93 (MQ=255)
                                                                  tGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTAt                           >  1:5494561/1‑93 (MQ=255)
                                                                  ggCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTAt                           <  1:2423703/92‑1 (MQ=255)
                                                                   gCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATc                          <  1:2580243/93‑1 (MQ=255)
                                                                   gCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATc                          >  2:3229885/1‑93 (MQ=255)
                                                                   gCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATc                          <  1:5968903/93‑1 (MQ=255)
                                                                   gCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATc                          <  1:5388807/93‑1 (MQ=255)
                                                                     ggTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATcgc                        >  1:1366384/1‑93 (MQ=255)
                                                                      gTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATcgcc                       <  1:2073626/93‑1 (MQ=255)
                                                                      gTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATcgcc                       <  1:2096982/93‑1 (MQ=255)
                                                                       tttCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATcgccg                      >  1:4374014/1‑93 (MQ=255)
                                                                       tttCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATcgccg                      <  1:314170/93‑1 (MQ=255)
                                                                        ttCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATcgccgc                     <  2:6638317/93‑1 (MQ=255)
                                                                        ttCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATcgccgc                     >  1:2833991/1‑93 (MQ=255)
                                                                           aTTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGATAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGaa                  >  2:1695431/1‑93 (MQ=255)
                                                                           aTTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGGAATTATCGCCGCGaa                  >  1:7201390/1‑93 (MQ=255)
                                                                           aTTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGaa                  <  1:3033897/93‑1 (MQ=255)
                                                                           aTTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGaa                  <  2:3561324/93‑1 (MQ=255)
                                                                            ttAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAAc                 <  2:1861206/93‑1 (MQ=255)
                                                                            ttAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAAc                 <  2:5034607/93‑1 (MQ=255)
                                                                            ttAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGGGAAATTATCGCCGCGacc                 >  2:4879705/1‑91 (MQ=255)
                                                                              aaCGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGc               >  2:7388585/1‑93 (MQ=255)
                                                                              aaCGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCAGCCAGGGTGAAATTATCGCCGCGAACGc               <  1:752670/93‑1 (MQ=255)
                                                                                cGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGcac             >  2:2261816/1‑93 (MQ=255)
                                                                                 gTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGcaca            >  2:4416354/1‑93 (MQ=255)
                                                                                  ttCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACaa           <  1:5000859/93‑1 (MQ=255)
                                                                                  ttCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACaa           <  2:6244252/93‑1 (MQ=255)
                                                                                  ttCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACaa           <  1:156877/93‑1 (MQ=255)
                                                                                        tgtgGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGat     >  2:2632615/1‑93 (MQ=255)
                                                                                         gtgGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGata    <  1:6311209/93‑1 (MQ=255)
                                                                                         gtgGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGata    <  1:269013/93‑1 (MQ=255)
                                                                                         gtgGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGata    >  2:1147104/1‑93 (MQ=255)
                                                                                            |                                                                                           
GGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATC  >  NC_000913/1533806‑1533989

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