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,0970CG100.0% 564.4 / NA 177T275S (ACG→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (86/91);  total (86/91)
*NC_0009133,221,0980GC100.0% 565.2 / NA 177T275S (ACG→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (86/91);  total (86/91)

Only 100 of 180 total aligned reads displayed.

GCTGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCACGTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAA  >  NC_000913/3221011‑3221184
                                                                                      ||                                                                                      
gCTGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGc                                                                                   >  2:4635951/1‑93 (MQ=255)
  tGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCtt                                                                                 >  2:2197658/1‑93 (MQ=255)
  tGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCtt                                                                                 >  2:3593054/1‑93 (MQ=255)
  tGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCtt                                                                                 >  2:3518262/1‑93 (MQ=255)
   ggAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAgctttgcttt                                                                                >  2:3154848/1‑93 (MQ=255)
    gAATATACCCTGTTTGATGGCGAACGCGTGGTGCATAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCtttt                                                                               <  2:2335669/93‑1 (MQ=255)
    gAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCtttt                                                                               <  1:3361193/93‑1 (MQ=255)
     aaTATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTa                                                                              <  1:1259033/93‑1 (MQ=255)
      atatACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTAc                                                                             <  2:1713545/93‑1 (MQ=255)
       tataCCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACt                                                                            <  1:2829854/93‑1 (MQ=255)
       tataCCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACt                                                                            <  1:2905095/93‑1 (MQ=255)
        ataCCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTg                                                                           <  1:269678/93‑1 (MQ=255)
        ataCCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTg                                                                           <  1:2008205/93‑1 (MQ=255)
         taCCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGt                                                                          >  2:3536440/1‑93 (MQ=255)
         taCCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGt                                                                          <  1:672049/93‑1 (MQ=255)
         taCCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGt                                                                          >  1:3483072/1‑93 (MQ=255)
             cTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGaa                                                                      <  1:3573886/93‑1 (MQ=255)
             cTGTTTGATGGCGAACGCGTGGGGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGGGCCAGCGTTGCTTTTACTGTCGaa                                                                      >  2:2912951/1‑93 (MQ=255)
              cgTTTGATGGCGAACGCGTGGTGCACAGCAGCTCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAAc                                                                     <  1:4455835/92‑1 (MQ=255)
               tttttATGGCGAATGCCGGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACa                                                                    <  2:2307585/92‑1 (MQ=255)
               gTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACa                                                                    <  2:3848697/93‑1 (MQ=255)
                 ttGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGc                                                                  >  1:2423365/1‑93 (MQ=255)
                 ttGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGc                                                                  <  2:342629/93‑1 (MQ=255)
                  tGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGcc                                                                 >  2:581655/1‑93 (MQ=255)
                   gATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCg                                                                >  2:2886818/1‑93 (MQ=255)
                    aTGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCgc                                                               >  1:778171/1‑93 (MQ=255)
                      ggCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCgcag                                                             >  2:3985677/1‑93 (MQ=255)
                      ggCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCgcag                                                             >  2:959658/1‑93 (MQ=255)
                       gCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCgcagc                                                            <  1:1120546/93‑1 (MQ=255)
                         gAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCaa                                                          <  2:4184186/93‑1 (MQ=255)
                         gAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCaa                                                          <  2:1285041/93‑1 (MQ=255)
                          aaCGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAAt                                                         >  2:2797141/1‑93 (MQ=255)
                          aaCGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAAt                                                         <  2:1337575/93‑1 (MQ=255)
                          aaCGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAAt                                                         >  1:1056450/1‑93 (MQ=255)
                          aaCGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAAt                                                         <  1:1523980/93‑1 (MQ=255)
                           aCGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATg                                                        <  2:424089/93‑1 (MQ=255)
                             gcgTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGt                                                      <  2:163355/93‑1 (MQ=255)
                              cgTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTc                                                     >  2:1449549/1‑93 (MQ=255)
                               gtggtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTca                                                    <  2:3351675/93‑1 (MQ=255)
                               gtggtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTca                                                    >  1:4664002/1‑93 (MQ=255)
                               gtggtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTca                                                    >  1:1029994/1‑93 (MQ=255)
                                tggtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcag                                                   >  2:4465953/1‑93 (MQ=255)
                                 ggtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcagc                                                  <  2:4145289/93‑1 (MQ=255)
                                  gtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcagca                                                 <  1:1704154/93‑1 (MQ=255)
                                   tgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcagcag                                                >  1:3824898/1‑93 (MQ=255)
                                     cacaGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGaa                                              >  1:1187878/1‑93 (MQ=255)
                                     cacaGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGaa                                              <  2:2582822/93‑1 (MQ=255)
                                      acaGCAGCGCTATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAAt                                             <  2:4306714/93‑1 (MQ=255)
                                      acaGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAAt                                             >  1:1364898/1‑93 (MQ=255)
                                         gcagcGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATccc                                          >  1:1844148/1‑93 (MQ=255)
                                         gcagcGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATccc                                          <  1:2809190/93‑1 (MQ=255)
                                          cagcGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATcccc                                         <  1:3329285/93‑1 (MQ=255)
                                           agcGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCt                                        <  1:800557/93‑1 (MQ=255)
                                            gcgcCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCtt                                       <  2:2858699/93‑1 (MQ=255)
                                            gcgcCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCtt                                       <  1:4242886/93‑1 (MQ=255)
                                            gcgcCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCtt                                       <  1:4155659/93‑1 (MQ=255)
                                               ccATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATc                                    >  1:3601465/1‑93 (MQ=255)
                                               ccATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATc                                    <  2:1861474/93‑1 (MQ=255)
                                                cATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCt                                   >  1:1374969/1‑93 (MQ=255)
                                                 aTTGATCATTTGGCAATTGAAAAACTGACCAGGGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCtt                                  >  1:1605720/1‑93 (MQ=255)
                                                 aTTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCtt                                  >  2:3715891/1‑93 (MQ=255)
                                                 aTTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCtt                                  <  1:30276/93‑1 (MQ=255)
                                                 aTTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCtt                                  >  2:762855/1‑93 (MQ=255)
                                                 aTTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCtt                                  >  1:3644717/1‑93 (MQ=255)
                                                  ttGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCttt                                 >  1:1782208/1‑93 (MQ=255)
                                                   tGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCtttc                                >  1:3779672/1‑92 (MQ=255)
                                                   tGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTa                                >  2:46348/1‑93 (MQ=255)
                                                    gATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTAc                               <  1:3939307/93‑1 (MQ=255)
                                                     aTCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTAcc                              >  1:3655890/1‑93 (MQ=255)
                                                      tCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCa                             <  1:3261166/93‑1 (MQ=255)
                                                        aTTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATc                           >  2:1311542/1‑93 (MQ=255)
                                                        aTTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATc                           >  1:4500649/1‑93 (MQ=255)
                                                        aTTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACACCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATc                           >  2:3316659/1‑93 (MQ=255)
                                                         tttGGCAATTTAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCt                          <  1:4014251/93‑1 (MQ=255)
                                                         tttGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCt                          <  1:4299105/93‑1 (MQ=255)
                                                         tttGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCt                          >  1:704584/1‑93 (MQ=255)
                                                             gCAATTGAAAAACTGACCAGCTCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTc                      <  2:76627/93‑1 (MQ=255)
                                                              cAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCa                     <  2:1782867/93‑1 (MQ=255)
                                                               aaTTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCAt                    >  1:2108345/1‑93 (MQ=255)
                                                                aTTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATg                   >  1:3811934/1‑93 (MQ=255)
                                                                 ttGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGa                  >  1:440119/1‑93 (MQ=255)
                                                                 ttGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGa                  >  1:316882/1‑93 (MQ=255)
                                                                 ttGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGa                  <  1:860216/93‑1 (MQ=255)
                                                                  tGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGAc                 <  2:2407111/93‑1 (MQ=255)
                                                                    aaaaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGc               <  2:520402/93‑1 (MQ=255)
                                                                    aaaaaCTGACCAGCGCCAGCTTTGCTTTGACTGGCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGc               <  2:831691/93‑1 (MQ=255)
                                                                     aaaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCt              <  2:1568886/93‑1 (MQ=255)
                                                                      aaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTg             >  2:3222086/1‑93 (MQ=255)
                                                                       aaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGa            >  2:1499582/1‑93 (MQ=255)
                                                                       aaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGa            <  2:3427910/93‑1 (MQ=255)
                                                                       aaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGa            <  1:1938766/93‑1 (MQ=255)
                                                                        aCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGaa           >  2:3941076/1‑93 (MQ=255)
                                                                          tGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAg         >  2:366813/1‑93 (MQ=255)
                                                                          tGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAg         >  1:198128/1‑93 (MQ=255)
                                                                           gACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGa        <  2:1972089/93‑1 (MQ=255)
                                                                           gACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGa        <  1:2106262/93‑1 (MQ=255)
                                                                            aCCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGAc       <  2:872304/93‑1 (MQ=255)
                                                                              cAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGc     <  1:2981623/93‑1 (MQ=255)
                                                                              cAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGc     <  2:718816/93‑1 (MQ=255)
                                                                                 cgcCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCaa  <  2:2010782/93‑1 (MQ=255)
                                                                                      ||                                                                                      
GCTGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCACGTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAA  >  NC_000913/3221011‑3221184

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