Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,877,994 C→T R709H (CGT→CAT)  gyrB ← DNA gyrase subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,877,9940CT100.0% 1129.4 / NA 268R709H (CGT→CAT) gyrBDNA gyrase subunit B
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (115/153);  total (115/153)

Only 100 of 268 total aligned reads displayed.

GACCTTTATAACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCACGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGCCACCGGTGATAAACTCGTGATCCAG  >  NC_000913/3877896‑3878091
                                                                                                  |                                                                                                 
gACCTTTATAACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGt                                                                                                 >  1:6606450/1‑101 (MQ=255)
   cTTTATAACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCg                                                                                              >  1:6606455/1‑101 (MQ=255)
    ttgataACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGa                                                                                             >  1:6606458/4‑101 (MQ=255)
     ttATAACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGAt                                                                                            >  2:6606461/1‑101 (MQ=255)
      tataACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATa                                                                                           >  2:6606465/1‑101 (MQ=255)
      tataACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATa                                                                                           >  1:6606463/1‑101 (MQ=255)
      tataACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATa                                                                                           <  1:6605820/101‑1 (MQ=255)
       ataACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATaa                                                                                          <  1:6605686/101‑1 (MQ=255)
        taACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATaaa                                                                                         <  1:6606065/101‑1 (MQ=255)
        taACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATaaa                                                                                         >  1:6606470/1‑101 (MQ=255)
         aaCGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAAc                                                                                        >  1:6606472/1‑101 (MQ=255)
           cGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGc                                                                                      <  1:6606197/101‑1 (MQ=255)
               ggatggaGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATaa                                                                                          <  1:6605925/93‑1 (MQ=255)
                gatggaGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCAtctt                                                                                 >  2:6606481/1‑101 (MQ=255)
                  tggaGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCAtcttct                                                                               >  2:6606483/1‑101 (MQ=255)
                    gagaGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCttcttc                                                                             <  1:6606001/101‑1 (MQ=255)
                    gagaGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCttcttc                                                                             <  1:6606165/101‑1 (MQ=255)
                    gagaGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCttcttc                                                                             >  1:6606484/1‑101 (MQ=255)
                         gCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCa                                                                        >  2:6606491/1‑101 (MQ=255)
                           cgcgACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAg                                                                      >  1:6606494/1‑101 (MQ=255)
                            gcgACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGc                                                                     >  1:6606497/1‑101 (MQ=255)
                              gACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCaa                                                                       >  2:6606501/1‑97 (MQ=255)
                               aCGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCAc                                                                  >  1:6606502/1‑101 (MQ=255)
                                  ggACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCa                                                               >  1:6606508/1‑101 (MQ=255)
                                   gACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAg                                                              >  1:6606511/1‑101 (MQ=255)
                                   gACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAg                                                              >  2:6606510/1‑101 (MQ=255)
                                     ctctTTCCCCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATa                                                                                           <  2:6606238/70‑1 (MQ=255)
                                     ctctTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAAc                                                                                        <  2:6605902/73‑1 (MQ=255)
                                     ctctTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAg                                                                          <  2:6606133/87‑1 (MQ=255)
                                     ctctTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGtt                                                            <  1:6605721/101‑1 (MQ=255)
                                     ctctTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGtt                                                            <  1:6606070/101‑1 (MQ=255)
                                      tctTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCAtctt                                                                                 <  1:6605889/79‑1 (MQ=255)
                                      tctTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCa                                                                           <  1:6605849/85‑1 (MQ=255)
                                        tttCCCCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTtc                                                          <  2:6606297/100‑1 (MQ=255)
                                        tttCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCTAGCAAGCCACGCAGTTtct                                                         >  1:6606518/1‑101 (MQ=255)
                                        tttCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCaa                                                                       <  1:6605942/87‑1 (MQ=255)
                                        tttCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCAc                                                                  <  1:6605900/92‑1 (MQ=255)
                                        tttCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTtct                                                         >  2:6606517/1‑101 (MQ=255)
                                         ttCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTtct                                                         >  1:6606521/1‑100 (MQ=255)
                                          tccccagccagTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTc                                                                                               <  2:6605909/58‑1 (MQ=255)
                                            aCCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAg                                                              <  1:6605767/92‑1 (MQ=255)
                                            aCCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCAcc                                                     <  1:6605970/101‑1 (MQ=255)
                                             ccagccagTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAAc                                                                                        <  2:6605789/65‑1 (MQ=255)
                                             ccagccagTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTcc                                                                            <  2:6605991/77‑1 (MQ=255)
                                             ccagccagTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCa                                                                           <  2:6605930/78‑1 (MQ=255)
                                             ccagccagTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGc                                                                         <  2:6605756/80‑1 (MQ=255)
                                             ccagccagTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCa                                                                   <  1:6605978/86‑1 (MQ=255)
                                             ccagccagTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCAccc                                                    <  1:6606028/101‑1 (MQ=255)
                                              cagccagTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGc                                                                     <  2:6606104/83‑1 (MQ=255)
                                                gccagTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCa                                                               <  1:6605934/87‑1 (MQ=255)
                                                gccagTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGc                                                 <  2:6606045/101‑1 (MQ=255)
                                                 ccagTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGc                                                                     >  1:6606527/1‑80 (MQ=255)
                                                 ccagTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCAc                                                      <  2:6606150/95‑1 (MQ=255)
                                                 ccagTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCg                                                >  2:6606526/1‑101 (MQ=255)
                                                  cagTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAg                                                              <  2:6606121/86‑1 (MQ=255)
                                                    gTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGc                                             <  1:6606134/101‑1 (MQ=255)
                                                      ccAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAg                                           <  2:6605962/101‑1 (MQ=255)
                                                      ccAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAg                                           <  2:6606173/101‑1 (MQ=255)
                                                      ccAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAg                                           <  1:6606081/101‑1 (MQ=255)
                                                      ccAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAg                                           <  1:6606153/101‑1 (MQ=255)
                                                       cAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGa                                          >  1:6606531/1‑101 (MQ=255)
                                                        aGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGAt                                         <  1:6605487/101‑1 (MQ=255)
                                                        aGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGAt                                         >  2:6606532/1‑101 (MQ=255)
                                                           gcCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATacg                                      <  2:6606253/101‑1 (MQ=255)
                                                            ccTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTg                                              <  2:6606052/92‑1 (MQ=255)
                                                            ccTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATacga                                     >  1:6606536/1‑101 (MQ=255)
                                                             cTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATacgac                                    <  2:6606111/101‑1 (MQ=255)
                                                              tGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATacgacg                                   >  1:6606539/1‑101 (MQ=255)
                                                               gCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATaaaag                                                                                       >  2:6606541/1‑46 (MQ=255)
                                                                 tCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGc                                             <  1:6605948/88‑1 (MQ=255)
                                                                 tCGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGata                                >  1:6606543/1‑101 (MQ=255)
                                                                  cGAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGatat                               >  2:6606546/1‑101 (MQ=255)
                                                                   gAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATAtt                              <  2:6606014/101‑1 (MQ=255)
                                                                   gAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATAtt                              <  2:6606129/101‑1 (MQ=255)
                                                                   gAAGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATAtt                              >  1:6606549/1‑101 (MQ=255)
                                                                    aaGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCCCCCAGCGTGCAGATACGACGATATTc                             <  1:6606261/101‑1 (MQ=255)
                                                                    aaGCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCAAGCGTGCAGATACGACGATATTc                             <  1:6606313/101‑1 (MQ=255)
                                                                      gCTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGc                           >  2:6606553/1‑101 (MQ=255)
                                                                       cTGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGcc                          <  1:6606017/101‑1 (MQ=255)
                                                                        tGGCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGCCa                         <  2:6606211/101‑1 (MQ=255)
                                                                         ggCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGCCAc                        >  2:6606556/1‑101 (MQ=255)
                                                                         ggCTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGCCAc                        >  1:6606557/1‑101 (MQ=255)
                                                                           cTACCGGCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGCCACCg                      >  2:6606560/1‑101 (MQ=255)
                                                                                ggCTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGCCACCGGTGAt                 >  2:6606563/1‑101 (MQ=255)
                                                                                  cTGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGCCACCGGTGATaa               <  1:6606047/101‑1 (MQ=255)
                                                                                   tGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGCCACCGGTGATaaa              >  2:6606566/1‑101 (MQ=255)
                                                                                   tGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGCCACCGGTGATaaa              <  2:6605683/101‑1 (MQ=255)
                                                                                   tGACGACGCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGCCACCGGTGATaaa              >  1:6606567/1‑101 (MQ=255)
                                                                                     acgacgCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGCCACCGGTGATAAACt            <  1:6606182/101‑1 (MQ=255)
                                                                                        acgCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGCCACCGGTGATAAACTCGt         <  1:6606226/101‑1 (MQ=255)
                                                                                        acgCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGCCACCGGTGATAAACTCGt         >  1:6606572/1‑101 (MQ=255)
                                                                                         cgCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGCCACCGGTGATAAACTCGTg        >  2:6606573/1‑101 (MQ=255)
                                                                                          gCTCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCg                                                <  2:6605845/60‑1 (MQ=255)
                                                                                            tCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGCCACCGGTGATAAACTCGTGATc     <  2:6605883/101‑1 (MQ=255)
                                                                                            tCGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGCCACCGGTGATAAACTCGTGATc     <  2:6606055/101‑1 (MQ=255)
                                                                                             cGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGCCACCGGTGATAAACTCGTGATcc    <  1:6606102/101‑1 (MQ=255)
                                                                                             cGCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGCCACCGGTGATAAACTCGTGATcc    <  2:6605974/101‑1 (MQ=255)
                                                                                              gCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGCCACCGGTGATAAACTCGTGATCCa   >  1:6606576/1‑101 (MQ=255)
                                                                                              gCCATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGCCACCGGTGATAAACTCGTGATCCa   <  2:6606418/101‑1 (MQ=255)
                                                                                               ccATGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGCCACCGGTGATAAACTCGTGATCCAg  >  1:6606579/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
GACCTTTATAACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCACGTTCGATAAACGCATCTTCTTCCAGCAAGCCACGCAGTTTCTCACCCAGCGTGCAGATACGACGATATTCGCCACCGGTGATAAACTCGTGATCCAG  >  NC_000913/3877896‑3878091

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