Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 338,763 C→T P147S (CCG→TCG)  yahG → DUF1116 domain‑containing protein YahG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913338,7630CT100.0% 970.6 / NA 239P147S (CCG→TCG) yahGDUF1116 domain‑containing protein YahG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (77/162);  total (77/162)

Only 100 of 241 total aligned reads displayed.

TGAAAGGCGCGGTCACCGGAGCGCTGGTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAAAAACAAAACCTACGGCAACATC  >  NC_000913/338665‑338858
                                                                                                  |                                                                                               
  aaaGGCGCGGTCACCGGAGCGCTGGTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCgtg                                                                                             >  1:580516/1‑101 (MQ=255)
  aaaGGCGCGGTCACCGGAGCGCTGGTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCgtg                                                                                             <  2:580051/101‑1 (MQ=255)
   aaGGCGCGGTCACCGGAGCGCTGGTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCgtgt                                                                                            <  1:579901/101‑1 (MQ=255)
     ggCGCGGTCACCGGAGCGCTGGTGTTGGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCa                                                                                          <  2:580162/101‑1 (MQ=255)
        gcgGTCACCGGAGCGCTGGTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGa                                                                                       >  1:580523/1‑101 (MQ=255)
          ggTCACCGGAGCGCTGGTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGc                                                                                     >  2:580525/1‑101 (MQ=255)
          ggTCACCGGAGCGCTGGTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGc                                                                                     <  2:580062/101‑1 (MQ=255)
           gTCACCGGAGCGCTGGTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCAc                                                                                         <  2:580153/96‑1 (MQ=255)
            tCACCGGAGCGCTGGTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCAc                                                                                   <  2:579839/101‑1 (MQ=255)
             cACCGGAGCGCTGGTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCCCCTTCTCGTCGTGTCACGAGCACg                                                                                  <  2:580011/101‑1 (MQ=255)
             cACCGGAGCGCTGGTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACg                                                                                  <  2:580231/101‑1 (MQ=255)
             cACCGGAGCGCTGGTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACg                                                                                  >  2:580530/1‑101 (MQ=255)
                cGGAGCGCTGGTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACt                                                                               >  1:580535/1‑101 (MQ=255)
                  gAGCGCTGGTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGc                                                                             <  2:580288/101‑1 (MQ=255)
                   aGCGCTGGTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGAc                                                                                >  1:580540/1‑97 (MQ=255)
                    gcgcTGGTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCGCCTTCTCGTCGTGTCACGAGCACGACTGCGt                                                                           <  1:580096/101‑1 (MQ=255)
                       cTGGTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTggg                                                                        <  1:580309/101‑1 (MQ=255)
                       cTGGTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTggg                                                                        >  2:580543/1‑101 (MQ=255)
                        tGGTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGa                                                                       <  2:580092/101‑1 (MQ=255)
                         ggTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGAt                                                                      <  2:579970/101‑1 (MQ=255)
                         ggTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGAt                                                                      <  2:579981/101‑1 (MQ=255)
                         ggTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGAt                                                                      <  1:579780/101‑1 (MQ=255)
                            gtTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCAc                                                                                   >  2:580551/1‑85 (MQ=255)
                            gtTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGa                                                                   >  1:580550/1‑101 (MQ=255)
                             ttCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGAt                                                                  <  2:579897/101‑1 (MQ=255)
                               cGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATgg                                                                <  1:580257/101‑1 (MQ=255)
                               cGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATgg                                                                <  1:579958/101‑1 (MQ=255)
                                gAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGc                                                               >  1:580557/1‑101 (MQ=255)
                                  aGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCgg                                                             >  2:580561/1‑101 (MQ=255)
                                   ggACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCggg                                                            >  1:580562/1‑101 (MQ=255)
                                     aCTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCGCCTTCTCGTCGTGTCACGAGCACGAc                                                                                <  1:579976/79‑1 (MQ=255)
                                     aCTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGgtg                                                          <  1:580020/101‑1 (MQ=255)
                                     aCTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGgtg                                                          <  2:579949/101‑1 (MQ=255)
                                       tGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGtt                                                        <  2:580216/101‑1 (MQ=255)
                                          cGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTAcc                                                     <  1:580221/101‑1 (MQ=255)
                                           gAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCt                                                    <  1:580157/101‑1 (MQ=255)
                                            aaaGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTa                                                       <  2:580211/97‑1 (MQ=255)
                                            aaaGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTc                                                   <  2:580397/101‑1 (MQ=255)
                                              aGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTg                                                                              >  1:580571/1‑72 (MQ=255)
                                               gATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGc                                                >  2:580572/1‑101 (MQ=255)
                                                aTCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGcaccac                                                                                >  1:580574/1‑68 (MQ=255)
                                                  ctcGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGAc                                                                                <  2:580016/66‑1 (MQ=255)
                                                  ctcGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTc                                             >  1:580578/1‑101 (MQ=255)
                                                   tcGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGa                                                                       <  1:580081/74‑1 (MQ=255)
                                                    cgacgaGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTc                                                                                           >  2:580581/1‑53 (MQ=255)
                                                       cgaGGCGGCTGAACTGGCGGCTTCCGGGGAGATCCCCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGt                                        <  2:580102/101‑1 (MQ=255)
                                                       cgaGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGt                                        <  2:580174/101‑1 (MQ=255)
                                                       cgaGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGt                                        <  2:580111/101‑1 (MQ=255)
                                                        gaGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATg                                                                 >  2:580585/1‑75 (MQ=255)
                                                        gaGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGtt                                       <  2:580262/101‑1 (MQ=255)
                                                         aGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGttt                                      >  1:580586/1‑101 (MQ=255)
                                                          ggcggcTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTa                                     <  2:580129/101‑1 (MQ=255)
                                                          ggcggcTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTa                                     <  1:580212/101‑1 (MQ=255)
                                                           gcggcTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTAt                                    <  2:580186/101‑1 (MQ=255)
                                                           gcggcTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTAt                                    <  1:580190/101‑1 (MQ=255)
                                                           gcggcTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTAt                                    <  1:579701/101‑1 (MQ=255)
                                                           gcggcTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTAt                                    <  2:580312/101‑1 (MQ=255)
                                                           gcggcTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTAt                                    <  1:580077/101‑1 (MQ=255)
                                                            cggcTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATg                                   <  1:580112/101‑1 (MQ=255)
                                                              gcTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTgg                                                                         <  1:580008/61‑1 (MQ=255)
                                                              gcTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGca                                 <  1:580053/101‑1 (MQ=255)
                                                              gcTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGGGTTACCTCGGCCTCGATGTTTATGca                                 >  2:580589/1‑101 (MQ=255)
                                                               cTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGcac                                <  2:580107/101‑1 (MQ=255)
                                                                tGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGc                                                >  2:580592/1‑84 (MQ=255)
                                                                tGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGcaca                               <  2:580180/101‑1 (MQ=255)
                                                                 gAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACAt                              >  1:580595/1‑101 (MQ=255)
                                                                 gAACTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACAt                              <  1:580032/101‑1 (MQ=255)
                                                                  aaCTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATc                             <  2:580116/101‑1 (MQ=255)
                                                                   aCTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCg                            <  1:579941/101‑1 (MQ=255)
                                                                    cTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGca                                 <  1:579966/95‑1 (MQ=255)
                                                                    cTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATc                             <  2:579794/99‑1 (MQ=255)
                                                                    cTGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGt                           <  2:580144/101‑1 (MQ=255)
                                                                     tGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTa                          <  1:580235/101‑1 (MQ=255)
                                                                     tGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTa                          >  2:580598/1‑101 (MQ=255)
                                                                     tGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTa                          <  1:580121/101‑1 (MQ=255)
                                                                     tGGCGGCTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTa                          >  1:580599/1‑101 (MQ=255)
                                                                       gcggcTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTaaa                        <  1:580203/101‑1 (MQ=255)
                                                                        cggcTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTaaaa                       <  2:580041/101‑1 (MQ=255)
                                                                         ggcTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCg                                                                            >  2:580602/1‑47 (MQ=255)
                                                                         ggcTTCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTaaaaa                      <  2:580332/101‑1 (MQ=255)
                                                                            ttCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAaaaaca                   <  1:580358/101‑1 (MQ=255)
                                                                            ttCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAaaaaca                   >  2:580605/1‑101 (MQ=255)
                                                                             tCCGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGt                                        <  2:580250/79‑1 (MQ=255)
                                                                              ccGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCt                                                    <  1:580125/66‑1 (MQ=255)
                                                                              ccGGGGAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAaaaacaaa                 >  1:580607/1‑101 (MQ=255)
                                                                                ggggAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAaaaacaaaac               >  2:580609/1‑101 (MQ=255)
                                                                                  ggAGATCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAAAAACAAAACCt             <  1:580176/101‑1 (MQ=255)
                                                                                   gagaTCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAAAAACAAAACCTa            >  2:580612/1‑101 (MQ=255)
                                                                                    agaTCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAAAAACAAAACCTAc           <  2:579995/101‑1 (MQ=255)
                                                                                      aTCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATc                                                                     >  2:580615/1‑41 (MQ=255)
                                                                                       tCACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAAAAACAAAACCTACGGc        <  1:580334/101‑1 (MQ=255)
                                                                                        cACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAAAAACAAAACCTACGGCa       <  1:579820/101‑1 (MQ=255)
                                                                                        cACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAAAAACAAAACCTACGGCa       <  1:580427/101‑1 (MQ=255)
                                                                                        cACCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAAAAACAAAACCTACGGCa       >  1:580616/1‑101 (MQ=255)
                                                                                         aCCTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAAAAACAAAACCTACGGCaa      <  1:580027/101‑1 (MQ=255)
                                                                                          ccTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAAAAACAAAACCTACGGCAAc     <  2:580238/101‑1 (MQ=255)
                                                                                           cTTCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAAAAACAAAACCTACGGCa       <  2:580069/98‑1 (MQ=255)
                                                                                            ttCTCGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAAAAACAAAACCTACGGCAACAt   >  2:580618/1‑101 (MQ=255)
                                                                                             tctcGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGAt                                          <  2:580083/61‑1 (MQ=255)
                                                                                             tctcGTCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAAAAACAAAACCTACGGCAACATc  <  2:580482/101‑1 (MQ=255)
                                                                                                  |                                                                                               
TGAAAGGCGCGGTCACCGGAGCGCTGGTGTTCGAAGGACTGGCGAAAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAAAAACAAAACCTACGGCAACATC  >  NC_000913/338665‑338858

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