Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,138,339 G→A V247I (GTT→ATT)  flgK → hook‑filament junction protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,138,3390GA100.0% 446.5 / NA 126V247I (GTT→ATT) flgKhook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (60/66);  total (60/66)

Only 100 of 131 total aligned reads displayed.

GTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAAT  >  NC_000913/1138249‑1138428
                                                                                          |                                                                                         
gTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAAtt                                                                                         >  2:4080008/1‑93 (MQ=255)
 tCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGCACGGCGCGGCAACTGGCGGCAAttc                                                                                        >  1:4528369/1‑93 (MQ=255)
  cAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAAttcc                                                                                       >  2:2926190/1‑93 (MQ=255)
  cAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAAttcc                                                                                       >  2:3499767/1‑93 (MQ=255)
    gCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAAttcctt                                                                                     >  2:4121788/1‑93 (MQ=255)
     cGTTCAGGATGGCGTCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAAttccttc                                                                                    <  1:1668805/93‑1 (MQ=255)
      gTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAAttccttcc                                                                                   >  1:416506/1‑93 (MQ=255)
      gTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAAttccttcc                                                                                   >  1:3721120/1‑93 (MQ=255)
      gTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAAttccttcc                                                                                   >  1:1984123/1‑93 (MQ=255)
      gTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAAttccttcc                                                                                   <  2:2330747/93‑1 (MQ=255)
       ttCAGGATGGCTGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCa                                                                                  <  2:280880/93‑1 (MQ=255)
       ttCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCa                                                                                  >  2:1281409/1‑93 (MQ=255)
         cAGGATGGTGGCACTTATAACAGCGCGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAgc                                                                                <  2:331265/93‑1 (MQ=255)
         cAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAgc                                                                                <  2:304023/93‑1 (MQ=255)
         cAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAgc                                                                                <  1:430402/93‑1 (MQ=255)
             aTGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTg                                                                            >  2:4523879/1‑93 (MQ=255)
              tGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGa                                                                           >  2:3909169/1‑93 (MQ=255)
              tGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGa                                                                           <  2:3889905/93‑1 (MQ=255)
              tGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGa                                                                           <  1:2123934/93‑1 (MQ=255)
              tGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGa                                                                           >  1:1967370/1‑93 (MQ=255)
              tGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGa                                                                           >  2:1929281/1‑93 (MQ=255)
               ggcggcACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGAc                                                                          >  2:1997598/1‑93 (MQ=255)
                gcggcACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGAcc                                                                         <  2:245505/93‑1 (MQ=255)
                    cACTTATAACATCACGTTGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTtc                                                                     <  2:664781/93‑1 (MQ=255)
                     aCTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTtct                                                                    >  2:3848405/1‑93 (MQ=255)
                     aCTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTtct                                                                    <  1:1557337/93‑1 (MQ=255)
                     aCTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTtct                                                                    <  2:186837/93‑1 (MQ=255)
                      cTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTtctc                                                                   >  1:2492943/1‑93 (MQ=255)
                      cTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTtctc                                                                   >  1:1978104/1‑93 (MQ=255)
                      cTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTtctc                                                                   >  1:4216547/1‑93 (MQ=255)
                       ttATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCg                                                                  >  2:473424/1‑93 (MQ=255)
                        tataACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGt                                                                 >  1:1386445/1‑93 (MQ=255)
                          taACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTAc                                                               >  2:3106255/1‑93 (MQ=255)
                          taACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTAc                                                               <  2:3615060/93‑1 (MQ=255)
                            aCATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCGTCTCGTACGa                                                             >  2:1881792/1‑93 (MQ=255)
                            aCATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCCGACCCTTTTAGTACgc                                                             >  2:28390/1‑92 (MQ=255)
                             cATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGAc                                                            <  1:3930601/93‑1 (MQ=255)
                                cACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGt                                                         <  2:3032270/93‑1 (MQ=255)
                                     tGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCtt                                                    >  1:4527976/1‑93 (MQ=255)
                                     tGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTGCTCGTACGACTGTCGCtt                                                    <  1:784517/93‑1 (MQ=255)
                                      ggCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTa                                                   >  2:3288181/1‑93 (MQ=255)
                                      ggCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTa                                                   <  2:3005986/93‑1 (MQ=255)
                                        ccAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATg                                                 <  2:2073895/93‑1 (MQ=255)
                                        ccAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATg                                                 <  1:1452583/93‑1 (MQ=255)
                                         cAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGt                                                >  1:954478/1‑93 (MQ=255)
                                         cAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGt                                                <  1:3082423/93‑1 (MQ=255)
                                          aaTGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGtt                                               >  1:1773260/1‑93 (MQ=255)
                                          aaTGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGtt                                               <  2:123125/93‑1 (MQ=255)
                                           aTGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTg                                              >  2:3869891/1‑93 (MQ=255)
                                           aTGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTg                                              >  1:410704/1‑93 (MQ=255)
                                            tGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTgt                                             >  2:3880194/1‑92 (MQ=255)
                                            tGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGa                                             <  1:1528175/93‑1 (MQ=255)
                                              gTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATg                                           <  2:905/93‑1 (MQ=255)
                                              gTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATg                                           >  1:1653300/1‑93 (MQ=255)
                                              gTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATg                                           >  1:4538451/1‑93 (MQ=255)
                                               ttACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATgg                                          <  1:3935705/93‑1 (MQ=255)
                                               ttACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATgg                                          <  2:1020307/93‑1 (MQ=255)
                                                tACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATggg                                         >  2:157058/1‑93 (MQ=255)
                                                tACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATggg                                         >  1:1515598/1‑93 (MQ=255)
                                                 aCTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGa                                        <  1:4461321/93‑1 (MQ=255)
                                                 aCTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGa                                        >  2:2237201/1‑93 (MQ=255)
                                                 aCTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGa                                        <  1:4376029/93‑1 (MQ=255)
                                                  cTCACTGGTTCAGTGAAGTACGGCGCGTCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGAc                                       <  1:3033153/93‑1 (MQ=255)
                                                  cTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGAc                                       <  2:4447485/93‑1 (MQ=255)
                                                   tCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGGTTATGTTGATGGGatg                                      >  2:3614321/1‑91 (MQ=255)
                                                    cACTGGTTCAGGGAAGTACGGCGTGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACgg                                     <  2:2701515/93‑1 (MQ=255)
                                                     aCTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACggc                                    <  2:376688/93‑1 (MQ=255)
                                                      cTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACggca                                   <  1:4055736/93‑1 (MQ=255)
                                                      cTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACggca                                   >  2:1787239/1‑93 (MQ=255)
                                                       tGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACggcag                                  <  1:274420/93‑1 (MQ=255)
                                                        ggTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACggcagg                                 <  1:2268217/93‑1 (MQ=255)
                                                             aGGGAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAata                            <  1:4411532/93‑1 (MQ=255)
                                                              gggAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAatat                           <  2:1873341/93‑1 (MQ=255)
                                                              gggAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAatat                           <  1:2726361/93‑1 (MQ=255)
                                                              gggAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAatat                           <  2:2483391/93‑1 (MQ=255)
                                                              gggAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAatat                           <  2:1277264/93‑1 (MQ=255)
                                                               ggAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATAtt                          <  2:1445961/93‑1 (MQ=255)
                                                                gAAGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTg                         <  1:2828551/93‑1 (MQ=255)
                                                                  aGTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTgag                       >  2:3996591/1‑93 (MQ=255)
                                                                   gTACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTgaga                      <  1:3685002/93‑1 (MQ=255)
                                                                    tACGGCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGAt                     <  2:4389256/93‑1 (MQ=255)
                                                                       ggCGCGGCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATccc                  >  1:2326386/1‑93 (MQ=255)
                                                                            ggCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGACGGAAAAATTGAGATCCCGggga             >  2:3646814/1‑93 (MQ=255)
                                                                             gCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGaa            <  1:2898997/93‑1 (MQ=255)
                                                                             gCAACTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGaa            <  1:2392557/93‑1 (MQ=255)
                                                                                 cTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTa        <  1:3089155/93‑1 (MQ=255)
                                                                                 cTGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTa        >  1:366050/1‑93 (MQ=255)
                                                                                  tGGCGGCAATTCCTTCCAGCTCTGACCCTTCGCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTAc       <  2:790494/93‑1 (MQ=255)
                                                                                  tGGCGGCAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTAc       <  1:1859110/93‑1 (MQ=255)
                                                                                   ggcggcAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACt      <  1:1681374/93‑1 (MQ=255)
                                                                                   ggcgccAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACt      <  2:1654733/93‑1 (MQ=255)
                                                                                    gcggcAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTg     >  1:1231573/1‑93 (MQ=255)
                                                                                     cggcAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGa    >  2:3799182/1‑93 (MQ=255)
                                                                                      ggcAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGaa   >  2:4006398/1‑93 (MQ=255)
                                                                                      ggcAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGaa   >  1:1455416/1‑93 (MQ=255)
                                                                                      ggcAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGaa   <  2:2006278/93‑1 (MQ=255)
                                                                                       gcAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAAt  >  1:2755959/1‑93 (MQ=255)
                                                                                       gcAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAAt  >  1:3178416/1‑93 (MQ=255)
                                                                                       gcAATTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAAt  <  1:3454484/93‑1 (MQ=255)
                                                                                       gcAATTCCGTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAAt  >  2:4581675/1‑93 (MQ=255)
                                                                                          |                                                                                         
GTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAAT  >  NC_000913/1138249‑1138428

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