Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,616,692 C→T G81S (GGC→AGC)  ydeB ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,616,6920CT100.0% 904.6 / NA 241G81S (GGC→AGC) ydeBhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (133/108);  total (133/108)

Only 100 of 260 total aligned reads displayed.

GCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATA  >  NC_000913/1616602‑1616781
                                                                                          |                                                                                         
gCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgctgc                                                                                         >  1:6393053/1‑93 (MQ=255)
gCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgctgc                                                                                         <  1:1714436/93‑1 (MQ=255)
gCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgctgc                                                                                         <  1:386314/93‑1 (MQ=255)
gCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgctgc                                                                                         >  1:2269548/1‑93 (MQ=255)
 ctctTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgctgcc                                                                                        >  1:1485562/1‑93 (MQ=255)
 ctctTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAgctgcc                                                                                        <  1:3373972/93‑1 (MQ=255)
  tctTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCt                                                                                       >  2:1347300/1‑93 (MQ=255)
  tctTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCt                                                                                       >  2:5742267/1‑93 (MQ=255)
  tctTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCt                                                                                       >  2:3830507/1‑93 (MQ=255)
  tctTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCt                                                                                       >  2:4878387/1‑93 (MQ=255)
     tCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCa                                                                                    >  1:3972060/1‑93 (MQ=255)
     tCCGACTTGCTTTTCGCCTCCGGTGAAGCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCa                                                                                    <  2:294446/93‑1 (MQ=255)
      ccGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCaa                                                                                   >  1:5930224/1‑93 (MQ=255)
      ccGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCaa                                                                                   >  1:4607137/1‑93 (MQ=255)
      ccGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCaa                                                                                   <  1:1977569/93‑1 (MQ=255)
        gACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAAtt                                                                                 >  1:6275973/1‑93 (MQ=255)
        gACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAAtt                                                                                 >  1:3174846/1‑93 (MQ=255)
         aCTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTc                                                                                >  1:4061191/1‑93 (MQ=255)
          cTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAGCCGATAAACGCGACGATTAAGCTGCCTGCAATTCg                                                                               <  1:6307014/93‑1 (MQ=255)
          cTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCg                                                                               >  2:5961170/1‑93 (MQ=255)
           ttGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGc                                                                              <  2:3500741/93‑1 (MQ=255)
            tGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCa                                                                             <  2:5641863/93‑1 (MQ=255)
            tGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCa                                                                             <  2:2566025/93‑1 (MQ=255)
            tGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCa                                                                             <  1:4236552/93‑1 (MQ=255)
             gCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAg                                                                            <  1:6011242/93‑1 (MQ=255)
               ttttCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGAc                                                                          <  1:1775310/93‑1 (MQ=255)
               ttttCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGAc                                                                          <  2:2191145/93‑1 (MQ=255)
               ttttCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGAc                                                                          >  1:2193866/1‑93 (MQ=255)
                   cGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCggg                                                                      <  1:4348784/93‑1 (MQ=255)
                    gCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGa                                                                     >  1:1087160/1‑93 (MQ=255)
                     ccTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGaa                                                                    <  2:4976541/93‑1 (MQ=255)
                     ccTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGaa                                                                    <  2:6620418/93‑1 (MQ=255)
                      cTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAAt                                                                   >  2:1700971/1‑93 (MQ=255)
                      cTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAAt                                                                   >  2:2861738/1‑93 (MQ=255)
                      cTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAAt                                                                   >  2:362584/1‑93 (MQ=255)
                       tCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATc                                                                  >  1:1220090/1‑93 (MQ=255)
                         cGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAACCGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGa                                                                <  1:673170/93‑1 (MQ=255)
                          ggTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGaa                                                               <  1:1214958/93‑1 (MQ=255)
                           gTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGaaa                                                              >  2:3379058/1‑93 (MQ=255)
                           gTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGaaa                                                              >  2:3731794/1‑93 (MQ=255)
                            tGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGCCGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAAt                                                             >  2:5776559/1‑93 (MQ=255)
                             gAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATg                                                            >  2:4100112/1‑93 (MQ=255)
                             gAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATg                                                            >  1:2327294/1‑93 (MQ=255)
                                tcaatcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCa                                                         >  1:4791990/1‑93 (MQ=255)
                                tcaatcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCa                                                         <  2:2101129/93‑1 (MQ=255)
                                 caatcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCaa                                                        <  2:3240734/93‑1 (MQ=255)
                                   atcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAAt                                                      >  1:234451/1‑93 (MQ=255)
                                   atcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAAt                                                      >  1:2884101/1‑93 (MQ=255)
                                   atcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAAt                                                      >  2:2491442/1‑93 (MQ=255)
                                   atcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAAt                                                      >  1:4458327/1‑93 (MQ=255)
                                      gCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTa                                                   <  1:490706/93‑1 (MQ=255)
                                      gCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTa                                                   >  2:2014317/1‑93 (MQ=255)
                                      gCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTa                                                   <  1:4772033/93‑1 (MQ=255)
                                       cTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTAt                                                  >  2:393032/1‑93 (MQ=255)
                                        tttCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATc                                                 >  2:885513/1‑93 (MQ=255)
                                        tttCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATc                                                 <  1:5161506/93‑1 (MQ=255)
                                        tttCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATc                                                 <  2:5054374/93‑1 (MQ=255)
                                        tttCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATc                                                 >  1:3105071/1‑93 (MQ=255)
                                         ttCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATcc                                                <  1:1014159/93‑1 (MQ=255)
                                            tGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATc                                             <  2:5968867/93‑1 (MQ=255)
                                            tGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATc                                             >  1:847779/1‑93 (MQ=255)
                                                gCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCa                                         <  2:6280287/93‑1 (MQ=255)
                                                 cGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAg                                        <  2:4511512/93‑1 (MQ=255)
                                                   gAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGtt                                      <  2:2116995/93‑1 (MQ=255)
                                                    aaaGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTg                                     >  1:1295863/1‑93 (MQ=255)
                                                    aaaGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTg                                     <  1:3944518/93‑1 (MQ=255)
                                                     aaGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGc                                    >  2:6481981/1‑93 (MQ=255)
                                                       gagCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGccc                                  >  2:2659128/1‑93 (MQ=255)
                                                       gagCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGccc                                  >  2:3038340/1‑93 (MQ=255)
                                                             ccGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTa                            <  1:3712186/93‑1 (MQ=255)
                                                             ccGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTa                            <  1:3443967/93‑1 (MQ=255)
                                                              cGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTaa                           <  1:3897284/93‑1 (MQ=255)
                                                               gAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAAt                          <  1:3300942/93‑1 (MQ=255)
                                                               gAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAAt                          >  1:5989951/1‑93 (MQ=255)
                                                                aaaaCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATa                         <  2:6111963/93‑1 (MQ=255)
                                                                  aaCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACg                       >  1:5421501/1‑93 (MQ=255)
                                                                  aaCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACg                       >  1:5637682/1‑93 (MQ=255)
                                                                   aCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGc                      >  2:1025664/1‑93 (MQ=255)
                                                                     cGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGcca                    >  2:2780773/1‑93 (MQ=255)
                                                                      gATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGccac                   >  1:538115/1‑93 (MQ=255)
                                                                      gATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGccac                   >  2:3659328/1‑93 (MQ=255)
                                                                      gATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGccac                   >  2:3166419/1‑93 (MQ=255)
                                                                        tAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGccacca                 >  1:558389/1‑93 (MQ=255)
                                                                        tAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGccacca                 >  2:1653729/1‑93 (MQ=255)
                                                                          aaCGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACcatc               <  1:2455375/93‑1 (MQ=255)
                                                                           aCGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACcatca              <  2:1398576/93‑1 (MQ=255)
                                                                           aCGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACcatca              >  1:2576759/1‑93 (MQ=255)
                                                                             gcgACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatg            <  1:3752398/93‑1 (MQ=255)
                                                                              cgacgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatga           <  2:5455871/93‑1 (MQ=255)
                                                                               gacgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatgat          <  2:6368962/93‑1 (MQ=255)
                                                                               gacgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGCATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatgat          <  2:4697035/93‑1 (MQ=255)
                                                                                acgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatgatg         >  2:1432845/1‑93 (MQ=255)
                                                                                 cgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATgg        <  2:1231055/93‑1 (MQ=255)
                                                                                  gaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGc       >  2:4071783/1‑93 (MQ=255)
                                                                                      aaGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAAt   <  1:4269262/93‑1 (MQ=255)
                                                                                      aaGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAAt   <  2:2343225/93‑1 (MQ=255)
                                                                                      aaGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAAt   <  1:5112694/93‑1 (MQ=255)
                                                                                       aGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATa  <  1:1536126/93‑1 (MQ=255)
                                                                                       aGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATa  >  2:4436453/1‑93 (MQ=255)
                                                                                       aGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATa  <  1:2674795/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATA  >  NC_000913/1616602‑1616781

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: