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% 895.6 / NA 280T275S (ACG→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (147/133);  total (147/133)
*NC_0009133,221,0980GC99.6% 887.8 / ‑5.0 281T275S (ACG→AGCb3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (147/133);  minor base T (0/1);  total (147/134)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-01

Only 100 of 292 total aligned reads displayed.

GCTGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCACGTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAA  >  NC_000913/3221011‑3221184
                                                                                      ||                                                                                      
gCTGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTgt                                                                                   >  1:4725069/1‑92 (MQ=255)
gCTGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGc                                                                                   >  1:7105885/1‑93 (MQ=255)
gCTGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGc                                                                                   <  2:6945924/93‑1 (MQ=255)
   ggAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAgctttgcttt                                                                                <  1:7618371/93‑1 (MQ=255)
      atatACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTAc                                                                             <  2:5971903/93‑1 (MQ=255)
      atatACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTAc                                                                             >  1:3875562/1‑93 (MQ=255)
       tataCCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACt                                                                            >  1:3049959/1‑93 (MQ=255)
       gataCCCGGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACt                                                                            <  2:4444415/92‑1 (MQ=255)
          aCCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTc                                                                         <  1:603742/93‑1 (MQ=255)
            ccTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGa                                                                       >  1:414925/1‑93 (MQ=255)
              tGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAAc                                                                     <  2:3233793/93‑1 (MQ=255)
              tGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAAc                                                                     >  2:5141327/1‑93 (MQ=255)
              tGTTTGATGGCGAACGCGGGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAAc                                                                     <  2:5937959/93‑1 (MQ=255)
               gTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGTTTTGCTTTTACTGTCGAACa                                                                    <  1:3219837/93‑1 (MQ=255)
                tttGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACac                                                                   >  1:890722/1‑92 (MQ=255)
                tttGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAg                                                                   >  1:5255215/1‑93 (MQ=255)
                tttGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAg                                                                   >  2:4781563/1‑93 (MQ=255)
                 ttGATGGCGAACGCGTGCTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGc                                                                  <  1:4807233/93‑1 (MQ=255)
                  tGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGcc                                                                 <  1:7382775/93‑1 (MQ=255)
                    aTGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCgc                                                               <  1:2658101/93‑1 (MQ=255)
                     tGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCgca                                                              >  2:4479019/1‑93 (MQ=255)
                          aaCGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAAt                                                         <  1:6114638/93‑1 (MQ=255)
                           aCGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATg                                                        <  1:5820102/93‑1 (MQ=255)
                           aCGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATg                                                        >  2:7323732/1‑93 (MQ=255)
                            cgcgTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATgg                                                       >  2:3154438/1‑93 (MQ=255)
                            cgcgTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATgg                                                       >  1:1545429/1‑93 (MQ=255)
                             gcgTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGt                                                      <  1:4541126/93‑1 (MQ=255)
                             gcgTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACAAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGt                                                      <  1:6215503/93‑1 (MQ=255)
                               gtggtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTca                                                    >  2:5342442/1‑93 (MQ=255)
                               gtggtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTca                                                    >  2:1387702/1‑93 (MQ=255)
                                  gtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcagca                                                 <  2:4303948/93‑1 (MQ=255)
                                  gtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcagca                                                 <  2:1086848/93‑1 (MQ=255)
                                  gtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcagaa                                                 >  1:2834058/1‑91 (MQ=255)
                                    gCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGa                                               <  2:4164662/93‑1 (MQ=255)
                                    gCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGa                                               >  2:4561955/1‑93 (MQ=255)
                                     cacaGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGaa                                              <  2:2992511/93‑1 (MQ=255)
                                      acaGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGaac                                             >  1:1138403/1‑92 (MQ=255)
                                      acaGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGaac                                             >  1:2047169/1‑92 (MQ=255)
                                       cagcagCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCNGCAGCAATGGTCAGCAGAATc                                            >  1:6349636/1‑93 (MQ=255)
                                       cagcagCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATc                                            <  2:5652621/93‑1 (MQ=255)
                                         gcagcGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATccc                                          >  1:49632/1‑93 (MQ=255)
                                          gagcGCCATTGATCATATGACAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGAAATGGTCAGCAGAATcccc                                         <  1:6999584/92‑1 (MQ=255)
                                          cagcGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATcccc                                         <  2:2558806/93‑1 (MQ=255)
                                          cagcGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATcccc                                         <  2:3356771/93‑1 (MQ=255)
                                           agcGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCt                                        <  1:3920244/93‑1 (MQ=255)
                                            gcgcCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCtt                                       <  2:4879196/93‑1 (MQ=255)
                                             cgcCATTGATCATTTGGCAATTGAAAAACTGACCAGCGTCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTa                                      <  2:3889015/93‑1 (MQ=255)
                                             cgcCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTa                                      <  2:7790288/93‑1 (MQ=255)
                                             cgcCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTa                                      >  2:2043564/1‑93 (MQ=255)
                                              gcCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTAt                                     >  2:3655619/1‑93 (MQ=255)
                                              gcCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTAt                                     <  2:7660268/93‑1 (MQ=255)
                                               ccATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATc                                    <  1:1077442/93‑1 (MQ=255)
                                               ccATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATc                                    <  1:1324835/93‑1 (MQ=255)
                                               ccATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATc                                    >  1:1835849/1‑93 (MQ=255)
                                                cATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCt                                   >  2:3546034/1‑93 (MQ=255)
                                                cATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCt                                   >  1:3424661/1‑93 (MQ=255)
                                                cATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCt                                   >  1:6722069/1‑93 (MQ=255)
                                                cATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCt                                   >  2:2722111/1‑93 (MQ=255)
                                                 aTTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCtt                                  >  1:6611743/1‑93 (MQ=255)
                                                 aTTGATCATTTGGCAATTGAAAAACTGACCAGAGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCtt                                  >  2:6585705/1‑93 (MQ=255)
                                                  ttGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCttt                                 >  2:517923/1‑93 (MQ=255)
                                                    gATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTAc                               >  1:5993919/1‑93 (MQ=255)
                                                     aTCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTAcc                              <  2:5607497/93‑1 (MQ=255)
                                                     aTCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTAcc                              <  1:2437972/93‑1 (MQ=255)
                                                     aTCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTAcc                              >  1:7194963/1‑93 (MQ=255)
                                                     aTCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTAcc                              <  1:570843/93‑1 (MQ=255)
                                                      tCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCa                             >  2:2363292/1‑93 (MQ=255)
                                                      tCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCa                             >  2:6412647/1‑93 (MQ=255)
                                                       cATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCAt                            <  1:7459899/93‑1 (MQ=255)
                                                       cATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCAt                            <  1:518441/93‑1 (MQ=255)
                                                         tttGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCt                          >  1:1417107/1‑93 (MQ=255)
                                                         tttGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCt                          >  1:856237/1‑93 (MQ=255)
                                                         tttGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCt                          >  1:2599690/1‑93 (MQ=255)
                                                             gCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTc                      <  2:7600582/93‑1 (MQ=255)
                                                              cAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCa                     <  1:2258228/93‑1 (MQ=255)
                                                              cAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCa                     >  1:5321137/1‑93 (MQ=255)
                                                               aaTTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCAt                    >  1:791626/1‑93 (MQ=255)
                                                               aaTTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCAt                    >  1:4325701/1‑93 (MQ=255)
                                                               aaTTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCAt                    >  2:1968934/1‑93 (MQ=255)
                                                                 ttGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGa                  >  2:871788/1‑93 (MQ=255)
                                                                 ttGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGa                  <  1:1522297/93‑1 (MQ=255)
                                                                 ttGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGa                  <  2:6835472/93‑1 (MQ=255)
                                                                 ttGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGa                  <  2:2465311/93‑1 (MQ=255)
                                                                   gAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACg                >  1:1993186/1‑93 (MQ=255)
                                                                    aaaaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGc               <  1:5528035/93‑1 (MQ=255)
                                                                    aaaaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGc               <  2:2819050/93‑1 (MQ=255)
                                                                     aaaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCt              <  1:6215545/93‑1 (MQ=255)
                                                                     aaaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCt              <  2:7485361/93‑1 (MQ=255)
                                                                       aaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGa            >  2:1825054/1‑93 (MQ=255)
                                                                       aaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGa            >  1:1224648/1‑93 (MQ=255)
                                                                         cTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGaaa          <  2:7713121/93‑1 (MQ=255)
                                                                         cTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGaaa          <  1:3143592/93‑1 (MQ=255)
                                                                           gACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGa        <  2:155240/93‑1 (MQ=255)
                                                                           gACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGa        >  2:101651/1‑93 (MQ=255)
                                                                           gACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGa        <  1:5084017/93‑1 (MQ=255)
                                                                            aCCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGAc       <  1:5451138/93‑1 (MQ=255)
                                                                            aCCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGAc       <  2:1249769/93‑1 (MQ=255)
                                                                             ccAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACg      >  1:1775559/1‑93 (MQ=255)
                                                                                gcgcCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCa   <  2:3770031/93‑1 (MQ=255)
                                                                                 cgcCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCaa  <  2:3087493/93‑1 (MQ=255)
                                                                                      ||                                                                                      
GCTGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCACGTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAA  >  NC_000913/3221011‑3221184

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