Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,388,932 A→G A52A (GCT→GCC b4161 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,388,9320AG100.0% 936.3 / NA 290A52A (GCT→GCCb4161hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (109/181);  total (109/181)

Only 100 of 307 total aligned reads displayed.

ACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTT  >  NC_000913/4388842‑4389021
                                                                                          |                                                                                         
  cccAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGt                                                                                       <  1:3314017/92‑1 (MQ=255)
   ccAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTg                                                                                      >  1:5852640/1‑93 (MQ=255)
     aGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCa                                                                                    >  1:803325/1‑93 (MQ=255)
      gACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCAt                                                                                   <  2:5361530/93‑1 (MQ=255)
      gACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCAt                                                                                   >  2:6080254/1‑93 (MQ=255)
      gACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCAt                                                                                   <  1:4765550/93‑1 (MQ=255)
        cTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATAc                                                                                 <  1:6384328/93‑1 (MQ=255)
         tACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATAcc                                                                                >  2:3346524/1‑93 (MQ=255)
         tACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATAcc                                                                                >  1:4282063/1‑93 (MQ=255)
          aCGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCa                                                                               >  1:5436161/1‑93 (MQ=255)
          aCGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGCGCATACCa                                                                               >  2:5802601/1‑93 (MQ=255)
          aCGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGAGCATACCa                                                                               >  2:3300369/1‑93 (MQ=255)
           cgcgGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCaa                                                                              >  2:5605229/1‑93 (MQ=255)
              ggTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGc                                                                           >  1:3234099/1‑93 (MQ=255)
              ggTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGc                                                                           >  1:3455062/1‑93 (MQ=255)
               gTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCg                                                                          >  1:5721215/1‑93 (MQ=255)
                  gCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCt                                                                       <  1:2161189/93‑1 (MQ=255)
                  gCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCt                                                                       <  1:3921690/93‑1 (MQ=255)
                    cGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCCAAGCGGAtga                                                                     >  1:5077941/1‑93 (MQ=255)
                     ggTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgat                                                                    >  1:1421878/1‑93 (MQ=255)
                      gTTACCAGCGAACGGATGGTACGGCTAATAGTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatg                                                                   <  1:5548864/93‑1 (MQ=255)
                      gTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatg                                                                   <  1:4446439/93‑1 (MQ=255)
                      gTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatg                                                                   <  1:1042786/93‑1 (MQ=255)
                      gTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatg                                                                   >  2:3943438/1‑93 (MQ=255)
                      gTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatg                                                                   <  2:6239425/93‑1 (MQ=255)
                       ttACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatga                                                                  <  1:3361894/93‑1 (MQ=255)
                       ttACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatga                                                                  <  2:2594358/93‑1 (MQ=255)
                        tACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatggg                                                                 >  1:3675839/1‑91 (MQ=255)
                        tACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGAc                                                                 <  2:3036069/93‑1 (MQ=255)
                           cAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGAt                                                              <  1:960527/93‑1 (MQ=255)
                             gCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATAc                                                            <  1:6490579/93‑1 (MQ=255)
                             gCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATAc                                                            <  1:1314917/93‑1 (MQ=255)
                              cGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATAcc                                                           <  2:2050818/93‑1 (MQ=255)
                               gagcGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCt                                                          <  1:604683/90‑1 (MQ=255)
                               gAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCt                                                          <  2:2886838/93‑1 (MQ=255)
                                   ggATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcat                                                      >  1:6610203/1‑93 (MQ=255)
                                   ggATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcat                                                      >  2:1718167/1‑93 (MQ=255)
                                   ggATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcat                                                      >  1:4883514/1‑93 (MQ=255)
                                   ggATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcat                                                      <  2:5225277/93‑1 (MQ=255)
                                     aTGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcatca                                                    <  2:2627636/93‑1 (MQ=255)
                                      tGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAg                                                   >  2:3109411/1‑93 (MQ=255)
                                      tGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAg                                                   <  2:50548/93‑1 (MQ=255)
                                      tGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAg                                                   >  1:2405896/1‑93 (MQ=255)
                                      tGGTACGGCGAATATTGCAGCGGTGAACGTCGCAATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAg                                                   <  1:2933643/93‑1 (MQ=255)
                                        gTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGc                                                 <  1:2108333/93‑1 (MQ=255)
                                        gTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGc                                                 <  2:1184799/93‑1 (MQ=255)
                                        gTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGc                                                 <  2:6375612/93‑1 (MQ=255)
                                          aCGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTc                                               <  2:3162716/93‑1 (MQ=255)
                                          aCGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTc                                               <  2:6688542/93‑1 (MQ=255)
                                           cGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCg                                              >  1:4943734/1‑93 (MQ=255)
                                            ggCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTctc                                             <  1:2601263/93‑3 (MQ=255)
                                            ggCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGc                                             <  2:81158/93‑1 (MQ=255)
                                            ggCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGc                                             >  1:4570390/1‑93 (MQ=255)
                                               gAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCaa                                          >  2:6511318/1‑93 (MQ=255)
                                               gAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCCGATGACGATACCTTCATCAGGCTCGccca                                          >  1:4143181/1‑91 (MQ=255)
                                                  tatTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACa                                       >  2:147204/1‑93 (MQ=255)
                                                  tatTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACa                                       >  2:4174295/1‑93 (MQ=255)
                                                    ttGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGa                                     <  1:2180202/93‑1 (MQ=255)
                                                    ttGCAGCGGTGAACGACGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGa                                     <  2:4999848/93‑1 (MQ=255)
                                                     tGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGAt                                    <  2:5414773/93‑1 (MQ=255)
                                                     tGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGAt                                    <  2:2343248/93‑1 (MQ=255)
                                                     tGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGAt                                    <  2:1279836/93‑1 (MQ=255)
                                                     tGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGAt                                    <  1:1916757/93‑1 (MQ=255)
                                                      gCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGAtt                                   <  2:5923976/93‑1 (MQ=255)
                                                      gCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGAtt                                   <  2:2812412/93‑1 (MQ=255)
                                                       cAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTg                                  <  2:2742698/93‑1 (MQ=255)
                                                          cGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtcg                               >  2:1544359/1‑93 (MQ=255)
                                                           cgTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtcgt                              <  2:4878705/92‑1 (MQ=255)
                                                             tGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGtcgtcg                            >  1:5671289/1‑93 (MQ=255)
                                                             tGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGtcgtcg                            <  2:3668218/93‑1 (MQ=255)
                                                              gAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTcgtcgt                           >  2:6594294/1‑93 (MQ=255)
                                                                 cGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGt                        <  1:6216957/93‑1 (MQ=255)
                                                                 cGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGt                        <  2:5118963/93‑1 (MQ=255)
                                                                 cGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGt                        >  2:1059937/1‑93 (MQ=255)
                                                                   tCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCg                      <  1:72919/93‑1 (MQ=255)
                                                                   tCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACACATTGTCGTCGTAGTCg                      >  2:57093/1‑93 (MQ=255)
                                                                     gCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCggg                    <  2:2226953/93‑1 (MQ=255)
                                                                      ccATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGc                   <  2:2720218/93‑1 (MQ=255)
                                                                      ccATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGc                   <  1:300666/93‑1 (MQ=255)
                                                                       cATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCt                  <  2:4918787/93‑1 (MQ=255)
                                                                         tCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTtc                <  1:27205/93‑1 (MQ=255)
                                                                           ggcggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTtctc              >  2:1913741/1‑93 (MQ=255)
                                                                           ggcggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTtctc              <  1:232046/93‑1 (MQ=255)
                                                                            gcggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTcc             <  2:4507229/93‑1 (MQ=255)
                                                                            gcggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTcc             <  2:1784858/93‑1 (MQ=255)
                                                                             cggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCt            <  1:3191634/93‑1 (MQ=255)
                                                                             cggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCt            <  2:432197/93‑1 (MQ=255)
                                                                             cggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCt            <  2:350262/93‑1 (MQ=255)
                                                                              ggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCtt           >  1:2521459/1‑93 (MQ=255)
                                                                              ggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCtt           <  1:1666584/93‑1 (MQ=255)
                                                                               gATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTa          <  2:6017555/93‑1 (MQ=255)
                                                                               gATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTa          <  2:6634202/93‑1 (MQ=255)
                                                                                 ttCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGa        >  2:6475395/1‑93 (MQ=255)
                                                                                    cacaTCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGt     <  2:135561/93‑1 (MQ=255)
                                                                                    cacaTCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGt     <  2:1875195/93‑1 (MQ=255)
                                                                                    cacaTCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGt     >  1:5295048/1‑93 (MQ=255)
                                                                                    cacaTCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGt     >  1:1869436/1‑93 (MQ=255)
                                                                                       aTCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGtttt  >  2:46526/1‑93 (MQ=255)
                                                                                       aTCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGtttt  <  1:3902552/93‑1 (MQ=255)
                                                                                       aTCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGtttt  <  1:1232815/93‑1 (MQ=255)
                                                                                          |                                                                                         
ACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTT  >  NC_000913/4388842‑4389021

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