Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,566,875 A→G R29R (CGT→CGC glgC ← glucose‑1‑phosphate adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,566,8750AG100.0% 453.4 / NA 149R29R (CGT→CGCglgCglucose‑1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (105/44);  total (105/44)

Only 100 of 236 total aligned reads displayed.

AGACAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCACGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGAC  >  NC_000913/3566785‑3566964
                                                                                          |                                                                                         
 gACAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGt                                                                                        >  1:2805069/1‑93 (MQ=255)
 gACAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGt                                                                                        >  2:3421126/1‑93 (MQ=255)
  aCAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGtc                                                                                       >  2:2416133/1‑93 (MQ=255)
  aCAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATACTTCAGGCTGGGACCGCGtc                                                                                       >  2:4412468/1‑93 (MQ=255)
   cAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGtcc                                                                                      <  2:4196913/93‑1 (MQ=255)
     gcgcAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGtcctc                                                                                    <  1:2956708/93‑1 (MQ=255)
       gcAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTccc                                                                                  >  1:4488230/1‑93 (MQ=255)
       gcAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTccc                                                                                  <  2:704006/93‑1 (MQ=255)
        cAAAGTCGCTAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCg                                                                                 >  1:3471716/1‑93 (MQ=255)
        cAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCg                                                                                 >  2:1223264/1‑93 (MQ=255)
           aGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCa                                                                              <  1:4536584/93‑1 (MQ=255)
            gTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAg                                                                             >  1:361562/1‑93 (MQ=255)
            gTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAg                                                                             >  1:2828358/1‑93 (MQ=255)
              cGATAATGCGGAACTTCCCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTa                                                                           <  1:3281447/93‑1 (MQ=255)
              cGATAATGCGGAACTTCCCGCCAAATTGTCGGGCCGGTTTTGCTCGCTAATTGGTGAAACCCTCCAGGCGGTTACCGCGTCCTCCCGCCAGTa                                                                           <  2:1181377/93‑1 (MQ=255)
              cGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGTTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTa                                                                           <  2:1749898/93‑1 (MQ=255)
                 tAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCa                                                                        >  1:1180715/1‑93 (MQ=255)
                 aaaTGCGGAACTTCCCGCCGAAGTGTCCGGCCGGTTTTGCTCGCTTATTGGTTAAACCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCa                                                                        <  2:3392992/92‑1 (MQ=255)
                  aaTGGGGAACTTACCGCCGAAGTGTACGGCCGTTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAg                                                                       <  1:3051249/93‑1 (MQ=255)
                   aTGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAgg                                                                      <  1:607657/93‑1 (MQ=255)
                   aTGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAgg                                                                      >  2:4277304/1‑93 (MQ=255)
                   aTGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAgg                                                                      >  2:3622056/1‑93 (MQ=255)
                    tGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAggg                                                                     <  1:809001/93‑1 (MQ=255)
                     gCGGAACTTACCGCCGAAGTGTACGGCCGGGTTTGCTCGCTTATTGGTTACATCCTTCAGGCGGGTACCGCGTGCTCCCGCCCGTATCAGGGc                                                                    >  2:2706849/1‑93 (MQ=255)
                      cGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCa                                                                   >  1:716990/1‑93 (MQ=255)
                       ggAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACAGCGTCCTCCCGCCAGTATCAGGGCaa                                                                  >  1:2364870/1‑93 (MQ=255)
                        gaccTTCCCGCCGAGGTGTACGGCTGGTTTTGCGCGCTTATTGGTTAAATCCTTCAGGCGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAAc                                                                 <  1:3525513/90‑1 (MQ=255)
                        gAACTTCCCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAAc                                                                 <  1:2104408/93‑1 (MQ=255)
                        gAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAAc                                                                 >  1:1062508/1‑93 (MQ=255)
                          aCTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAg                                                               >  2:2841983/1‑93 (MQ=255)
                          aCTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAg                                                               >  1:2205830/1‑93 (MQ=255)
                           cTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGa                                                              <  2:1613299/93‑1 (MQ=255)
                            ttACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGAt                                                             <  1:1796351/93‑1 (MQ=255)
                             tACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGAtt                                                            >  2:1536650/1‑93 (MQ=255)
                             tACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGAtt                                                            >  1:4034563/1‑93 (MQ=255)
                              aCCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGAttt                                                           >  1:2120688/1‑93 (MQ=255)
                               cccccgAAGTGTCCGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTc                                                          <  1:3373030/90‑1 (MQ=255)
                                 gccgAATTGTACGGCCGGTTTTGCTCGTTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCaa                                                        <  2:2136327/93‑1 (MQ=255)
                                 gccgAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCaa                                                        >  1:2057481/1‑93 (MQ=255)
                                   cgAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATg                                                      >  2:721781/1‑93 (MQ=255)
                                   cgAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATg                                                      <  1:4289201/93‑1 (MQ=255)
                                    gAAGTGTACGGCCGGTTTTGCTCGCTTTTTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATgg                                                     >  2:2234374/1‑93 (MQ=255)
                                    gAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATgg                                                     >  2:4665276/1‑93 (MQ=255)
                                         gtACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCt                                                >  1:1262921/1‑93 (MQ=255)
                                         gtACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCt                                                <  1:1441610/93‑1 (MQ=255)
                                         gtACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCt                                                >  1:2295494/1‑93 (MQ=255)
                                         gtACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCt                                                >  2:1710960/1‑93 (MQ=255)
                                          tACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTg                                               >  1:1202880/1‑93 (MQ=255)
                                          tACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTg                                               >  2:3763596/1‑93 (MQ=255)
                                          tACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTg                                               >  1:467550/1‑93 (MQ=255)
                                           aCGGCCGTTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTgg                                              <  2:2604004/93‑1 (MQ=255)
                                             ggCCGGGTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGGCCTCCCGCCAGTATCAGGGCGACAGATTTCAATGGCAGCTgggg                                            >  2:1766971/1‑91 (MQ=255)
                                              gCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGgcgc                                           >  1:1954782/1‑93 (MQ=255)
                                                 ggTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGcc                                        <  2:3073727/93‑1 (MQ=255)
                                                  gTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCa                                       >  1:3812650/1‑93 (MQ=255)
                                                  gTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCa                                       >  2:1439344/1‑93 (MQ=255)
                                                    tttGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAAc                                     >  1:2740492/1‑93 (MQ=255)
                                                      tGCTCTCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACAt                                   <  1:3387294/93‑1 (MQ=255)
                                                      tGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACAt                                   >  2:239615/1‑93 (MQ=255)
                                                      tGCTCGCTTATTGGTAAAATCCTTCAGGCGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACAt                                   <  1:3843743/93‑1 (MQ=255)
                                                       gCTCTCTTATTGTTTAAATCCTACGGGCGAGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACAtt                                  <  2:1075970/93‑1 (MQ=255)
                                                         tCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTaa                                >  1:3798437/1‑93 (MQ=255)
                                                         tCGCTAATTGGTAAAATCCTTCAGGCGTGCACCGCGTCCCCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTaa                                <  1:2541772/93‑1 (MQ=255)
                                                          ctcTTTTTGTTTAAATCCTTCAGGCGGGTACCCCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAg                               <  1:912884/91‑1 (MQ=255)
                                                          cGCTTTTTGGTAAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAg                               <  2:3485659/93‑1 (MQ=255)
                                                           tcTTATTGTTAAAATCCTTCAGGCGGTTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGt                              <  2:2496821/92‑1 (MQ=255)
                                                            cTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTg                             >  2:4119606/1‑93 (MQ=255)
                                                            cTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCCCCAACATCAAGTg                             >  2:2450225/1‑93 (MQ=255)
                                                              ttttGATTAAATCCTTCGGGCGGTCCCCGCGCCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGAt                           <  1:322305/91‑1 (MQ=255)
                                                              tATTGTTTAAATCCTTCAGGGGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACAATAAGTGAt                           <  2:4571278/93‑1 (MQ=255)
                                                              tATTGTTTAAATCCTCCGGGCGGGTCCCGCGCCCCCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGAt                           <  1:1628866/93‑1 (MQ=255)
                                                              tATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGAt                           >  2:3003321/1‑93 (MQ=255)
                                                               tttGGTTAAATCCTTCGGGGGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATc                          <  1:3769727/92‑1 (MQ=255)
                                                               atgggTAAAATCTTTCAGGCGGCTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATc                          <  1:1315888/90‑1 (MQ=255)
                                                               aTTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATc                          >  1:2422930/1‑93 (MQ=255)
                                                                 tGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGGGATCGt                        >  2:2327825/1‑93 (MQ=255)
                                                                  ggTTAAATCCTTCAGGCGGGTACCGCCCCCTACCGCCAGTATCAGGGCAACAGCTTTCAATGGCAGCTGGCGCGCCAACCCTTAGTGATCgat                       >  2:1901373/1‑91 (MQ=255)
                                                                   gTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGttc                      <  2:3377248/93‑1 (MQ=255)
                                                                    ttAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGttct                     >  1:4072837/1‑93 (MQ=255)
                                                                    taaaaCCCTGCAGGCGGTTCCCCCGCCCGCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGttct                     <  2:2597010/91‑1 (MQ=255)
                                                                     tAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGttctt                    >  1:98597/1‑93 (MQ=255)
                                                                      aaaTCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGttcttc                   >  2:2647180/1‑93 (MQ=255)
                                                                      aaaCCCTCCAGGCGGGTGCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGttcttc                   <  1:420820/93‑1 (MQ=255)
                                                                       aaTCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATCGCAGCTGTCGCGCCAACATTAAGTGATCGTtcttct                  >  1:4509020/1‑93 (MQ=255)
                                                                            ctCAGGCGGTTCCCCCGCCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTaaa             <  2:1940033/92‑1 (MQ=255)
                                                                             tCAGGCGGGTCCCGCGCCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAAc            <  1:691844/93‑1 (MQ=255)
                                                                              cAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACAATAAGTGATCGGTCTTCTCTAAACt           >  1:3938251/1‑93 (MQ=255)
                                                                               aGGGGGGTGCCGCGCCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTa          <  2:4060255/93‑1 (MQ=255)
                                                                               aGGCGGGTACCGCGCCCCCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTa          <  2:1357036/93‑1 (MQ=255)
                                                                                ggCGGGTACCGCGTCCTCCCGTCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTaa         >  2:3931282/1‑93 (MQ=255)
                                                                                ggCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTaa         >  2:3302456/1‑93 (MQ=255)
                                                                                 gCGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAAc        <  2:2868285/93‑1 (MQ=255)
                                                                                 gCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAAc        >  2:147279/1‑93 (MQ=255)
                                                                                   gggTCCCGCGCCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCa      <  2:1689812/93‑1 (MQ=255)
                                                                                   gggTCCCGCGCCCCCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCa      <  2:3130595/93‑1 (MQ=255)
                                                                                   gggTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCa      <  1:540445/93‑1 (MQ=255)
                                                                                    ggTCCCCCGCCCCCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCAt     <  2:3588435/93‑1 (MQ=255)
                                                                                      tACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGa   >  2:424121/1‑93 (MQ=255)
                                                                                       aCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGAc  >  1:2936306/1‑93 (MQ=255)
                                                                                       aCCGCGCCCACCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGAc  <  2:4054707/93‑1 (MQ=255)
                                                                                          |                                                                                         
AGACAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCACGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGAC  >  NC_000913/3566785‑3566964

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