Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 167,743 C→T T87I (ACC→ATC)  fhuA → ferrichrome outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913167,7430CT100.0% 248.9 / NA 76T87I (ACC→ATC) fhuAferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/44);  total (32/44)

TCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGA  >  NC_000913/167608‑167875
                                                                                                                                       |                                                                                                                                    
tCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGAtcgc                                                                                                                                   <  1:1508572/139‑1 (MQ=255)
          cGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATg                                                                                                                         >  1:55472/1‑139 (MQ=255)
          cGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATg                                                                                                                         >  2:273235/1‑139 (MQ=255)
            cTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGgc                                                                                                                       >  2:1539224/1‑139 (MQ=255)
             tGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGgcg                                                                                                                      <  1:629582/139‑1 (MQ=255)
                    gcgcCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATc                                                                                                               >  1:904481/1‑139 (MQ=255)
                          cAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGa                                                                                                         >  2:1494296/1‑139 (MQ=255)
                            aGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAg                                                                                                       <  1:683554/139‑1 (MQ=255)
                                      tGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAaag                                                                                             <  2:77046/139‑1 (MQ=255)
                                      tGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAaag                                                                                             <  2:32088/139‑1 (MQ=255)
                                       gggggCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAaaga                                                                                            <  2:904481/139‑1 (MQ=255)
                                       gggggCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAaaga                                                                                            <  2:1325603/139‑1 (MQ=255)
                                         gggCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAaagaag                                                                                          <  1:1539224/139‑1 (MQ=255)
                                         gggCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAaagaag                                                                                          >  1:1539176/1‑139 (MQ=255)
                                         gggCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAaagaag                                                                                          <  1:1143729/139‑1 (MQ=255)
                                         gggCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTC‑‑‑TGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAgcgc                                                                                       >  2:1568653/1‑139 (MQ=255)
                                                  gcAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCt                                                                                 >  2:976620/1‑139 (MQ=255)
                                                     aCTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTaca                                                                              <  1:892537/139‑1 (MQ=255)
                                                       tATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACg                                                                            <  1:8719/139‑1 (MQ=255)
                                                         ttGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGcc                                                                          <  2:667770/139‑1 (MQ=255)
                                                          tGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCg                                                                         <  1:1494296/139‑1 (MQ=255)
                                                             ggcgCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgt                                                                      <  1:976620/139‑1 (MQ=255)
                                                             ggcgCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgt                                                                      <  2:710275/139‑1 (MQ=255)
                                                             ggcgCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgt                                                                      <  1:631608/139‑1 (MQ=255)
                                                               cgcgACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgtgt                                                                    <  2:55472/139‑1 (MQ=255)
                                                               cgcgACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgtgt                                                                    <  2:312499/139‑1 (MQ=255)
                                                                gcgACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGgcg                                                                                                                      <  2:815328/88‑1 (MQ=255)
                                                                gcgACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGgcg                                                                                                                      >  1:815328/1‑88 (MQ=255)
                                                                 cgACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGtct                                                                  >  2:1069668/1‑139 (MQ=255)
                                                                 cgACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGtct                                                                  >  1:464656/1‑139 (MQ=255)
                                                                 cgACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGtct                                                                  <  2:454011/139‑1 (MQ=255)
                                                                 cgACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGtct                                                                  >  1:675112/1‑139 (MQ=255)
                                                                         tGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGt                                                          >  2:710398/1‑139 (MQ=255)
                                                                              ccGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAaagaag                                                                                          <  2:469184/102‑1 (MQ=255)
                                                                              ccGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAaagaag                                                                                          >  1:469184/1‑102 (MQ=255)
                                                                                 gCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGaa                                                                                                        <  2:1357540/85‑1 (MQ=255)
                                                                                 gCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGaa                                                                                                        >  1:1357540/1‑85 (MQ=255)
                                                                                 gCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTa                                                                                >  2:1404514/1‑109 (MQ=255)
                                                                                 gCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTa                                                                                <  1:1404514/109‑1 (MQ=255)
                                                                                 gCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTgg                                                  >  1:1510694/1‑139 (MQ=255)
                                                                                  cACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGgc                                                 <  1:575279/139‑1 (MQ=255)
                                                                                   aCTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGgcg                                                <  1:236884/139‑1 (MQ=255)
                                                                                   aCTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGgcg                                                <  1:862165/139‑1 (MQ=255)
                                                                                   aCTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGgcg                                                >  2:862165/1‑139 (MQ=255)
                                                                                     tAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCa                                              <  1:273235/139‑1 (MQ=255)
                                                                                       aaaCCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATc                                            <  2:675112/139‑1 (MQ=255)
                                                                                        aaCCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATcc                                           <  2:202664/139‑1 (MQ=255)
                                                                                         aCCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCa                                          <  1:241795/139‑1 (MQ=255)
                                                                                          ccGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCaa                                         <  2:464656/139‑1 (MQ=255)
                                                                                          ccGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCaa                                         <  2:1539176/139‑1 (MQ=255)
                                                                                             aTACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGc                                                                                  <  1:170571/95‑1 (MQ=255)
                                                                                             aTACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGc                                                                                  >  2:170571/1‑95 (MQ=255)
                                                                                             aTACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAacac                                      <  1:1069668/139‑1 (MQ=255)
                                                                                                   cGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGa                                >  2:1474243/1‑139 (MQ=255)
                                                                                                               tGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGgcg                                                                                                                      >  1:288852/1‑41 (MQ=255)
                                                                                                               tGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGgcg                                                                                                                      <  2:288852/41‑1 (MQ=255)
                                                                                                               tGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAgcg                                                                                        <  2:223903/71‑1 (MQ=255)
                                                                                                               tGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAgcg                                                                                        >  1:223903/1‑71 (MQ=255)
                                                                                                                gCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCg                                                                                                          >  1:68185/1‑52 (MQ=255)
                                                                                                                gCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCg                                                                                                          <  2:68185/52‑1 (MQ=255)
                                                                                                                gCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTAt                                  <  2:1321151/124‑1 (MQ=255)
                                                                                                                gCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTAt                                  >  1:1321151/1‑124 (MQ=255)
                                                                                                                   acaGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTcgc                <  1:710398/139‑1 (MQ=255)
                                                                                                                   acaGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTcgc                >  1:1145674/1‑139 (MQ=255)
                                                                                                                       tCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCt            >  2:147178/1‑139 (MQ=255)
                                                                                                                           tttCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAg                                                                                                       <  2:210865/44‑1 (MQ=255)
                                                                                                                           tttCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAg                                                                                                       >  1:210865/1‑44 (MQ=255)
                                                                                                                                 ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgtgt                                                                    >  2:324137/1‑73 (MQ=255)
                                                                                                                                 ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgtgt                                                                    <  1:324137/73‑1 (MQ=255)
                                                                                                                                 ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGa                        >  1:53312/1‑117 (MQ=255)
                                                                                                                                 ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGa                        <  2:53312/117‑1 (MQ=255)
                                                                                                                                 ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGa  <  2:1145674/139‑1 (MQ=255)
                                                                                                                                 ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGa  >  1:1603210/1‑139 (MQ=255)
                                                                                                                                   gtgATCGCCGAAGAGATGGCGCTGCATCAGCCGAAg                                                                                                       <  2:1334589/36‑1 (MQ=255)
                                                                                                                                   gtgATCGCCGAAGAGATGGCGCTGCATCAGCCGAAg                                                                                                       >  1:1334589/1‑36 (MQ=37)
                                                                                                                                   gtgATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgtgt                                                                    <  2:1096776/71‑1 (MQ=255)
                                                                                                                                   gtgATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgtgt                                                                    >  1:1096776/1‑71 (MQ=255)
                                                                                                                                       |                                                                                                                                    
TCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGA  >  NC_000913/167608‑167875

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