Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,787,068 G→A A204T (GCC→ACC)  aroH → 3‑deoxy‑D‑arabinoheptulosonate‑7‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,787,0680GA100.0% 547.0 / NA 157A204T (GCC→ACC) aroH3‑deoxy‑D‑arabinoheptulosonate‑7‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (74/83);  total (74/83)

Only 100 of 179 total aligned reads displayed.

ATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGT  >  NC_000913/1786978‑1787157
                                                                                          |                                                                                         
aTGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCacc                                                                                         >  2:2130921/1‑93 (MQ=255)
 tGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCacca                                                                                        <  1:2213389/93‑1 (MQ=255)
   gCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCaccagc                                                                                      >  2:68437/1‑93 (MQ=255)
   gCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCaccagc                                                                                      <  1:4291882/93‑1 (MQ=255)
     ttCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCaccagcca                                                                                    >  2:3319010/1‑93 (MQ=255)
      tCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCat                                                                                   <  1:2831495/93‑1 (MQ=255)
       cGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCata                                                                                  <  2:4704175/93‑1 (MQ=255)
         gCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATg                                                                                <  1:4072624/93‑1 (MQ=255)
         gCACTCTCCTGTCCGGTAGGTTTGAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATg                                                                                <  1:2780266/93‑1 (MQ=255)
          cACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGt                                                                               >  1:1967350/1‑93 (MQ=255)
          cACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGt                                                                               <  1:3153975/93‑1 (MQ=255)
           cctctCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGtt                                                                              <  2:6097679/92‑1 (MQ=255)
            ctctCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTc                                                                             <  2:1635116/93‑1 (MQ=255)
             tctcCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTcc                                                                            <  1:2418454/93‑1 (MQ=255)
               tcCTGTCCGGTAGGTTTTAAAAATGGTACCGGTGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCctc                                                                          <  1:6202863/93‑1 (MQ=255)
               tcCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCctc                                                                          >  2:1185263/1‑93 (MQ=255)
                ccTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCctct                                                                         <  1:5336525/93‑1 (MQ=255)
                  tGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCg                                                                       >  2:2498539/1‑93 (MQ=255)
                     ccGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCa                                                                    <  2:1437616/93‑1 (MQ=255)
                     ccGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCa                                                                    <  1:4217574/93‑1 (MQ=255)
                      cGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAg                                                                   >  1:622057/1‑93 (MQ=255)
                      cGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAg                                                                   >  1:3479608/1‑93 (MQ=255)
                       ggTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGa                                                                  <  1:1725707/93‑1 (MQ=255)
                        gTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGAc                                                                 >  1:3045006/1‑93 (MQ=255)
                          aGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaa                                                               >  2:3619565/1‑93 (MQ=255)
                           ggTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaaa                                                              <  1:5055867/93‑1 (MQ=255)
                            gTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaaaa                                                             <  2:3340130/93‑1 (MQ=255)
                            gTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaaaa                                                             >  1:5763791/1‑93 (MQ=255)
                             ttttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaaaaa                                                            <  2:2432168/93‑1 (MQ=255)
                              tttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAAt                                                           <  1:3173284/93‑1 (MQ=255)
                              tttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAAt                                                           >  1:5079796/1‑93 (MQ=255)
                                tAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATgg                                                         >  1:4455834/1‑93 (MQ=255)
                                  aaaaTGTTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTc                                                       <  1:3980827/93‑1 (MQ=255)
                                  aaaaTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTc                                                       <  1:1586195/93‑1 (MQ=255)
                                  aaaaTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTc                                                       <  1:2935963/93‑1 (MQ=255)
                                   aaaTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCa                                                      <  2:4877186/93‑1 (MQ=255)
                                    aaTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAg                                                     >  2:5652436/1‑93 (MQ=255)
                                     aTGGTACCGATGGCAATAAGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGa                                                    <  2:1848010/93‑1 (MQ=255)
                                      tGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGAt                                                   >  1:3732705/1‑93 (MQ=255)
                                      tGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGAt                                                   <  2:4598914/93‑1 (MQ=255)
                                        gTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGa                                                 <  1:2653250/93‑1 (MQ=255)
                                         tACCGATGGAAATAAGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGAc                                                <  1:2086616/93‑1 (MQ=255)
                                          aCCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGAcc                                               <  2:4513473/93‑1 (MQ=255)
                                          aCCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGAcc                                               >  2:5777879/1‑93 (MQ=255)
                                           ccGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCa                                              <  2:5847544/93‑1 (MQ=255)
                                           ccGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCa                                              >  2:1672385/1‑93 (MQ=255)
                                             gagggCAATACGGGGATGGTTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATc                                            <  2:3703093/90‑1 (MQ=255)
                                             gATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATc                                            >  2:264643/1‑93 (MQ=255)
                                             gATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATc                                            >  1:544321/1‑93 (MQ=255)
                                             gATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATc                                            >  2:5793954/1‑93 (MQ=255)
                                                ggCAATACGCGGATTGTTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTAt                                         <  2:5718971/93‑1 (MQ=255)
                                                 gCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATc                                        <  1:5557603/93‑1 (MQ=255)
                                                  cAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCa                                       <  2:5752402/93‑1 (MQ=255)
                                                   aaTACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAg                                      >  2:4084442/1‑93 (MQ=255)
                                                    aTACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGa                                     >  1:1876435/1‑93 (MQ=255)
                                                     tACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGAc                                    >  2:3126616/1‑93 (MQ=255)
                                                     tACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGAc                                    <  1:1939100/93‑1 (MQ=255)
                                                        gcgGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAg                                 <  2:1698575/93‑1 (MQ=255)
                                                        gcgGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAg                                 >  2:4408098/1‑93 (MQ=255)
                                                        gcgGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAg                                 <  1:6042654/93‑1 (MQ=255)
                                                        gcgGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAg                                 >  1:1643915/1‑93 (MQ=255)
                                                         cgGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGc                                <  1:1608650/93‑1 (MQ=255)
                                                         cgGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGc                                <  1:2372807/93‑1 (MQ=255)
                                                         cgGATTGCTGTGGATGCTATCCGCGCACCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGc                                <  2:2770243/93‑1 (MQ=255)
                                                          ggATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCg                               <  2:4627839/93‑1 (MQ=255)
                                                          ggATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCg                               <  1:903853/93‑1 (MQ=255)
                                                          ggATTGCTGTGGATGCTATCCGCACAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCg                               <  1:5634961/93‑1 (MQ=255)
                                                            aTTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCTTCTATCAGACCAGCGGc                             >  1:4804673/1‑93 (MQ=255)
                                                              tGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCaa                           >  2:5685583/1‑93 (MQ=255)
                                                              tGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCaa                           <  2:354131/93‑1 (MQ=255)
                                                               gCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAAc                          <  1:3594210/93‑1 (MQ=255)
                                                                cTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAAcc                         >  1:5579938/1‑93 (MQ=255)
                                                                 tgtgGATGCTATCCGCGCAGCCCGGACCAGCCATATGGTCATCTCGGCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCAGCAAccg                        >  2:3979339/1‑92 (MQ=255)
                                                                 tgtgGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAAccc                        >  2:5894017/1‑93 (MQ=255)
                                                                 tgtgGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAAccc                        <  2:1000508/93‑1 (MQ=255)
                                                                  gtgGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCg                       >  1:5646154/1‑93 (MQ=255)
                                                                  gtgGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCg                       <  2:396545/93‑1 (MQ=255)
                                                                     gATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTAt                    <  2:4281632/93‑1 (MQ=255)
                                                                     gATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTAt                    >  1:2743442/1‑93 (MQ=255)
                                                                      aTGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATg                   <  1:1270624/93‑1 (MQ=255)
                                                                      aTGCTATCCGCCCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATg                   >  1:3816688/1‑93 (MQ=255)
                                                                         cTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGcc                <  2:4264625/93‑1 (MQ=255)
                                                                         cTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGcc                <  2:538705/93‑1 (MQ=255)
                                                                           aTCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCcac              <  2:5942084/93‑1 (MQ=255)
                                                                            tctgcgcAGCCCGCACCATCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCcaca             <  1:1587378/90‑1 (MQ=255)
                                                                            tCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCcaca             >  2:1268970/1‑93 (MQ=255)
                                                                            tCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCcaca             >  1:820653/1‑93 (MQ=255)
                                                                            cccGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCcaca             <  1:3148413/92‑1 (MQ=255)
                                                                             ccGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACat            <  1:686528/93‑1 (MQ=255)
                                                                               gcgcAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACatta          >  1:4939878/1‑93 (MQ=255)
                                                                                cgcAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACattat         <  1:1449357/93‑1 (MQ=255)
                                                                                cgcAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACattat         >  1:2569498/1‑93 (MQ=255)
                                                                                 gcAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACGATCTATCAGACCAGCGGCAACCCGTATGGCCACattatt        <  1:4550111/93‑1 (MQ=255)
                                                                                  cAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGGCCATCTATCAGACCAGCGGCAACCCGTATGGCCACAttatta       <  2:2254957/93‑1 (MQ=255)
                                                                                  cAGCCCCCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACAttatta       <  2:1366217/93‑1 (MQ=255)
                                                                                   aGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATtattat      <  1:1224667/93‑1 (MQ=255)
                                                                                   aGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATtattat      >  2:884927/1‑93 (MQ=255)
                                                                                    acccGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATg     <  1:6115339/92‑1 (MQ=255)
                                                                                      ccGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCg   >  2:108121/1‑93 (MQ=255)
                                                                                       cgcaCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGt  <  1:1548996/93‑1 (MQ=255)
                                                                                          |                                                                                         
ATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGT  >  NC_000913/1786978‑1787157

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