Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,224,192 A→G 5.4% G520G (GGA→GGG ebgA → evolved beta‑D‑galactosidase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,224,1920AG 5.4% 233.2 / 3.1 74G520G (GGA→GGGebgAevolved beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (27/43);  new base G (2/2);  total (29/45)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

GGTGCATTACGAAGAAGATCGCGATGCTGAAGTGGTCGATATTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAGTGGTGCGACCACGGGATCCAGGCACAGGACGACCACGGCAATGTCTGGTATAAATTCGGCGGCGACTACGG  >  NC_000913/3224045‑3224335
                                                                                                                                                   |                                                                                                                                               
ggTGCATTACGAAGAAGATCGCGATGCTGAAGTGGTCGATATTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGa                                                                                                                                                 >  1:85677/1‑148 (MQ=255)
   gCATTACGAAGAAGATCGCGATGCTGAAGTGGTCGATATTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACcg                                                                                                                                              >  2:307982/1‑148 (MQ=255)
       tACGAAGAAGATCGCGATGCTGAAGTGGTCGATATTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACcgggcg                                                                                                                                          >  7:206800/1‑148 (MQ=255)
         cGAAGAAGATCGCGATGCTGAAGTGGTCGATATTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACcgggcggg                                                                                                                                        <  6:105112/148‑1 (MQ=255)
         cGAAGAAGATCGCGATGCTGAAGTGGTCGATATTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACcgggcgg                                                                                                                                         <  7:50723/147‑1 (MQ=255)
         cGAAGAAGATCGCGATGCTGAAGTGGTCGATATTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACcgggcgg                                                                                                                                         <  1:307982/147‑1 (MQ=255)
            agaagaTCGCGATGCTGAAGTGGTCGATATTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCt                                                                                                                                      >  8:239351/1‑147 (MQ=255)
            agaagaTCGCGATGCTGAAGTGGTCGATATTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCt                                                                                                                                      >  4:292268/1‑147 (MQ=255)
            agaacaCCGCTATGCTCAAGTGGTCGATATTATTTCCACCATCTACACCCGCGAGCCGCTGAGGACTGAGTTTGGTGAATACCCGCATCCGAATCCGCGCATCATCTGTGAATATGCTCATCCGATGGGGAACGGACCGGGCGGGCTg                                                                                                                                     <  4:106375/148‑1 (MQ=255)
             gaagaTCGCGATGCTGAAGTGGTCGATATTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGGGTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATCTGCTCATGCGATGGGGAACGGACCGGGCGGGCTGa                                                                                                                                    <  4:251322/148‑1 (MQ=255)
               agaTCGCGATGCTGAAGTGGTCGATATTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCACATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGAc                                                                                                                                   <  6:227683/147‑1 (MQ=255)
                gaTCGCGATGCTGAAGTGGTCGATATTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACgg                                                                                                                                 <  7:367615/148‑1 (MQ=255)
                    gcgATGCTGAAGTGGTCGATATTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGt                                                                                                                              <  6:423044/147‑1 (MQ=255)
                    gcgATGCTGAAGTGGTCGATATTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGAGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTa                                                                                                                             >  3:33483/1‑148 (MQ=255)
                      gATGCTGAAGTGGTCGATATTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTAc                                                                                                                            >  1:189830/1‑147 (MQ=255)
                              aGTGGTCGATATTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACgg                                                                                                                                 <  2:256503/134‑1 (MQ=255)
                              aGTGGTCGATATTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACgg                                                                                                                                 >  1:256503/1‑134 (MQ=255)
                               gTGGTCGATATTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGtct                                                                                                                  <  3:61973/148‑1 (MQ=255)
                                   tCGATATTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCta                                                                                                              <  7:29701/148‑1 (MQ=255)
                                   tCGATATTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCta                                                                                                              <  4:14780/148‑1 (MQ=255)
                                    cGATATTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCta                                                                                                              >  1:257143/1‑147 (MQ=255)
                                     gATATTATTTCCACCATGTACACCCGCTTGCCGCTGATTAATGAGTTTGTTGAATAGCCGCATCCGAAGCCGAGCATCATCTGTTAATTTGCTGATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTAGCAGAACGTCTTCtata                                                                                                            <  2:96436/148‑1 (MQ=255)
                                      atatTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACgg                                                                                                                                 <  4:184119/126‑1 (MQ=255)
                                      atatTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACgg                                                                                                                                 >  3:184119/1‑126 (MQ=255)
                                      atatTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCtat                                                                                                             >  1:13485/1‑146 (MQ=255)
                                       tattatTTCCACCATGTACACCCGCGTGCCGATGATTAATGAGTTTGGTGATTACCCGCATCCGAAGCCGGGCTTCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATaa                                                                                                           <  7:99389/147‑1 (MQ=255)
                                         ttattTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGc                                                                                                         <  2:85677/147‑1 (MQ=255)
                                            tttCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGa                                                                                                     <  7:239351/148‑1 (MQ=255)
                                              tCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGAt                                                                                                    <  2:110168/147‑1 (MQ=255)
                                                caccaTGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGc                                                                                                 <  7:17922/148‑1 (MQ=255)
                                                caccaTGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGc                                                                                                 <  2:7262/148‑1 (MQ=255)
                                                  ccaTGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCAt                                                                                               >  7:434261/1‑148 (MQ=255)
                                                   caTGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACAGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAATCAAGATTGCag                                                                                               >  4:325031/1‑146 (MQ=255)
                                                     tGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCa                                                                                            <  7:37668/148‑1 (MQ=255)
                                                     tGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCa                                                                                            <  7:253322/148‑1 (MQ=255)
                                                       tACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAg                                                                                           <  3:346293/147‑1 (MQ=255)
                                                         cacCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCtata                                                                                                            >  7:281720/1‑128 (MQ=255)
                                                         cacCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCtata                                                                                                            <  8:281720/128‑1 (MQ=255)
                                                             cgcgTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCAtt                                                                                    <  8:103743/148‑1 (MQ=255)
                                                                 tGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATGAGCACGATTGCATTCAGGGTCATTATGt                                                                                >  4:154597/1‑148 (MQ=255)
                                                                 tGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGt                                                                                <  8:434261/148‑1 (MQ=255)
                                                                  gCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTc                                                                               >  8:235656/1‑148 (MQ=255)
                                                                   ccGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATAATCTCTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCt                                                                              <  8:87477/148‑1 (MQ=255)
                                                                      cTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTggg                                                                           <  3:356879/148‑1 (MQ=255)
                                                                      cTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTggg                                                                           <  7:235656/148‑1 (MQ=255)
                                                                                aGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAGTGGTGCg                                                                  <  2:189830/147‑1 (MQ=255)
                                                                                  tttGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAGTGGTGCGAc                                                                >  7:89555/1‑147 (MQ=255)
                                                                                         aaTACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAGTGGTGCGACCACGGGAt                                                        <  5:78645/148‑1 (MQ=255)
                                                                                           tACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAgtg                                                                       <  3:174998/131‑1 (MQ=255)
                                                                                           tACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAgtg                                                                       >  4:174998/1‑131 (MQ=255)
                                                                                             cccGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAGTGGTGCGACCACGGGATCCa                                                     <  3:292268/147‑1 (MQ=255)
                                                                                               cGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCCGGGTCATTATGTCTGGGCGTGGTGCGACCCCGGGATCCAgg                                                   >  3:215944/1‑147 (MQ=255)
                                                                                                               catcatCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAGTGGTGCGACCACGGGATCCAGGCACAGGACGACCACgg                                   >  6:412797/1‑147 (MQ=255)
                                                                                                                 tcatcTGTGAATATGCTCATGCGATGGGGAACGGGCCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAg                                                                         >  2:215317/1‑107 (MQ=255)
                                                                                                                 tcatcTGTGAATATGCTCATGCGATGGGGAACGGGCCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAg                                                                         <  1:209264/107‑1 (MQ=255)
                                                                                                                 tcatcTGTGAATATGCTCATGCGATGGGGAACGGGCCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAg                                                                         <  1:215317/107‑1 (MQ=255)
                                                                                                                 tcatcTGTGAATATGCTCATGCGATGGGGAACGGGCCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAg                                                                         >  2:209264/1‑107 (MQ=255)
                                                                                                                 tcatcTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCATAACGGCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAGTTTTGCGACCACGGGATCCAGGCACAGTACGACCACGGCaa                                <  4:33483/148‑1 (MQ=255)
                                                                                                                  catcTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAATCACGATTGCATTAAGGGTCATTATGTCTGGGAGTGGTGCGACCACTGGATCCAGGCACAGGACGACCACGGCaa                                >  7:379882/1‑147 (MQ=255)
                                                                                                                  catcTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCt                                                                              <  4:154744/101‑1 (MQ=255)
                                                                                                                  catcTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCt                                                                              >  3:154744/1‑101 (MQ=255)
                                                                                                                  catcTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAGTGGTGCGACCACGGGATCCAGGCACAGGACGACCACGGCAAt                               >  4:211422/1‑148 (MQ=255)
                                                                                                                   atcTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAGTGGTGCGACCACGGGATCCAGGCACAGGACGACCACGGCAATg                              >  7:190099/1‑148 (MQ=255)
                                                                                                                     cTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAGTGGTGCGACCACGGGATCCAGGCACAGGACGACCACGGCAATGTc                            >  1:9196/1‑148 (MQ=255)
                                                                                                                        tgAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAGTGGTGCGACCACGGGATCCAGGCACAGGACGACCACGGCAATGTCt                           <  2:257143/146‑1 (MQ=255)
                                                                                                                        tgAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAGTGGTGCGACCACGGGATCCAGGCACAGGACGACCACGGCAATGTCt                           <  8:206800/146‑1 (MQ=255)
                                                                                                                         gAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAGTGGTGCGACCACGGGATCCAGGCACAGGACGACCACGGCAATGTCTgg                         >  3:274149/1‑147 (MQ=255)
                                                                                                                          aaTATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAGTGGTGCGACCACGGGATCCAGGCACAGGACGACCACGGCAATGTCTGGt                        <  8:190099/147‑1 (MQ=255)
                                                                                                                             atGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCtat                                                                                                             <  7:332111/59‑1 (MQ=255)
                                                                                                                             atGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCtat                                                                                                             >  8:332111/1‑59 (MQ=255)
                                                                                                                                 tCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAGTGGTGCGACCACGGGATCCAGGCACAGGACGACCACGGCAATGTCTGGTATAAATTc                <  7:190225/148‑1 (MQ=255)
                                                                                                                                      cGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAGTGGTGCGACCACGGGATCCAGGCACAGGACGACCACGGCAATGTCTGGTATAAATTcggcgg           <  6:428434/148‑1 (MQ=255)
                                                                                                                                        aTGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAGTGGTGCGACCACGGGATCCAGGCACAGGACGACCACGGCAATGTc                            <  3:163184/129‑1 (MQ=255)
                                                                                                                                        aTGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAGTGGTGCGACCACGGGATCCAGGCACAGGACGACCACGGCAATGTc                            >  4:163184/1‑129 (MQ=255)
                                                                                                                                               aCGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAGTGGTGCGACCACGGGATCCAGGCACAGGACGACCACGGCAATGTCTGGTATAAATTCGGcggcgactacgg  <  2:13485/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GGTGCATTACGAAGAAGATCGCGATGCTGAAGTGGTCGATATTATTTCCACCATGTACACCCGCGTGCCGCTGATGAATGAGTTTGGTGAATACCCGCATCCGAAGCCGCGCATCATCTGTGAATATGCTCATGCGATGGGGAACGGACCGGGCGGGCTGACGGAGTACCAGAACGTCTTCTATAAGCACGATTGCATTCAGGGTCATTATGTCTGGGAGTGGTGCGACCACGGGATCCAGGCACAGGACGACCACGGCAATGTCTGGTATAAATTCGGCGGCGACTACGG  >  NC_000913/3224045‑3224335

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: