Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 369,076 A→G K156E (AAA→GAA)  mhpA → 3‑(3‑hydroxyphenyl)propionate hydroxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913369,0760AG100.0% 273.8 / NA 76K156E (AAA→GAA) mhpA3‑(3‑hydroxyphenyl)propionate hydroxylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (38/38);  total (38/38)

GTAACGCCTTTATTCAGCCGCAGGTCGATGCGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGATTGTGGTA  >  NC_000913/368942‑369210
                                                                                                                                      |                                                                                                                                      
gTAACGCCTTTATTCAGCCGCAGGTCGATGCGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAAc                                                                                                                                    <  1:48569/139‑1 (MQ=255)
gTAACGCCTTTATTCAGCCGCAGGTCGATGCGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAAc                                                                                                                                    <  1:1560822/139‑1 (MQ=255)
               aGCCGCAGGTCGATGCGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCg                                                                                                                     >  1:300138/1‑139 (MQ=255)
                 ccGCAGGTCGATGCGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAAc                                                                                                                                    <  1:847672/122‑1 (MQ=255)
                 ccGCAGGTCGATGCGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAAc                                                                                                                                    >  2:847672/1‑122 (MQ=255)
                 ccGCAGGTCGATGCGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCggg                                                                                                                   >  2:467405/1‑139 (MQ=255)
                      ggTCGATGCGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATa                                                                                                              >  2:611824/1‑139 (MQ=255)
                           aTGCGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCaa                                                                                                         >  1:1397682/1‑139 (MQ=255)
                           aTGCGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCAAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCaa                                                                                                         >  2:921652/1‑139 (MQ=255)
                                  gATGCTGGAAGGGGTGTCGCATTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAg                                                                                                  <  2:1522419/139‑1 (MQ=255)
                                   aTGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGt                                                                                                 <  1:467405/139‑1 (MQ=255)
                                         gAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAAt                                                                                                               <  1:1942963/119‑1 (MQ=255)
                                         gAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAAt                                                                                                               >  2:1942963/1‑119 (MQ=255)
                                          aaGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGt                                                                                          >  2:262454/1‑139 (MQ=255)
                                              ggTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGcc                                                                                      >  1:389013/1‑139 (MQ=255)
                                              ggTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGcc                                                                                      >  1:549912/1‑139 (MQ=255)
                                               gtgtCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAggg                                                                                                                          >  2:1770304/1‑102 (MQ=255)
                                               gtgtCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAggg                                                                                                                          <  1:1770304/102‑1 (MQ=255)
                                                  tCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCtgtg                                                                                  >  1:1639150/1‑139 (MQ=255)
                                                  tCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCtgtg                                                                                  <  2:1397682/139‑1 (MQ=255)
                                                      gTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACgg                                                                              <  1:262454/139‑1 (MQ=255)
                                                        tttCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAggg                                                                                                                          >  1:1476651/1‑93 (MQ=255)
                                                        tttCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAggg                                                                                                                          <  2:1476651/93‑1 (MQ=255)
                                                                    cgcTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATTACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCAGGAAATAGTCAAAGCCCAGTGGCTgg                                                                                           <  1:2039929/112‑1 (MQ=255)
                                                                    cgcTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACgg                                                                                                                                  <  2:296822/73‑1 (MQ=255)
                                                                    cgcTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACgg                                                                                                                                  >  1:296822/1‑73 (MQ=255)
                                                                    cgcTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGaaa                                                                                                                >  2:1705113/1‑91 (MQ=255)
                                                                    cgcTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGaaa                                                                                                                <  1:1705113/91‑1 (MQ=255)
                                                                    cgcTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTgg                                                                                           >  2:2039929/1‑112 (MQ=255)
                                                                    cgcTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTg                                                                >  1:1577273/1‑139 (MQ=255)
                                                                    cgcTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTg                                                                >  1:14322/1‑139 (MQ=255)
                                                                    cgcTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTg                                                                >  2:1220169/1‑139 (MQ=255)
                                                                        gctTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCg                                                            <  2:536621/139‑1 (MQ=255)
                                                                         ctTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGt                                                           >  2:2013155/1‑139 (MQ=255)
                                                                           tGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCg                                                         <  1:2026307/139‑1 (MQ=255)
                                                                           tGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCg                                                         <  2:549912/139‑1 (MQ=255)
                                                                            gTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGc                                                        >  1:965325/1‑139 (MQ=255)
                                                                              tttCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCAc                                                      <  2:14322/139‑1 (MQ=255)
                                                                                         cTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCa                                                                                                   >  1:1942407/1‑83 (MQ=255)
                                                                                         cTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCa                                                                                                   <  2:1942407/83‑1 (MQ=255)
                                                                                              ggCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGttt                                      <  1:1220169/139‑1 (MQ=255)
                                                                                               gCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTg                                     <  1:611824/139‑1 (MQ=255)
                                                                                                  ttCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCa                                                                                                          <  1:1046560/67‑1 (MQ=255)
                                                                                                  ttCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCa                                                                                                          >  2:1046560/1‑67 (MQ=255)
                                                                                                  ttCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAg                                  >  2:1917477/1‑139 (MQ=255)
                                                                                                         agCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAAc                           >  1:1119918/1‑139 (MQ=255)
                                                                                                             aGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTgcg                       <  1:809399/139‑1 (MQ=255)
                                                                                                             aGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTgcg                       <  2:392531/139‑1 (MQ=255)
                                                                                                              gATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTgcgc                      >  2:1227647/1‑139 (MQ=255)
                                                                                                                tGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCa                    <  1:2013155/139‑1 (MQ=255)
                                                                                                                 gACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTgcgc                      >  2:335663/1‑136 (MQ=255)
                                                                                                                 gACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTgcgc                      <  1:335663/136‑1 (MQ=255)
                                                                                                                 gACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGCAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCaa                   >  1:74198/1‑139 (MQ=255)
                                                                                                                   cGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAAt                 >  1:1141789/1‑139 (MQ=255)
                                                                                                                        tGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCa                                                                                                   <  1:57787/52‑1 (MQ=255)
                                                                                                                        tGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCa                                                                                                   >  2:57787/1‑52 (MQ=255)
                                                                                                                        tGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTg            <  2:74198/139‑1 (MQ=255)
                                                                                                                        tGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTg            >  2:24438/1‑139 (MQ=255)
                                                                                                                         gACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGttt                                      <  1:1400708/112‑1 (MQ=255)
                                                                                                                         gACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGttt                                      >  2:1400708/1‑112 (MQ=255)
                                                                                                                         gACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTgg           >  2:940328/1‑139 (MQ=255)
                                                                                                                         gACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTgg           >  1:1478392/1‑139 (MQ=255)
                                                                                                                          aCCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCg                                                         >  2:1673518/1‑92 (MQ=255)
                                                                                                                          aCCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCg                                                         <  1:1673518/92‑1 (MQ=255)
                                                                                                                          aCCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGa          <  2:1706135/139‑1 (MQ=255)
                                                                                                                          aCCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGa          <  2:300138/139‑1 (MQ=255)
                                                                                                                          aCCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGa          <  2:1478392/139‑1 (MQ=255)
                                                                                                                          aCCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGa          <  2:133818/139‑1 (MQ=255)
                                                                                                                          aCCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGa          <  2:1119918/139‑1 (MQ=255)
                                                                                                                          aCCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGa          <  1:84548/139‑1 (MQ=255)
                                                                                                                          aCCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGa          >  1:133818/1‑139 (MQ=255)
                                                                                                                          aCCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGa          <  1:1077200/139‑1 (MQ=255)
                                                                                                                            cTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGAtt        <  1:1227647/139‑1 (MQ=255)
                                                                                                                              tGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTc                                                          >  2:1524176/1‑87 (MQ=255)
                                                                                                                              tGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTc                                                          <  1:1524176/87‑1 (MQ=255)
                                                                                                                                  ccTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGATTGTGGTa  >  2:660484/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                      
GTAACGCCTTTATTCAGCCGCAGGTCGATGCGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGATTGTGGTA  >  NC_000913/368942‑369210

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