Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,221,097 2 bp→GC coding (824‑825/3093 nt) b3076 → evolved beta‑D‑galactosidase alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,221,0970CG99.4% 524.4 / NA 162T275S (ACG→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (88/73);  minor base T (1/0);  total (89/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,221,0980GC100.0% 527.2 / NA 162T275S (ACG→AGCb3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (89/73);  total (89/73)

Only 100 of 168 total aligned reads displayed.

CGCTGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCACGTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAA  >  NC_000913/3221010‑3221184
                                                                                       ||                                                                                      
cgCTGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTg                                                                                    >  2:3875340/1‑93 (MQ=255)
   tGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCtt                                                                                 >  1:4176483/1‑93 (MQ=255)
   tGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCtt                                                                                 >  1:1606707/1‑93 (MQ=255)
   tGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCtt                                                                                 >  2:4517545/1‑93 (MQ=255)
    ggAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAgctttgcttt                                                                                >  1:4506510/1‑93 (MQ=255)
    ggAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGtttt                                                                                >  1:263685/1‑93 (MQ=255)
        tataCCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACt                                                                            <  2:4165350/93‑1 (MQ=255)
        tataCCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACt                                                                            >  1:3260006/1‑93 (MQ=255)
         ataCCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTg                                                                           <  1:4247526/93‑1 (MQ=255)
         ataCCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTg                                                                           <  1:2757987/93‑1 (MQ=255)
          taCCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGt                                                                          >  1:2077878/1‑93 (MQ=255)
          taCCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGt                                                                          >  1:4561547/1‑93 (MQ=255)
          taCCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGt                                                                          >  1:1789890/1‑93 (MQ=255)
          taCCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGt                                                                          <  2:3313724/93‑1 (MQ=255)
           aCCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTc                                                                         >  2:1591423/1‑93 (MQ=255)
            cccTGTTTGATGGCGAACGCGTGGGGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCg                                                                        >  2:1847545/1‑93 (MQ=255)
              cTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGaa                                                                      <  2:2056613/93‑1 (MQ=255)
               tGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAAc                                                                     >  2:547617/1‑93 (MQ=255)
                gTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACa                                                                    >  2:1484010/1‑93 (MQ=255)
                gTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACa                                                                    <  1:3288654/93‑1 (MQ=255)
                 tttGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAg                                                                   <  2:4623964/93‑1 (MQ=255)
                   tGATGGCGAACGCGTGGTGCACAGCATCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGcc                                                                 <  1:2489375/93‑1 (MQ=255)
                   tGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGcc                                                                 <  1:3822878/93‑1 (MQ=255)
                    gATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCg                                                                <  2:369223/93‑1 (MQ=255)
                     aTGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAAATGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCgc                                                               >  2:2318251/1‑93 (MQ=255)
                      tGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCgca                                                              >  2:4711483/1‑93 (MQ=255)
                      tGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCgca                                                              <  1:787429/93‑1 (MQ=255)
                        gCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCgcagc                                                            >  2:743349/1‑93 (MQ=255)
                          gAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCaa                                                          <  2:178757/93‑1 (MQ=255)
                           aaCGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAAt                                                         >  2:2174917/1‑93 (MQ=255)
                           aaCGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAAt                                                         <  2:3927391/93‑1 (MQ=255)
                           aaCGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAAt                                                         >  2:4298414/1‑93 (MQ=255)
                           aaCGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAAt                                                         <  1:4128514/93‑1 (MQ=255)
                           aaCGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAAt                                                         <  1:4074131/93‑1 (MQ=255)
                            aCGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATg                                                        <  1:2002983/93‑1 (MQ=255)
                            aCGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCACACAGCCGCAGCCATg                                                        >  1:4078967/1‑93 (MQ=255)
                              gcgTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGt                                                      <  2:3065722/93‑1 (MQ=255)
                                gtggtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTca                                                    <  2:301703/93‑1 (MQ=255)
                                 tggtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcag                                                   >  1:3440007/1‑93 (MQ=255)
                                 tggtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcag                                                   <  1:3159721/93‑1 (MQ=255)
                                  ggtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcagc                                                  <  1:2632411/93‑1 (MQ=255)
                                  ggtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcagc                                                  >  2:2777933/1‑93 (MQ=255)
                                   gtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcagca                                                 <  2:4486495/93‑1 (MQ=255)
                                    tgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcagcag                                                >  1:3389202/1‑93 (MQ=255)
                                     gCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGa                                               <  1:3929602/93‑1 (MQ=255)
                                        cagcagCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATc                                            <  1:1619909/93‑1 (MQ=255)
                                        cagcagCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATc                                            >  2:1123958/1‑93 (MQ=255)
                                        cagcagCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATc                                            <  2:775111/93‑1 (MQ=255)
                                         agcagcGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATcc                                           >  2:1675758/1‑93 (MQ=255)
                                          gcagcGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATccc                                          >  1:2766105/1‑93 (MQ=255)
                                          gcagcGCCATTGATCATTTGGCAATTGAAAAACTGACCAGAGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATCGTCAGCAGAATccc                                          >  1:3297118/1‑93 (MQ=255)
                                           cagcGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATcccc                                         <  1:797702/93‑1 (MQ=255)
                                               gcCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTAt                                     >  1:1102516/1‑93 (MQ=255)
                                               gcCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTAt                                     >  2:3317447/1‑93 (MQ=255)
                                               gcCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTAt                                     <  2:704477/93‑1 (MQ=255)
                                                ccATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATc                                    >  2:2864821/1‑93 (MQ=255)
                                                ccATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATc                                    >  2:1412621/1‑93 (MQ=255)
                                                ccATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATc                                    <  2:875762/93‑1 (MQ=255)
                                                  aTTGATCATTTGGCAATTGAAAACCTGACCAGCGCCAGCTTTGCTTTTACTTTCGAACAGTCGCAGGAATTGGCAGCACAATCCCCTTATcat                                  >  2:1673835/1‑91 (MQ=255)
                                                  aTTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCtt                                  >  2:4449340/1‑93 (MQ=255)
                                                  aTTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCtt                                  <  1:2415418/93‑1 (MQ=255)
                                                   ttGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCttt                                 <  2:3960159/93‑1 (MQ=255)
                                                    tGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTa                                >  1:3106085/1‑93 (MQ=255)
                                                    tGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTa                                >  1:4490564/1‑93 (MQ=255)
                                                      aTCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTAcc                              <  1:2340702/93‑1 (MQ=255)
                                                      aTCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTAcc                              >  1:2613650/1‑93 (MQ=255)
                                                       tCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATTTTTAcct                             >  1:4548374/1‑92 (MQ=255)
                                                        cATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCAt                            >  1:1341858/1‑93 (MQ=255)
                                                         aTTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATc                           >  1:4598052/1‑93 (MQ=255)
                                                         aTTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATc                           >  2:201954/1‑93 (MQ=255)
                                                          tttGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCt                          >  1:3045435/1‑93 (MQ=255)
                                                          tttGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCt                          >  2:1225248/1‑93 (MQ=255)
                                                           ttGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTg                         >  1:2394402/1‑93 (MQ=255)
                                                           ttGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTg                         >  1:2300696/1‑93 (MQ=255)
                                                            tGGCAATTATAAAACTGACCAGCGCCAGCTTTGCTTTTAATGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATATTTACCATCTgg                        <  2:4219975/93‑1 (MQ=255)
                                                             ggCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGt                       <  1:144359/93‑1 (MQ=255)
                                                              gCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTc                      >  2:3220842/1‑93 (MQ=255)
                                                                aaTTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCAt                    <  2:248847/93‑1 (MQ=255)
                                                                 aTTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATg                   >  2:1010594/1‑93 (MQ=255)
                                                                  ttGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGa                  >  1:3038763/1‑93 (MQ=255)
                                                                  ttGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGa                  >  1:3908995/1‑93 (MQ=255)
                                                                   tGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGAc                 <  2:3950675/93‑1 (MQ=255)
                                                                    gAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACg                >  1:3618479/1‑93 (MQ=255)
                                                                     aaaaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGc               >  2:3443293/1‑93 (MQ=255)
                                                                     aaaaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGc               <  2:4392536/93‑1 (MQ=255)
                                                                     aaaaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGc               <  2:3782327/93‑1 (MQ=255)
                                                                     aaaaaCTAACCAGAGCCAGCTTTGATTTTCATGACGATCAGCCGCAAAAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGc               <  2:1220428/93‑1 (MQ=255)
                                                                      aaaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCt              >  2:3145662/1‑93 (MQ=255)
                                                                      aaaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCt              >  2:1975011/1‑93 (MQ=255)
                                                                       aaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTg             >  2:243398/1‑93 (MQ=255)
                                                                       aaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTg             <  1:1129430/93‑1 (MQ=255)
                                                                            gACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGa        <  2:2719014/93‑1 (MQ=255)
                                                                            gACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGa        >  1:290070/1‑93 (MQ=255)
                                                                            gACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGa        <  1:1253791/93‑1 (MQ=255)
                                                                               cAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGc     <  2:252331/93‑1 (MQ=255)
                                                                               cAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGc     <  2:1267024/93‑1 (MQ=255)
                                                                                aGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGcc    <  1:4435808/93‑1 (MQ=255)
                                                                                  cgcCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCaa  <  2:434915/93‑1 (MQ=255)
                                                                                  cgcCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCaa  <  1:598613/93‑1 (MQ=255)
                                                                                  cgcCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCaa  <  2:572492/93‑1 (MQ=255)
                                                                                       ||                                                                                      
CGCTGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCACGTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAA  >  NC_000913/3221010‑3221184

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