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.8% 2051.3 / ‑5.0 528T275S (ACG→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (276/251);  minor base T (1/0);  total (277/251)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.39e-01
*NC_0009133,221,0980GC99.2% 2041.8 / ‑4.6 528T275S (ACG→AGCb3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (273/250);  minor base A (4/0);  total (277/251)
Fisher's exact test for biased strand distribution p-value = 1.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 543 total aligned reads displayed.

CGACGCTGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCACGTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGG  >  NC_000913/3221007‑3221187
                                                                                          ||                                                                                         
cgacgCTGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTg                                                                                       >  2:4981211/1‑96 (MQ=255)
  acgCTGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCt                                                                                     <  2:4310705/96‑1 (MQ=255)
     cTGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCtttt                                                                                  <  2:1649442/96‑1 (MQ=255)
     cTGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCtttt                                                                                  >  1:12841845/1‑96 (MQ=255)
      tcgAATATACCCTGTTTGATGGCGAACGCGTGGGGCACAGCAGCGCCATTGAGCATTTGGCAATTGAAAAACTAACCAGCGCCAGCTTTGCTTTTa                                                                                 >  1:10600288/3‑96 (MQ=255)
       ggAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTAc                                                                                <  2:10358198/96‑1 (MQ=255)
        gAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACt                                                                               >  1:5583863/1‑96 (MQ=255)
        gAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACt                                                                               <  1:6370138/96‑1 (MQ=255)
         aaTATACCCTGTTTGATGGCGAACGCGTGGTGCACCGCAGCCCCATTGATCATTTGGCAATTGAAAAACCGACCAGCGCCAGCTTTGCTTTTACTg                                                                              <  1:4003712/96‑1 (MQ=255)
          atatACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGt                                                                             >  2:4064086/1‑96 (MQ=255)
          atatACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGt                                                                             >  2:5451797/1‑96 (MQ=255)
           tataCCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTc                                                                            >  2:961545/1‑96 (MQ=255)
            ataCCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCg                                                                           >  1:2651087/1‑96 (MQ=255)
              aCCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGaa                                                                         >  1:8345101/1‑96 (MQ=255)
               cccTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAAc                                                                        <  1:12460000/96‑1 (MQ=255)
               cccGGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAAc                                                                        >  2:10101323/1‑96 (MQ=255)
                    tttGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCg                                                                   >  1:1563890/1‑96 (MQ=255)
                    tttGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCg                                                                   <  1:566285/96‑1 (MQ=255)
                     ttGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCgc                                                                  <  2:7689443/96‑1 (MQ=255)
                     ttGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCgc                                                                  <  1:6009722/96‑1 (MQ=255)
                      tGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCgca                                                                 <  1:7573971/96‑1 (MQ=255)
                           gCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAAt                                                            <  2:12055471/96‑1 (MQ=255)
                            cGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATg                                                           <  1:4841729/96‑1 (MQ=255)
                             gAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATgg                                                          >  1:6836959/1‑96 (MQ=255)
                             gAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATgg                                                          >  2:9783931/1‑96 (MQ=255)
                             gAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATgg                                                          >  1:11659998/1‑96 (MQ=255)
                              aaCGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGt                                                         >  2:10743551/1‑96 (MQ=255)
                              aaCGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGt                                                         >  1:6216419/1‑96 (MQ=255)
                               aCGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTc                                                        <  2:2843419/96‑1 (MQ=255)
                                cgcgTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTca                                                       <  1:9836944/96‑1 (MQ=255)
                                cgcgTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTca                                                       <  1:2783515/96‑1 (MQ=255)
                                 gcgTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcag                                                      <  1:1303659/96‑1 (MQ=255)
                                 gcgTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcag                                                      >  2:8493715/1‑96 (MQ=255)
                                  cgTGGTGCACAGCAGCGCCTTTGATCATTTGGCAATTGAAAAACTGACTAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcagc                                                     >  2:12773246/1‑96 (MQ=255)
                                  cgTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcagc                                                     >  1:8984805/1‑96 (MQ=255)
                                   gtggtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcagca                                                    <  1:12086368/96‑1 (MQ=255)
                                     ggtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGa                                                  >  1:2299966/1‑96 (MQ=255)
                                      gtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGaa                                                 <  1:4538807/96‑1 (MQ=255)
                                      gtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGaa                                                 <  1:2446327/96‑1 (MQ=255)
                                       tgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAAt                                                <  2:6866603/96‑1 (MQ=255)
                                       tgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAAt                                                <  2:2447424/96‑1 (MQ=255)
                                       tgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAAt                                                <  2:7253873/96‑1 (MQ=255)
                                        gCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATc                                               >  2:2612366/1‑96 (MQ=255)
                                        gCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATc                                               <  1:9519118/96‑1 (MQ=255)
                                         cacaGCAGCGCCATTGATCATTTGGCAATTCAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAAttc                                              <  1:10283590/96‑3 (MQ=255)
                                          acaGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATccc                                             >  1:6568772/1‑96 (MQ=255)
                                             gcagcGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCtt                                          <  2:10212899/96‑1 (MQ=255)
                                              cagcGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTa                                         <  2:11228378/96‑1 (MQ=255)
                                               aacgcCATTGATCATTTGGCAATGGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTAt                                        >  1:14333637/3‑96 (MQ=255)
                                                  gcCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCtt                                     >  2:6407473/1‑96 (MQ=255)
                                                  gcCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGATTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCATCAGAATCCCCTTATCtt                                     >  1:11354279/1‑96 (MQ=255)
                                                   ccATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCttt                                    <  2:5204726/96‑1 (MQ=255)
                                                   ccATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCttt                                    >  2:7070768/1‑96 (MQ=255)
                                                     aTTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTAc                                  <  1:471614/96‑1 (MQ=255)
                                                      ttGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTAcc                                 <  1:9174171/96‑1 (MQ=255)
                                                      ttGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTAcc                                 >  1:11008982/1‑96 (MQ=255)
                                                        gATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCAt                               >  1:13378760/1‑96 (MQ=255)
                                                        gATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCAt                               <  1:5266334/96‑1 (MQ=255)
                                                        gATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCAt                               >  2:14282697/1‑96 (MQ=255)
                                                         agcATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATc                              >  1:7365381/3‑96 (MQ=255)
                                                         aTCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCGGAATCCCCTTATCTTAACCATc                              >  2:11768452/1‑96 (MQ=255)
                                                         aTCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGGATCCCCTTATCTTTACCATc                              >  2:13970968/1‑96 (MQ=255)
                                                         aTCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATc                              <  2:12283308/96‑1 (MQ=255)
                                                          tCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCt                             <  2:6774235/96‑1 (MQ=255)
                                                          tCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCt                             <  2:13315226/96‑1 (MQ=255)
                                                          tCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCt                             <  1:6751819/96‑1 (MQ=255)
                                                          tCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCt                             >  1:348755/1‑96 (MQ=255)
                                                           cATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTg                            <  1:14005686/96‑1 (MQ=255)
                                                            aTTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAAGAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTgg                           >  2:317085/1‑96 (MQ=255)
                                                            aTTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTgg                           >  1:1843847/1‑96 (MQ=255)
                                                            aTTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTgg                           <  2:11601064/96‑1 (MQ=255)
                                                            aTTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGAAATGGTCAGCAGAATCCCCTTATCTTTACCATCTgg                           >  1:7252951/1‑96 (MQ=255)
                                                             tttGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGt                          <  1:822974/96‑1 (MQ=255)
                                                              ttGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTc                         >  2:8371243/1‑96 (MQ=255)
                                                              ttGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTc                         <  2:11027904/96‑1 (MQ=255)
                                                               tGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCa                        <  2:7383219/96‑1 (MQ=255)
                                                               tGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCa                        <  2:6568772/96‑1 (MQ=255)
                                                                 gCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATg                      >  2:1834256/1‑96 (MQ=255)
                                                                 gCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATg                      >  2:13769220/1‑96 (MQ=255)
                                                                   aaTTGAAAAACTGGCCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGAc                    <  1:802810/96‑1 (MQ=255)
                                                                   aaTTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGAc                    >  2:3822790/1‑96 (MQ=255)
                                                                    aTTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACg                   <  1:7119624/96‑1 (MQ=255)
                                                                     ttGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGc                  <  2:13611702/96‑1 (MQ=255)
                                                                      tGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCt                 >  1:3776742/1‑96 (MQ=255)
                                                                       gAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTg                >  2:9219306/1‑96 (MQ=255)
                                                                        ataaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGa               >  1:11091890/3‑96 (MQ=255)
                                                                         aaaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGaa              >  2:3600099/1‑96 (MQ=255)
                                                                         aaaaCTGACCAGCGCCAGCTTTGCTTTTACTGGCGAACAGCCGCAGCAATAGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGaa              <  2:8211687/96‑1 (MQ=255)
                                                                          aaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCATCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGaaa             >  1:7752140/1‑96 (MQ=255)
                                                                          aaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGaaa             >  2:940813/1‑96 (MQ=255)
                                                                           aaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAg            <  2:4598074/96‑1 (MQ=255)
                                                                             cTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGAc          <  2:584663/96‑1 (MQ=255)
                                                                              tGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACg         <  2:148589/96‑1 (MQ=255)
                                                                                aCCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGcc       >  2:6082593/1‑96 (MQ=255)
                                                                                 ccAGCGCCAGCTTTGCTTTTACTGTCGATCAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCa      <  1:12986947/96‑1 (MQ=255)
                                                                                  cAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCaa     >  1:10028759/1‑96 (MQ=255)
                                                                                   aGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAAc    <  1:3270073/96‑1 (MQ=255)
                                                                                   aGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAAc    <  2:12990436/96‑1 (MQ=255)
                                                                                    gcgcCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACg   >  1:5008178/1‑96 (MQ=255)
                                                                                     cgcCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACgg  >  2:3436202/1‑96 (MQ=255)
                                                                                          ||                                                                                         
CGACGCTGGAATATACCCTGTTTGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCACGTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGG  >  NC_000913/3221007‑3221187

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: