Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,222,944 4 bp→AGCA coding (2671‑2674/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,222,9440CA99.8% 1828.3 / NA 472Q891S (CAA→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (229/242);  total (229/243)
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,222,9450AG99.8% 1788.0 / NA 472Q891S (CAA→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (229/241);  total (229/243)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
*NC_0009133,222,9460AC99.8% 1789.9 / NA 472Q891S (CAA→AGCb3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (229/242);  minor base T (0/1);  total (229/243)
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,222,9470GA98.3% 1712.5 / NA 471A892T (GCC→ACC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (228/234);  minor base C (0/8);  total (228/243)
Fisher's exact test for biased strand distribution p-value = 7.54e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 483 total aligned reads displayed.

TAACGGCGAATACGATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCCAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACTAAC  >  NC_000913/3222850‑3223030
                                                                                              ||||                                                                                   
          tACGATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGTCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGAt                                                                             >  1:12104540/1‑96 (MQ=255)
          tACGATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGAt                                                                             <  1:2162814/96‑1 (MQ=255)
            cGATCAGGTGGCGTATTACGGTCGTGGACCGGGGGAAGACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGAGCAGCACCGTCGATGc                                                                           >  1:12822976/1‑96 (MQ=255)
            cGATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGc                                                                           <  1:11582504/96‑1 (MQ=255)
             gATCAGGTGGCGTATTACGGTCGTGGACCGGGCGATAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGcc                                                                          <  2:11768327/96‑1 (MQ=255)
             gATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGcc                                                                          <  2:7215058/96‑1 (MQ=255)
             gATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGcc                                                                          <  2:6066624/96‑1 (MQ=255)
             gATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGcc                                                                          >  1:2437393/1‑96 (MQ=255)
              aTCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCa                                                                         >  1:6757572/1‑96 (MQ=255)
                cAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATg                                                                       <  1:4236448/96‑1 (MQ=255)
                 aGGTGGCGTATTACGGTCGTGGACCGGTCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGt                                                                      <  2:12531865/96‑1 (MQ=255)
                 aGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGt                                                                      >  2:7686947/1‑96 (MQ=255)
                 aGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGt                                                                      <  2:2338450/96‑1 (MQ=255)
                  ggTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGtt                                                                     <  1:7075990/96‑1 (MQ=255)
                   gTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTc                                                                    >  1:594421/1‑96 (MQ=255)
                    tGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCg                                                                   <  2:6625383/96‑1 (MQ=255)
                    tGGCGTATTACGGTCGGGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCg                                                                   >  1:1891040/1‑96 (MQ=255)
                        gTATTACGGTCGTGGACCGGGCGCAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGaa                                                               <  1:10137713/96‑1 (MQ=255)
                        gTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGaa                                                               >  1:453790/1‑96 (MQ=255)
                         tATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAAc                                                              <  1:2715612/96‑1 (MQ=255)
                          aTTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAGCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACt                                                             >  2:11933863/1‑96 (MQ=255)
                          aTTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCCCCGTCGATGCCATGTTCGAGAACt                                                             <  2:315107/96‑1 (MQ=255)
                          aTTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACt                                                             >  2:8197953/1‑96 (MQ=255)
                           ttACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTa                                                            >  1:5031240/1‑96 (MQ=255)
                            tACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTAt                                                           <  2:8010915/96‑1 (MQ=255)
                             aCGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATc                                                          >  2:6469431/1‑96 (MQ=255)
                              cGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATcc                                                         <  1:403588/96‑1 (MQ=255)
                               ggTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATccc                                                        >  2:2835481/1‑96 (MQ=255)
                                gTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCATCGTCGATGCCATGTTCGAGAACTATctct                                                       >  1:9149168/1‑96 (MQ=255)
                                gTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCt                                                       >  1:1975953/1‑96 (MQ=255)
                                gTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCt                                                       <  1:6458242/96‑1 (MQ=255)
                                ctCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCCTGTTCGAGAACTATCCCt                                                       <  1:3150750/95‑1 (MQ=255)
                                 tCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCtt                                                      <  2:9850360/96‑1 (MQ=255)
                                 tCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCtt                                                      <  1:5593862/96‑1 (MQ=255)
                                  cGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTc                                                     >  1:13523099/1‑96 (MQ=255)
                                  cGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTc                                                     <  1:8747704/96‑1 (MQ=255)
                                   gTGGACCGGGCGAAAACTCCGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTcc                                                    <  2:2899709/96‑1 (MQ=255)
                                    tGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTccc                                                   >  1:7310062/1‑96 (MQ=255)
                                     ggACCGGGCGATAACTACGCCGCCAGCCAGCAGGCTATCCTCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCg                                                  <  1:14331311/96‑1 (MQ=255)
                                     ggACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCg                                                  <  1:5163260/96‑1 (MQ=255)
                                     ggACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCg                                                  >  1:9004201/1‑96 (MQ=255)
                                      gACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGc                                                 >  2:1124693/1‑96 (MQ=255)
                                      gACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGc                                                 <  2:872115/96‑1 (MQ=255)
                                       aCCGGCCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCa                                                <  2:4064097/96‑1 (MQ=255)
                                         cGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGa                                              >  2:10278678/1‑96 (MQ=255)
                                          gggCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaa                                             <  1:3611480/96‑1 (MQ=255)
                                          gggCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaa                                             >  2:5255145/1‑96 (MQ=255)
                                           ggCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaac                                            <  2:7343001/96‑1 (MQ=255)
                                            gCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaaca                                           <  2:5670054/96‑1 (MQ=255)
                                             cGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaa                                          >  1:10999077/1‑96 (MQ=255)
                                              gAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaac                                         >  1:7961899/1‑96 (MQ=255)
                                               aaaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCCCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACg                                        <  2:5015697/96‑1 (MQ=255)
                                               aaaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACg                                        <  1:4936343/96‑1 (MQ=255)
                                                 aaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGt                                      <  1:7357058/96‑1 (MQ=255)
                                                  aCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTa                                     >  2:12310178/1‑96 (MQ=255)
                                                   cTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGtga                                    <  1:11121198/96‑3 (MQ=255)
                                                   cTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTaa                                    <  1:6905014/96‑1 (MQ=255)
                                                    tACGGCGACAGCCAGCAGGCTAACATCAACGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAc                                   <  2:3655834/96‑1 (MQ=255)
                                                    tACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAc                                   <  1:8294434/96‑1 (MQ=255)
                                                     aCGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGTACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAcc                                  <  1:9340065/96‑1 (MQ=255)
                                                      cGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCg                                 >  2:13801870/1‑96 (MQ=255)
                                                       gCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGt                                >  2:409484/1‑96 (MQ=255)
                                                         cGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCa                              >  1:9691547/1‑96 (MQ=255)
                                                          gACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAg                             >  1:13219458/1‑96 (MQ=255)
                                                            cagccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCa                           >  2:11455158/1‑96 (MQ=255)
                                                            cagccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCa                           >  2:13300930/1‑96 (MQ=255)
                                                             agccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTACCCGTCAGcct                          <  1:10828334/96‑3 (MQ=255)
                                                             agccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCTTCAGCAt                          >  1:634643/1‑96 (MQ=255)
                                                             agccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCAt                          >  1:7736274/1‑96 (MQ=255)
                                                              gccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATg                         <  2:7451040/96‑1 (MQ=255)
                                                              gccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATg                         <  1:10006881/96‑1 (MQ=255)
                                                              gccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATg                         >  2:8579441/1‑96 (MQ=255)
                                                               ccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGt                        >  2:10511932/1‑96 (MQ=255)
                                                               ccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGt                        <  1:7826186/96‑1 (MQ=255)
                                                                cagcagGCTAACATCATCGATATCTGGCGCAGCCCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTc                       <  2:9471920/96‑1 (MQ=255)
                                                                cagcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAGCGGTAACCGTCAGCATGTc                       >  2:6801896/1‑96 (MQ=255)
                                                                cagcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTc                       >  1:1592339/1‑96 (MQ=255)
                                                                cagcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTc                       >  2:4787411/1‑96 (MQ=255)
                                                                 agcagGCTAACATCATCGATATCTGGCGCATCACCGTCGATGCCATGTTCGAGAACTATCCCATCCCGCAGAACAACGGTAACCGTCAGCACGTcc                      <  2:14087766/96‑1 (MQ=255)
                                                                 agcagGCTAACATCATCGATATCTGGCGCAGCACCGTTGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGttc                      <  1:10652932/96‑3 (MQ=255)
                                                                 agcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTcc                      <  2:12058981/96‑1 (MQ=255)
                                                                 agcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTcc                      >  2:9595163/1‑96 (MQ=255)
                                                                  gcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGGTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCg                     >  2:6349641/1‑96 (MQ=255)
                                                                   cagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGc                    >  2:10894585/1‑96 (MQ=255)
                                                                   cagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGc                    >  2:12989173/1‑96 (MQ=255)
                                                                   cagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGc                    <  2:437132/96‑1 (MQ=255)
                                                                    agGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCt                   >  1:10396360/1‑96 (MQ=255)
                                                                     ggCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTg                  <  1:11997151/96‑1 (MQ=255)
                                                                      gCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTgg                 <  2:2720184/96‑1 (MQ=255)
                                                                       cTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGa                <  1:9964191/96‑1 (MQ=255)
                                                                         aaCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACg              <  1:8617178/96‑1 (MQ=255)
                                                                           catcatCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGc            >  2:3451776/1‑96 (MQ=255)
                                                                            atcatcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCa           <  1:3464458/96‑1 (MQ=255)
                                                                            atcatcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCa           <  2:8981097/96‑1 (MQ=255)
                                                                               atcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTg        >  1:12467537/1‑96 (MQ=255)
                                                                               atcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCACGTCCGCTGGACGGCACTg        >  1:13974153/1‑96 (MQ=255)
                                                                                  gATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACt     >  2:13730858/1‑96 (MQ=255)
                                                                                  gATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACt     <  2:3967757/96‑1 (MQ=255)
                                                                                  gATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACt     >  1:4747810/1‑96 (MQ=255)
                                                                                   atatCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACTa    >  2:1533644/1‑96 (MQ=255)
                                                                                              ||||                                                                                   
TAACGGCGAATACGATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCCAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACTAAC  >  NC_000913/3222850‑3223030

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