Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,132,875 T→C L360P (CTG→CCG)  flgE → hook protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,132,8750TC100.0% 307.5 / NA 100L360P (CTG→CCG) flgEhook protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (35/65);  total (35/65)

Only 100 of 103 total aligned reads displayed.

CAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGACCATCAAAA  >  NC_000913/1132786‑1132964
                                                                                         |                                                                                         
cAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGa                                                                                        >  1:2021575/1‑93 (MQ=255)
  aCGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAg                                                                                      <  1:3053411/93‑1 (MQ=255)
  aCGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAg                                                                                      <  2:4193845/93‑1 (MQ=255)
  aCGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAg                                                                                      >  2:4300094/1‑93 (MQ=255)
  aCGTCTGGTCTGCAACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAg                                                                                      >  1:1978547/1‑93 (MQ=255)
    gtctggtctgCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCg                                                                                    <  1:3771544/93‑1 (MQ=255)
     tctggtctgCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGt                                                                                   <  2:1162023/93‑1 (MQ=255)
       tggtctgCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTcc                                                                                 >  2:2158133/1‑93 (MQ=255)
            tgCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGt                                                                            >  2:574205/1‑93 (MQ=255)
             gCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTc                                                                           >  2:3667910/1‑93 (MQ=255)
             gCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTc                                                                           >  2:2115718/1‑93 (MQ=255)
               gACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGa                                                                         >  2:2644/1‑93 (MQ=255)
               gACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGa                                                                         <  2:626363/93‑1 (MQ=255)
                aCGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGAt                                                                        <  2:2230416/93‑1 (MQ=255)
                aCGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGAt                                                                        >  2:4023679/1‑93 (MQ=255)
                 cGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGAtc                                                                       <  1:3012049/93‑1 (MQ=255)
                 cGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGAtc                                                                       <  1:2968006/93‑1 (MQ=255)
                  gCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGAtct                                                                      <  2:4054340/93‑1 (MQ=255)
                    aaTCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCa                                                                    <  2:560566/93‑1 (MQ=255)
                    aaTCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCa                                                                    <  2:3588893/93‑1 (MQ=255)
                    aaTCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCa                                                                    <  2:1692864/93‑1 (MQ=255)
                      tcttctGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGt                                                                  <  1:4304298/93‑1 (MQ=255)
                       cttctGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAgga                                                                 >  2:510448/1‑91 (MQ=255)
                       cttctGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTa                                                                 <  1:1316249/93‑1 (MQ=255)
                        ttctGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTaa                                                                >  1:79008/1‑93 (MQ=255)
                        ttctGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTaa                                                                >  1:1071401/1‑93 (MQ=255)
                         tctGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTaaa                                                               <  2:3431737/93‑1 (MQ=255)
                          ctGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAg                                                              <  1:1382059/93‑1 (MQ=255)
                           tggcgtggcgCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGa                                                             <  1:193707/93‑1 (MQ=255)
                            ggcgtggcgCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGaa                                                            <  1:383206/93‑1 (MQ=255)
                            ggcgtggcgCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGaa                                                            >  1:1285506/1‑93 (MQ=255)
                            ggcgtggcgCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGaa                                                            <  1:3339341/93‑1 (MQ=255)
                             gcgtggcgCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAAc                                                           <  2:1808064/93‑1 (MQ=255)
                             gcgtggcgCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAAc                                                           <  2:3658310/93‑1 (MQ=255)
                             gcgtggcgCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAAc                                                           <  2:3626025/93‑1 (MQ=255)
                               gtggcgCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTg                                                         >  1:1057977/1‑93 (MQ=255)
                               gtggcgCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTg                                                         <  1:4009265/93‑1 (MQ=255)
                                tggcgCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTgg                                                        <  1:3162459/93‑1 (MQ=255)
                                tggcgCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTgg                                                        <  2:1107965/93‑1 (MQ=255)
                                 ggcgCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGt                                                       >  1:3216681/1‑93 (MQ=255)
                                  gcgcTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTc                                                      <  1:2993372/93‑1 (MQ=255)
                                  gcgcTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTc                                                      <  1:3652570/93‑1 (MQ=255)
                                    gcTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCaa                                                    <  1:3952915/93‑1 (MQ=255)
                                     cTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAat                                                   >  1:3588285/1‑93 (MQ=255)
                                     cTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAat                                                   <  1:1671140/93‑1 (MQ=255)
                                      tGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAata                                                  >  1:4676957/1‑93 (MQ=255)
                                       gTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAatat                                                 >  1:2253377/1‑93 (MQ=255)
                                        gtGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATg                                                <  1:379909/92‑1 (MQ=255)
                                         tGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGa                                               <  2:4075657/93‑1 (MQ=255)
                                         tGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGa                                               <  1:4356931/93‑1 (MQ=255)
                                          ggggACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGAt                                              <  1:1377048/93‑1 (MQ=255)
                                          ggggACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACATCGATCTCAGTAAAGAACTGGTCAATATGAt                                              >  2:2936395/1‑93 (MQ=255)
                                           gggACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATc                                             <  2:4440205/93‑1 (MQ=255)
                                           gggACAGCCGGCACGGAAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATc                                             <  2:2422081/93‑1 (MQ=255)
                                             gACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGt                                           <  2:3408801/93‑1 (MQ=255)
                                              aCAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGtt                                          <  2:3957922/93‑1 (MQ=255)
                                              aCAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGtt                                          <  2:4488654/93‑1 (MQ=255)
                                                  ccGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGccc                                      <  1:3459954/93‑1 (MQ=255)
                                                   cgggacgggaAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCa                                     <  2:452870/93‑1 (MQ=255)
                                                    gggacggtaAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAg                                    >  2:4609534/1‑93 (MQ=255)
                                                    gggacgggaAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAg                                    >  2:1986062/1‑93 (MQ=255)
                                                    gggacgggaAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAg                                    >  1:4518041/1‑93 (MQ=255)
                                                    gggacgggaAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAg                                    >  2:3523913/1‑93 (MQ=255)
                                                        cgggaAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTa                                <  2:3056371/93‑1 (MQ=255)
                                                         gggaAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTaa                               <  1:702877/93‑1 (MQ=255)
                                                         gggaAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTaa                               <  1:4507570/93‑1 (MQ=255)
                                                          ggaAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAAc                              <  1:1677729/93‑1 (MQ=255)
                                                          ggaAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAAc                              <  2:772707/93‑1 (MQ=255)
                                                          ggaAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAAc                              >  1:2361294/1‑93 (MQ=255)
                                                          ggaAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAAc                              <  2:3920952/93‑1 (MQ=255)
                                                          ggaAACTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAAc                              <  2:3957802/93‑1 (MQ=255)
                                                            aaaCTTTCGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTa                            <  1:4294869/93‑1 (MQ=255)
                                                             aaCTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTAt                           <  1:102508/93‑1 (MQ=255)
                                                             aaCTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTAt                           <  2:456691/93‑1 (MQ=255)
                                                              aCTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATc                          <  1:4331416/93‑1 (MQ=255)
                                                              aCTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATc                          >  2:1425738/1‑93 (MQ=255)
                                                               cTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCa                         <  1:2708084/93‑1 (MQ=255)
                                                               cTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCa                         <  2:262851/93‑1 (MQ=255)
                                                               cTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCa                         <  2:573437/93‑1 (MQ=255)
                                                               cTTTGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCa                         <  1:18697/93‑1 (MQ=255)
                                                                tttGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAg                        >  1:1926188/1‑93 (MQ=255)
                                                                tttGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAg                        >  2:3563534/1‑93 (MQ=255)
                                                                tttGGCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAg                        <  2:1621228/93‑1 (MQ=255)
                                                                 ttGGCACCCTGACCAACGGTGCTCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGt                       <  1:2388356/93‑1 (MQ=255)
                                                                   ggCACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCt                     >  2:2109102/1‑93 (MQ=255)
                                                                     cACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTaa                   <  1:1397384/93‑1 (MQ=255)
                                                                     cACCCTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTAACAGTCTaa                   <  1:1947722/93‑1 (MQ=255)
                                                                        ccTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGc                >  1:3561854/1‑93 (MQ=255)
                                                                        ccTGACCAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGc                <  1:4599355/93‑1 (MQ=255)
                                                                             ccAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGa           <  1:2372073/93‑1 (MQ=255)
                                                                             ccAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGa           >  2:934135/1‑93 (MQ=255)
                                                                              cAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGAc          <  1:3175899/93‑1 (MQ=255)
                                                                              cAACGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGAc          >  1:1449197/1‑93 (MQ=255)
                                                                               aaCGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGAcc         >  2:3056735/1‑93 (MQ=255)
                                                                               aaCGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGAcc         <  2:3943660/93‑1 (MQ=255)
                                                                               aaCGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGAcc         >  1:1350748/1‑93 (MQ=255)
                                                                                aCGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGACCa        >  1:4103559/1‑93 (MQ=255)
                                                                                 cGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGACCAt       <  1:3378844/93‑1 (MQ=255)
                                                                                 cGGTGCGCCGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGACCAt       <  1:1613172/93‑1 (MQ=255)
                                                                                      cgcCGGAAGCGTCCACCGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGACCATCaaaa  <  2:4395512/93‑1 (MQ=255)
                                                                                         |                                                                                         
CAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGACCATCAAAA  >  NC_000913/1132786‑1132964

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