Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,823,526 C→T 5.2% A82T (GCC→ACC)  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,5260CT 5.2% 302.2 / 4.3 95A82T (GCC→ACC) 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 C (46/44);  new base T (2/3);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.74e-01

TCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTC  >  NC_000913/2823378‑2823673
                                                                                                                                                      |                                                                                                                                                   
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)
    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                                                                                                                                    >  6:206435/1‑149 (MQ=255)
                 cTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGcagcgtca                                                                                                                                    <  7:465129/149‑1 (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                                                                                                                                    <  6:320218/134‑4 (MQ=255)
                                cGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGCAGGTGATCGCC‑Gcagcgcag                                                                                                                                    >  5:320218/1‑131 (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                                                                                                 >  7:444199/1‑136 (MQ=255)
                                                                 gCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGAtt                                                                                                 <  8:444199/136‑1 (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                                                                            <  1:90890/146‑1 (MQ=255)
                                                                            gTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGc                                                                            >  2:90890/1‑146 (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                                                                                        <  5:263444/129‑1 (MQ=255)
                                                                                 tAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATa                                                                                        >  6:263444/1‑129 (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                                                                                        >  2:317746/1‑103 (MQ=255)
                                                                                                           aTCGATAAACCCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATa                                                                                        <  1:317746/103‑1 (MQ=255)
                                                                                                               aTAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAgcgc                                                        >  7:96867/1‑131 (MQ=255)
                                                                                                               aTAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAgcgc                                                        <  8:96867/131‑1 (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                                         <  8:117315/137‑1 (MQ=255)
                                                                                                                        caGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAg                                         >  7:117315/1‑137 (MQ=255)
                                                                                                                         aGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTagag                            <  5:50016/149‑1 (MQ=255)
                                                                                                                            ttttACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGaaa                                          <  5:260679/132‑1 (MQ=255)
                                                                                                                            ttttACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGaaa                                          >  6:260679/1‑132 (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)
                                                                                                                                                     gcgATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCAGATGCATGGAACGGTc  >  6:373452/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
TCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTC  >  NC_000913/2823378‑2823673

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