Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,823,523:1 +A 5.6% coding (247/1062 nt) recA ← DNA recombination and repair protein; ssDNA‑dependent ATPase; synaptase; ssDNA and dsDNA binding protein; ATP‑dependent homologous DNA strand exchanger; recombinase A; LexA autocleavage cofactor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,823,5231.A 5.6% 357.1 / 17.1 89coding (247/1062 nt)recADNA recombination and repair protein; ssDNA‑dependent ATPase; synaptase; ssDNA and dsDNA binding protein; ATP‑dependent homologous DNA strand exchanger; recombinase A; LexA autocleavage cofactor
Reads supporting (aligned to +/- strand):  ref base . (41/43);  new base A (2/3);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

GTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACG  >  NC_000913/2823377‑2823670
                                                                                                                                                    |                                                                                                                                                   
gTCAGAGATTTCAAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTGTTTATCGTCGCCCAGTTTACGTTCGTTGATTGGGTCCAGCGCGTGTTGAGGTTCGATAAACGCTCAGGTTTTACCTTCACGCTG‑CGCTgc‑gg                                                                                                                                                   <  3:390432/149‑1 (MQ=255)
 tCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTgc‑ggc                                                                                                                                                  >  4:104315/1‑149 (MQ=255)
   acaGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGa                                                                                                                                                >  4:290717/1‑149 (MQ=255)
   acaGATTTACAGTTCCTGAGCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTT‑CGCTGC‑GGCGa                                                                                                                                                <  6:18917/149‑1 (MQ=255)
    caGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGAt                                                                                                                                               <  3:379613/149‑1 (MQ=255)
    caGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGAt                                                                                                                                               <  8:8961/149‑1 (MQ=255)
     aGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGctgacgctgcaggtga                                                                                                                                                >  8:67213/1‑143 (MQ=255)
     aGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATc                                                                                                                                              >  8:333350/1‑149 (MQ=255)
        tttCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCAcc                                                                                                                                           >  5:432538/1‑149 (MQ=255)
          tCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGGACAGCAGGTTGTCGATATAGACGCACAGTTTACGTGCGTAGATTGGGTCGAGCGCGTGTTCAGCTTCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTg                                                                                                                                         >  1:278712/1‑149 (MQ=255)
            cAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGca                                                                                                                                       <  7:333350/149‑1 (MQ=255)
                  cTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGACGCTGCAGGTGATCGCC‑Gcagcgca                                                                                                                                  <  7:67213/149‑3 (MQ=255)
                  cTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGcagcgtca                                                                                                                                 <  7:465129/149‑1 (MQ=255)
                  cTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGcagcgtca                                                                                                                                 >  6:206435/1‑149 (MQ=255)
                        gCCGGTGTCCTGCTGGGAGCACAGCAGGTTGTCGTTATCGACGCGCAGTTTACGTGTGTAGATTGGGTCCACCCCGTTTTCAGGATCGATAAACCCACAGGGATTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTg                                                                                                                            <  4:254923/148‑1 (MQ=255)
                          cGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATACACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGtt                                                                                                                         >  5:5619/1‑149 (MQ=255)
                               tCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTAcc                                                                                                                    >  5:327391/1‑149 (MQ=255)
                               tCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTAc                                                                                                                     <  8:470646/148‑1 (MQ=255)
                                 cGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGCAGGTGATCGCC‑Gcagcgcag                                                                                                                                 >  5:320218/1‑131 (MQ=255)
                                 cGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGCAGGTGATCGCC‑Gcagcgcag                                                                                                                                 <  6:320218/134‑4 (MQ=255)
                                  ggCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGa                                                                                                                 <  1:244912/149‑1 (MQ=255)
                                   gCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCCGCAGCGTCAGCGTGGTTTTACCGGaa                                                                                                                <  6:188893/149‑1 (MQ=255)
                                        gAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTc                                                                                                           >  3:1694/1‑149 (MQ=255)
                                          gCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCg                                                                                                         <  3:104315/149‑1 (MQ=255)
                                            acaGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGt                                                                                                       <  5:206435/149‑1 (MQ=255)
                                              agcagGTTGTCGATATCGACGCCCAGTTTATGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTcc                                                                                                     >  4:53458/1‑149 (MQ=255)
                                              agcagGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTcc                                                                                                     >  3:417598/1‑149 (MQ=255)
                                              agcagGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATCGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTcc                                                                                                     >  6:369731/1‑149 (MQ=255)
                                                 agGTTGGCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTa                                                                                                  <  2:114181/149‑1 (MQ=255)
                                                   gTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCg                                                                                                    >  8:16927/1‑145 (MQ=255)
                                                         gATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTcga                                                                                          <  6:432538/149‑1 (MQ=255)
                                                               gACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGAt                                                                                               >  3:421675/1‑138 (MQ=255)
                                                               gACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGAt                                                                                               <  4:421675/138‑1 (MQ=255)
                                                                  gCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGAtt                                                                                              <  8:444199/136‑1 (MQ=255)
                                                                  gCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGAtt                                                                                              >  7:444199/1‑136 (MQ=255)
                                                                  gCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGc                                                                                 <  8:241849/149‑1 (MQ=255)
                                                                    ccAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGcc                                                                                >  6:209615/1‑148 (MQ=255)
                                                                            tgTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGACGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCa                                                                        <  3:53458/147‑1 (MQ=255)
                                                                            cGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAg                                                                       <  3:187935/149‑1 (MQ=255)
                                                                            cGTGCGTAGATTGGGTCCAGCGCGGGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGGCCGTAGATTTCGACGATACGGCCCATCGGCAg                                                                       <  1:145934/149‑1 (MQ=255)
                                                                             gTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGc                                                                         >  2:90890/1‑146 (MQ=255)
                                                                             gTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGc                                                                         <  1:90890/146‑1 (MQ=255)
                                                                                cGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGacca                                                                   <  3:290717/149‑1 (MQ=255)
                                                                                 gTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGcagcgtca                                                                                                                                 >  4:331640/1‑86 (MQ=255)
                                                                                 gTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGcagcgtca                                                                                                                                 <  3:331640/86‑1 (MQ=255)
                                                                                 gTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGaccac                                                                  <  6:327391/149‑1 (MQ=255)
                                                                                  tAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATa                                                                                     >  6:263444/1‑129 (MQ=255)
                                                                                  tAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATa                                                                                     <  5:263444/129‑1 (MQ=255)
                                                                                  tAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGGCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGaccacc                                                                 >  5:123322/1‑149 (MQ=255)
                                                                                    gATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTg                                                               >  1:101419/1‑149 (MQ=255)
                                                                                    gATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTg                                                               >  6:50016/1‑149 (MQ=255)
                                                                                      ttGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGc                                                              <  4:1694/148‑1 (MQ=255)
                                                                                          gTCCAGCGCGTGTTAAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGcccc                                                           >  3:361485/1‑147 (MQ=255)
                                                                                            ccAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAgc                                                       >  7:357823/1‑149 (MQ=255)
                                                                                                cgcgTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGcgcg                                                    >  2:418965/1‑148 (MQ=255)
                                                                                                cgcgTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGa                                                   <  4:417598/149‑1 (MQ=255)
                                                                                                cgcgGGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCAG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGa                                                   <  2:278712/149‑1 (MQ=255)
                                                                                                 gcgTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGat                                                  <  1:418965/149‑1 (MQ=255)
                                                                                                    tgtTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATc                                               <  7:16927/149‑1 (MQ=255)
                                                                                                          gCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGacc                                                                    <  2:354573/122‑1 (MQ=255)
                                                                                                          gCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGacc                                                                    >  1:354573/1‑122 (MQ=255)
                                                                                                          gCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTg                                          >  7:472727/1‑148 (MQ=255)
                                                                                                           cATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGa                                         >  2:368833/1‑148 (MQ=255)
                                                                                                            aTCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGCTTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGAAATCCACTGa                                         >  2:393689/1‑147 (MQ=255)
                                                                                                            aTCGATAAACCCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATa                                                                                     <  1:317746/103‑1 (MQ=255)
                                                                                                            aTCGATAAACCCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATa                                                                                     >  2:317746/1‑103 (MQ=255)
                                                                                                                aTAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAgcgc                                                     <  8:96867/131‑1 (MQ=255)
                                                                                                                aTAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAgcgc                                                     >  7:96867/1‑131 (MQ=255)
                                                                                                                  aaaCGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCg                                                                                                         <  8:375417/77‑1 (MQ=255)
                                                                                                                  aaaCGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCg                                                                                                         >  7:375417/1‑77 (MQ=255)
                                                                                                                   aaCGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAAc                                 <  5:180156/148‑1 (MQ=255)
                                                                                                                    aCGCTCAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCTATCAGCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAgc                                                       <  4:326878/125‑1 (MQ=255)
                                                                                                                    aCGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACATGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGACCCAAgc                                                       >  3:326878/1‑125 (MQ=255)
                                                                                                                     cGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAAccag                              >  2:344142/1‑147 (MQ=255)
                                                                                                                      gCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGt                             <  5:369731/149‑1 (MQ=255)
                                                                                                                       cacaGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTa                            <  1:31800/149‑1 (MQ=255)
                                                                                                                        acaGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCgg                              >  7:125900/1‑146 (MQ=255)
                                                                                                                        acaGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCgg                              >  3:348807/1‑146 (MQ=255)
                                                                                                                         caGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAg                                      >  7:117315/1‑137 (MQ=255)
                                                                                                                         caGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAg                                      <  8:117315/137‑1 (MQ=255)
                                                                                                                          aGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTagag                         <  5:50016/149‑1 (MQ=255)
                                                                                                                             ttttACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGaaa                                       >  6:260679/1‑132 (MQ=255)
                                                                                                                             ttttACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGaaa                                       <  5:260679/132‑1 (MQ=255)
                                                                                                                                  ccTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGGAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTc                 >  1:326112/1‑149 (MQ=255)
                                                                                                                                   cTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTgg                                                                                                                           >  8:437970/1‑42 (MQ=255)
                                                                                                                                   cTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTgg                                                                                                                           <  7:437970/42‑1 (MQ=255)
                                                                                                                                   cTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTc                 >  7:106565/1‑148 (MQ=255)
                                                                                                                                    ttCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCca               >  2:145134/1‑149 (MQ=255)
                                                                                                                                       acgctgacgctgCAGGTGATCGCC‑GCAGCG‑CAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGAt                       <  7:411014/138‑1 (MQ=37)
                                                                                                                                       acgctgacgctgCAGGTGATCGCC‑GCAGCG‑CAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGAt                       >  8:411014/1‑138 (MQ=37)
                                                                                                                                        cgctg‑cgctgG‑GGCGATCACCTGCTGGGTGAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTTGAGATGGTTTCCACAt            <  4:361485/148‑1 (MQ=255)
                                                                                                                                        cgctg‑cgctgC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATc           <  1:368833/149‑1 (MQ=255)
                                                                                                                                           tg‑cgctgc‑gGCGATCACCTGCAGCGTCAGGGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCAt        >  4:401094/1‑149 (MQ=255)
                                                                                                                                           tg‑cgctgc‑gGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCAt        >  7:154816/1‑149 (MQ=255)
                                                                                                                                           tg‑cgctgc‑gGCGATCACCTGCAGCGTCAGCCTCGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCAGCCCCAAGCGCCATATCCAGTGAAAGCGCCCCGGTAGAGATGGTTTCCACATCCAt        >  2:458380/1‑149 (MQ=255)
                                                                                                                                               gctgc‑gGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAAc                                 >  1:207146/1‑121 (MQ=255)
                                                                                                                                               gctgc‑gGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAAc                                 <  2:207146/121‑1 (MQ=255)
                                                                                                                                                ctgc‑gGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGACCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGaa    >  5:167201/1‑149 (MQ=255)
                                                                                                                                                   c‑ggcgATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACg  <  6:123322/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACG  >  NC_000913/2823377‑2823670

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: ‑