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,8980AG99.3% 1022.7 / ‑3.3 298intergenic (‑16/+63)yddE/narVhypothetical protein/cryptic nitrate reductase 2 gamma subunit
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (148/148);  total (150/148)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

Only 100 of 304 total aligned reads displayed.

GAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATA  >  NC_000913/1533808‑1533987
                                                                                          |                                                                                         
 aaaaCAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTg                                                                                        <  1:3382518/93‑1 (MQ=255)
  aaaCAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTgg                                                                                       >  1:7083385/1‑93 (MQ=255)
   aacaacACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGga                                                                                      >  2:3156759/1‑93 (MQ=255)
   aacaacACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGga                                                                                      <  2:1632552/93‑1 (MQ=255)
   aacaacACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGga                                                                                      >  2:1196886/1‑93 (MQ=255)
   aacaacACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGga                                                                                      <  1:6951809/93‑1 (MQ=255)
    acaacaCCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTgggg                                                                                     >  2:5373309/1‑91 (MQ=255)
    acaacaCCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTggcc                                                                                     >  1:2833168/1‑91 (MQ=255)
    acaacaCCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGgag                                                                                     >  2:2343973/1‑93 (MQ=255)
    acaacaCCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGgag                                                                                     >  1:3941802/1‑93 (MQ=255)
     caacaCCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGgaga                                                                                    >  2:1544086/1‑93 (MQ=255)
      aacaCCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGaa                                                                                   <  1:4319403/93‑1 (MQ=255)
       acacCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAAc                                                                                  >  1:30421/1‑93 (MQ=255)
         acCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTa                                                                                <  1:1232009/93‑1 (MQ=255)
          ccGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTAt                                                                               <  1:384349/93‑1 (MQ=255)
            ggCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCa                                                                             <  2:6153183/93‑1 (MQ=255)
             gCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCAt                                                                            <  2:2784967/93‑1 (MQ=255)
             gCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCAt                                                                            >  2:8069443/1‑93 (MQ=255)
                gAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGc                                                                         <  2:3000236/93‑1 (MQ=255)
                gAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGc                                                                         >  1:4624498/1‑93 (MQ=255)
                gAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGc                                                                         <  2:3580168/93‑1 (MQ=255)
                gAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGc                                                                         >  1:5242381/1‑93 (MQ=255)
                 aaTTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGcc                                                                        >  2:4185704/1‑93 (MQ=255)
                  aTTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCt                                                                       <  1:7493001/93‑1 (MQ=255)
                  aTTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCt                                                                       <  1:4044104/93‑1 (MQ=255)
                  aTTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCt                                                                       <  1:32951/93‑1 (MQ=255)
                   ttGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTg                                                                      <  1:5128774/93‑1 (MQ=255)
                    tGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGc                                                                     <  2:7769460/93‑1 (MQ=255)
                     gCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCa                                                                    <  2:8504526/93‑1 (MQ=255)
                     gCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCa                                                                    <  1:6072926/93‑1 (MQ=255)
                     gCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCa                                                                    >  2:1879599/1‑93 (MQ=255)
                      ccGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCaa                                                                   <  2:3286467/93‑1 (MQ=255)
                        gcgAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGt                                                                 <  1:6795287/93‑1 (MQ=255)
                        ccgAAACGGTTGGAAGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTg                                                                <  1:8449670/92‑1 (MQ=255)
                        cagAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGt                                                                 <  1:6144689/91‑1 (MQ=255)
                         cgAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTg                                                                <  1:4959913/93‑1 (MQ=255)
                         cgAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTg                                                                <  2:5621215/93‑1 (MQ=255)
                         cgAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTg                                                                <  2:7253041/93‑1 (MQ=255)
                           aaaCGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGc                                                              <  2:4866597/93‑1 (MQ=255)
                           aaaCGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGc                                                              >  2:2875240/1‑93 (MQ=255)
                           aaaCGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGc                                                              >  2:3442189/1‑93 (MQ=255)
                           aaaCGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGc                                                              <  2:7889117/93‑1 (MQ=255)
                            aaCGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGcc                                                             >  2:5324078/1‑93 (MQ=255)
                            aaCGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGcc                                                             >  2:6385065/1‑93 (MQ=255)
                                  tgtgAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGagag                                                       <  1:6854594/93‑1 (MQ=255)
                                  tgtgAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGagag                                                       <  1:8613495/93‑1 (MQ=255)
                                  tgtgAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGagag                                                       <  1:4679009/93‑1 (MQ=255)
                                    tgAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCg                                                     <  2:2452092/93‑1 (MQ=255)
                                     gAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGGGCGa                                                    >  1:6399511/1‑93 (MQ=255)
                                     gAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGa                                                    <  2:5098160/93‑1 (MQ=255)
                                      aGGTAAAAGCATCGGCGTTGTACGCCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAATTATCATAGCCTGCAAGTGGCCGGAGAGCGaa                                                   <  2:2696967/93‑1 (MQ=255)
                                       ggTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAg                                                  <  1:5515440/93‑1 (MQ=255)
                                        gTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAgg                                                 <  1:801457/93‑1 (MQ=255)
                                         tAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAggg                                                >  1:901010/1‑93 (MQ=255)
                                          aaaaGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGc                                               <  1:8254951/93‑1 (MQ=255)
                                          aaaaGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGc                                               <  2:5787180/93‑1 (MQ=255)
                                          aaaaGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGGGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGc                                               <  2:1785686/93‑1 (MQ=255)
                                            aaGCATCGACGTGGTACACCTGCGGTTTCATTAACGTACTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTa                                             <  1:2431945/93‑1 (MQ=255)
                                             aGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTAt                                            <  2:2246185/93‑1 (MQ=255)
                                               cATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATcc                                          <  1:6667712/93‑1 (MQ=255)
                                               cATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATcc                                          <  2:7471082/93‑1 (MQ=255)
                                               cATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATcc                                          <  2:6045006/93‑1 (MQ=255)
                                                aTCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCg                                         >  1:1954162/1‑93 (MQ=255)
                                                   gACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCCGcc                                      >  2:6502704/1‑93 (MQ=255)
                                                    aCGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGTAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCa                                     >  1:5430125/1‑93 (MQ=255)
                                                       tGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAggg                                  <  1:7667952/93‑1 (MQ=255)
                                                        ggTACCCCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGt                                 <  1:5672349/93‑1 (MQ=255)
                                                        ggTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAgggg                                 >  1:8098418/1‑92 (MQ=255)
                                                         gTACACCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTg                                <  2:6701636/93‑1 (MQ=255)
                                                           acacCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGGTATCCGGCCAGGGTGaa                              >  2:4327180/1‑93 (MQ=255)
                                                            cacCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGTACTATCATAGCCTGCAAGTGGCCGGAGGGCGAAGGGCTTTCCGGCCAGGGAGaaa                             >  2:8578275/1‑93 (MQ=255)
                                                            cacCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGaaa                             <  2:7054407/93‑1 (MQ=255)
                                                             acCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAAt                            >  1:345737/1‑93 (MQ=255)
                                                             acCTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAAt                            >  2:8446928/1‑93 (MQ=255)
                                                              ccTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGTGAGCGAAGGGCTATCCGGCCAGGGTGAAAtt                           >  1:3187368/1‑93 (MQ=255)
                                                              ccTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGCGCGAAGGGGTATCCGGCCAGGGGGAAAtt                           >  1:1694589/1‑93 (MQ=255)
                                                              ccTGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAAtt                           >  2:2618771/1‑93 (MQ=255)
                                                                tGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTAt                         >  1:6311691/1‑93 (MQ=255)
                                                                tGCGGTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTAt                         <  1:7578136/93‑1 (MQ=255)
                                                                    gTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCTGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATcgcc                     <  2:3073035/93‑1 (MQ=255)
                                                                    gTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATcgcc                     >  1:2366617/1‑93 (MQ=255)
                                                                    gTTTCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATcgcc                     >  2:3654530/1‑93 (MQ=255)
                                                                       tCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGGCAGGGTGAAATTATCGCcgcg                  >  2:8284839/1‑93 (MQ=255)
                                                                       tCATTAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCcgcg                  >  2:889224/1‑93 (MQ=255)
                                                                         cttAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGaa                <  1:2072347/92‑1 (MQ=255)
                                                                          ttAACGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAAc               <  2:3889391/93‑1 (MQ=255)
                                                                             aCGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGca            >  1:1380329/1‑93 (MQ=255)
                                                                              cGTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGcac           <  1:6456560/93‑1 (MQ=255)
                                                                               gTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGcaca          >  2:6241557/1‑93 (MQ=255)
                                                                               gTTCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGcaca          >  1:4460476/1‑93 (MQ=255)
                                                                                ttCTCCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACaa         <  1:4773457/93‑1 (MQ=255)
                                                                                 tctcCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGGCGCGCACCCACAAt        >  2:1452498/1‑93 (MQ=255)
                                                                                 tctcCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAAt        <  1:7287157/93‑1 (MQ=255)
                                                                                   tcCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAAttt      >  2:7199771/1‑93 (MQ=255)
                                                                                   tcCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAAttt      <  2:4411391/93‑1 (MQ=255)
                                                                                   tcCTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAAttt      >  2:4968193/1‑93 (MQ=255)
                                                                                     cTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGa    <  1:2631847/93‑1 (MQ=255)
                                                                                     cTGTGGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGa    >  2:4620576/1‑93 (MQ=255)
                                                                                      tgtgGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGat   <  1:1050280/93‑1 (MQ=255)
                                                                                       gtgGCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGata  <  1:738583/93‑1 (MQ=255)
                                                                                          |                                                                                         
GAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATA  >  NC_000913/1533808‑1533987

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