Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 232,218 A→G D98G (GAC→GGC)  yafE → putative biotin synthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913232,2180AG100.0% 563.2 / NA 180D98G (GAC→GGC) yafEputative biotin synthesis protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (99/81);  total (99/81)

Only 100 of 187 total aligned reads displayed.

TTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACTACGCCA  >  NC_000913/232129‑232307
                                                                                         |                                                                                         
ttATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGt                                                                                        >  1:3291433/1‑93 (MQ=255)
 tATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAGTAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTa                                                                                       >  1:4297820/1‑93 (MQ=255)
   ttCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGAGCGTAAt                                                                                     <  2:4608212/93‑1 (MQ=255)
   ttCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGAGAATAGGATATTGAAACCTGGCGGCAGGCTGATTGTGATGGGCGTAAt                                                                                     >  1:1824508/1‑93 (MQ=255)
    tCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATg                                                                                    <  1:4433002/93‑1 (MQ=255)
    tCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATg                                                                                    <  1:2216191/93‑1 (MQ=255)
         ccATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTATGCTGATTGTGATGGGCGTAATGTCTcc                                                                               >  2:4265089/1‑93 (MQ=255)
         ccATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTcc                                                                               >  1:3969581/1‑93 (MQ=255)
         ccATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTcc                                                                               >  2:1932195/1‑93 (MQ=255)
         ccATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTcc                                                                               >  2:2375664/1‑93 (MQ=255)
           aTCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCgg                                                                             >  1:2768593/1‑93 (MQ=255)
            tCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGGCTCCggg                                                                            >  1:2551403/1‑93 (MQ=255)
              aCTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTc                                                                          >  2:4675454/1‑93 (MQ=255)
                 ggCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCAcc                                                                       >  1:2289597/1‑93 (MQ=255)
                 ggCATGATGTTGGTGCAGCACTGCGAGAAGTGAAAAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCAcc                                                                       <  1:782082/93‑1 (MQ=255)
                   cATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCa                                                                     >  1:4071661/1‑93 (MQ=255)
                    atgatgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAg                                                                    <  1:1150365/93‑1 (MQ=255)
                     tgatgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGt                                                                   <  1:1501105/93‑1 (MQ=255)
                      gatgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTg                                                                  <  1:2388005/93‑1 (MQ=255)
                      gatgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTg                                                                  >  1:3446820/1‑93 (MQ=255)
                         gTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTgcgc                                                               >  2:1172476/1‑93 (MQ=255)
                            ggTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGAc                                                            >  2:2229511/1‑93 (MQ=255)
                             gTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACa                                                           <  1:2789053/93‑1 (MQ=255)
                             gTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACa                                                           <  1:523227/93‑1 (MQ=255)
                               gcagcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATc                                                         <  1:2093860/93‑1 (MQ=255)
                               gcagcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATc                                                         <  2:1472872/93‑1 (MQ=255)
                               gcagcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATc                                                         <  2:4146230/93‑1 (MQ=255)
                               gcagcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATc                                                         <  1:422055/93‑1 (MQ=255)
                                cagcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCt                                                        <  1:241912/93‑1 (MQ=255)
                                cagcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCt                                                        >  2:1609477/1‑93 (MQ=255)
                                 agcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTg                                                       >  2:2840765/1‑93 (MQ=255)
                                 agcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTg                                                       <  1:642300/93‑1 (MQ=255)
                                  gcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTgg                                                      <  2:836959/93‑1 (MQ=255)
                                      tGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTAc                                                  >  2:4262514/1‑93 (MQ=255)
                                       gCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACa                                                 >  2:3544765/1‑93 (MQ=255)
                                        cGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAg                                                <  1:1523771/93‑1 (MQ=255)
                                         gagaAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGa                                               >  2:2479535/1‑93 (MQ=255)
                                           gaAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACg                                             <  2:1394461/93‑1 (MQ=255)
                                           gaAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACg                                             >  1:1303306/1‑93 (MQ=255)
                                           gaAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACg                                             >  2:1656673/1‑93 (MQ=255)
                                           gaAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACg                                             <  1:4483666/93‑1 (MQ=255)
                                           gaAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACg                                             >  1:2469422/1‑93 (MQ=255)
                                           gaAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACg                                             >  1:2291258/1‑93 (MQ=255)
                                            aaGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACgg                                            >  2:1091532/1‑93 (MQ=255)
                                            aaGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACgg                                            >  2:2955268/1‑93 (MQ=255)
                                            aaGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACgg                                            >  2:903051/1‑93 (MQ=255)
                                            aaGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACgg                                            <  1:389856/93‑1 (MQ=255)
                                             aGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGt                                           <  2:2699738/93‑1 (MQ=255)
                                              gTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTa                                          <  1:4147460/93‑1 (MQ=255)
                                              gTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTa                                          <  2:3890450/93‑1 (MQ=255)
                                              gTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTa                                          <  2:4350861/93‑1 (MQ=255)
                                              atGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTa                                          <  1:4398475/92‑1 (MQ=255)
                                               tGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAg                                         >  1:3965721/1‑93 (MQ=255)
                                                 aaTAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGaa                                       >  1:384015/1‑93 (MQ=255)
                                                 aaTAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGaa                                       >  1:3720643/1‑93 (MQ=255)
                                                  aTAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAg                                      <  1:3557425/93‑1 (MQ=255)
                                                   tAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGc                                     <  1:1248919/93‑1 (MQ=255)
                                                   tAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGc                                     >  1:3694292/1‑93 (MQ=255)
                                                    aGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCa                                    <  2:3768774/93‑1 (MQ=255)
                                                     ggATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCGAGTGCGCGACATCTGGTTACAGACGGTAGAAGCAt                                   <  1:2003422/93‑1 (MQ=255)
                                                     ggATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCAt                                   >  1:1435965/1‑93 (MQ=255)
                                                     ggATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCAt                                   <  1:1649045/93‑1 (MQ=255)
                                                     ggATATTGAAACCTGGCGGTAGGCTGATTGCGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGCAGAAGCAt                                   >  2:3168714/1‑93 (MQ=255)
                                                      gATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGCCGGTAGAAGCAtt                                  <  2:3452211/93‑1 (MQ=255)
                                                      gATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCAtt                                  <  2:3414435/93‑1 (MQ=255)
                                                      gATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCAtt                                  <  1:2707445/93‑1 (MQ=255)
                                                       atatTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTa                                 <  1:3011100/93‑1 (MQ=255)
                                                        tatTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAc                                >  2:3089707/1‑93 (MQ=255)
                                                        tatTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAc                                >  2:3063268/1‑93 (MQ=255)
                                                        tatTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAc                                <  1:604792/93‑1 (MQ=255)
                                                            gAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGa                            <  2:1785831/93‑1 (MQ=255)
                                                            gAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGa                            <  1:3516846/93‑1 (MQ=255)
                                                            gAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGa                            <  2:3618690/93‑1 (MQ=255)
                                                             aaaCCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGAt                           >  1:757481/1‑93 (MQ=255)
                                                             aaaCCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGAt                           >  1:1015512/1‑93 (MQ=255)
                                                              aaCCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATa                          >  2:362538/1‑93 (MQ=255)
                                                              aaCCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATa                          >  1:2860651/1‑93 (MQ=255)
                                                                ccTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATAcc                        <  2:3977234/93‑1 (MQ=255)
                                                                 cTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACct                       >  2:2568813/1‑93 (MQ=255)
                                                                  tGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACctc                      <  2:491729/93‑1 (MQ=255)
                                                                  cgggggTAGGCGGATTGTGATGGGCGTAATGTCTCCGGTTCAGCTGGTGCGTGGCATCTGGTTGCAGACGGTAGAAGCATTACGGGATACctc                      <  2:3945692/89‑1 (MQ=255)
                                                                    gCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCtctc                    >  2:1974629/1‑93 (MQ=255)
                                                                    gCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCtctc                    <  1:4090910/93‑1 (MQ=255)
                                                                     cGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCa                   <  1:3659460/93‑1 (MQ=255)
                                                                       gTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACg                 >  2:1870350/1‑93 (MQ=255)
                                                                          ggCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTAc              >  1:556136/1‑93 (MQ=255)
                                                                           gCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACg             <  2:1655849/93‑1 (MQ=255)
                                                                           gCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACg             <  1:3518988/93‑1 (MQ=255)
                                                                            cTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGa            >  1:236843/1‑93 (MQ=255)
                                                                             tGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGaa           >  2:3228697/1‑93 (MQ=255)
                                                                              gATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGaaa          >  2:2404720/1‑93 (MQ=255)
                                                                               aTTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAAc         >  2:4075936/1‑93 (MQ=255)
                                                                                ttGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACt        >  2:2166184/1‑93 (MQ=255)
                                                                                ttGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACt        >  1:989794/1‑93 (MQ=255)
                                                                                ttGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACt        >  1:4454033/1‑93 (MQ=255)
                                                                                ttGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACt        >  1:110218/1‑93 (MQ=255)
                                                                                ttGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACt        >  1:1061787/1‑93 (MQ=255)
                                                                                     aTGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACTACGcc   >  2:2029748/1‑93 (MQ=255)
                                                                                     aTGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACTACGcc   <  1:3176286/93‑1 (MQ=255)
                                                                                     aTGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACTACGcc   >  1:2176179/1‑93 (MQ=255)
                                                                                         |                                                                                         
TTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACTACGCCA  >  NC_000913/232129‑232307

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