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

GTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACG  >  NC_000913/2823377‑2823670
                                                                                                                                                    |                                                                                                                                                   
gTCAGAGATTTCAAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTGTTTATCGTCGCCCAGTTTACGTTCGTTGATTGGGTCCAGCGCGTGTTGAGGTTCGATAAACGCTCAGGTTTTACCTTCACGCTGCGCTgcgg                                                                                                                                                   <  3:390432/149‑1 (MQ=255)
 tCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTgcggc                                                                                                                                                  >  4:104315/1‑149 (MQ=255)
   acaGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGa                                                                                                                                                >  4:290717/1‑149 (MQ=255)
   acaGATTTACAGTTCCTGAGCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTTCGCTGCGGCGa                                                                                                                                                <  6:18917/149‑1 (MQ=255)
    caGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGAt                                                                                                                                               <  3:379613/149‑1 (MQ=255)
    caGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGAt                                                                                                                                               <  8:8961/149‑1 (MQ=255)
     aGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGctgacgctgcaggtga                                                                                                                                                >  8:67213/1‑143 (MQ=255)
     aGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATc                                                                                                                                              >  8:333350/1‑149 (MQ=255)
        tttCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCAcc                                                                                                                                           >  5:432538/1‑149 (MQ=255)
          tCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGGACAGCAGGTTGTCGATATAGACGCACAGTTTACGTGCGTAGATTGGGTCGAGCGCGTGTTCAGCTTCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTg                                                                                                                                         >  1:278712/1‑149 (MQ=255)
            cAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGca                                                                                                                                       <  7:333350/149‑1 (MQ=255)
                  cTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGACGCTGCAGGTGATCGCCGcagcgca                                                                                                                                  <  7:67213/149‑3 (MQ=255)
                  cTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGcagcgtca                                                                                                                                 <  7:465129/149‑1 (MQ=255)
                  cTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGcagcgtca                                                                                                                                 >  6:206435/1‑149 (MQ=255)
                        gCCGGTGTCCTGCTGGGAGCACAGCAGGTTGTCGTTATCGACGCGCAGTTTACGTGTGTAGATTGGGTCCACCCCGTTTTCAGGATCGATAAACCCACAGGGATTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTg                                                                                                                            <  4:254923/148‑1 (MQ=255)
                          cGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATACACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGtt                                                                                                                         >  5:5619/1‑149 (MQ=255)
                               tCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTAcc                                                                                                                    >  5:327391/1‑149 (MQ=255)
                               tCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTAc                                                                                                                     <  8:470646/148‑1 (MQ=255)
                                 cGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCAGGTGATCGCCGcagcgcag                                                                                                                                 >  5:320218/1‑131 (MQ=255)
                                 cGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCAGGTGATCGCCGcagcgcag                                                                                                                                 <  6:320218/134‑4 (MQ=255)
                                  ggCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGa                                                                                                                 <  1:244912/149‑1 (MQ=255)
                                   gCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCCGCAGCGTCAGCGTGGTTTTACCGGaa                                                                                                                <  6:188893/149‑1 (MQ=255)
                                        gAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTc                                                                                                           >  3:1694/1‑149 (MQ=255)
                                          gCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCg                                                                                                         <  3:104315/149‑1 (MQ=255)
                                            acaGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGt                                                                                                       <  5:206435/149‑1 (MQ=255)
                                              agcagGTTGTCGATATCGACGCCCAGTTTATGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTcc                                                                                                     >  4:53458/1‑149 (MQ=255)
                                              agcagGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTcc                                                                                                     >  3:417598/1‑149 (MQ=255)
                                              agcagGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATCGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTcc                                                                                                     >  6:369731/1‑149 (MQ=255)
                                                 agGTTGGCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTa                                                                                                  <  2:114181/149‑1 (MQ=255)
                                                   gTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCg                                                                                                    >  8:16927/1‑145 (MQ=255)
                                                         gATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTcga                                                                                          <  6:432538/149‑1 (MQ=255)
                                                               gACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGAt                                                                                               >  3:421675/1‑138 (MQ=255)
                                                               gACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGAt                                                                                               <  4:421675/138‑1 (MQ=255)
                                                                  gCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGAtt                                                                                              <  8:444199/136‑1 (MQ=255)
                                                                  gCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGAtt                                                                                              >  7:444199/1‑136 (MQ=255)
                                                                  gCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGc                                                                                 <  8:241849/149‑1 (MQ=255)
                                                                    ccAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGcc                                                                                >  6:209615/1‑148 (MQ=255)
                                                                            tgTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGACGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCa                                                                        <  3:53458/147‑1 (MQ=255)
                                                                            cGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAg                                                                       <  3:187935/149‑1 (MQ=255)
                                                                            cGTGCGTAGATTGGGTCCAGCGCGGGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGGCCGTAGATTTCGACGATACGGCCCATCGGCAg                                                                       <  1:145934/149‑1 (MQ=255)
                                                                             gTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGc                                                                         >  2:90890/1‑146 (MQ=255)
                                                                             gTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGc                                                                         <  1:90890/146‑1 (MQ=255)
                                                                                cGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGacca                                                                   <  3:290717/149‑1 (MQ=255)
                                                                                 gTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGcagcgtca                                                                                                                                 >  4:331640/1‑86 (MQ=255)
                                                                                 gTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGcagcgtca                                                                                                                                 <  3:331640/86‑1 (MQ=255)
                                                                                 gTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGaccac                                                                  <  6:327391/149‑1 (MQ=255)
                                                                                  tAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATa                                                                                     >  6:263444/1‑129 (MQ=255)
                                                                                  tAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATa                                                                                     <  5:263444/129‑1 (MQ=255)
                                                                                  tAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGGCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGaccacc                                                                 >  5:123322/1‑149 (MQ=255)
                                                                                    gATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTg                                                               >  1:101419/1‑149 (MQ=255)
                                                                                    gATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTg                                                               >  6:50016/1‑149 (MQ=255)
                                                                                      ttGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGc                                                              <  4:1694/148‑1 (MQ=255)
                                                                                          gTCCAGCGCGTGTTAAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGcccc                                                           >  3:361485/1‑147 (MQ=255)
                                                                                            ccAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAgc                                                       >  7:357823/1‑149 (MQ=255)
                                                                                                cgcgTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGcgcg                                                    >  2:418965/1‑148 (MQ=255)
                                                                                                cgcgTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGa                                                   <  4:417598/149‑1 (MQ=255)
                                                                                                cgcgGGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCAGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGa                                                   <  2:278712/149‑1 (MQ=255)
                                                                                                 gcgTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGat                                                  <  1:418965/149‑1 (MQ=255)
                                                                                                    tgtTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATc                                               <  7:16927/149‑1 (MQ=255)
                                                                                                          gCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGacc                                                                    <  2:354573/122‑1 (MQ=255)
                                                                                                          gCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGacc                                                                    >  1:354573/1‑122 (MQ=255)
                                                                                                          gCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTg                                          >  7:472727/1‑148 (MQ=255)
                                                                                                           cATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGa                                         >  2:368833/1‑148 (MQ=255)
                                                                                                            aTCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGCTTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGAAATCCACTGa                                         >  2:393689/1‑147 (MQ=255)
                                                                                                            aTCGATAAACCCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATa                                                                                     <  1:317746/103‑1 (MQ=255)
                                                                                                            aTCGATAAACCCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATa                                                                                     >  2:317746/1‑103 (MQ=255)
                                                                                                                aTAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAgcgc                                                     <  8:96867/131‑1 (MQ=255)
                                                                                                                aTAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAgcgc                                                     >  7:96867/1‑131 (MQ=255)
                                                                                                                  aaaCGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCg                                                                                                         <  8:375417/77‑1 (MQ=255)
                                                                                                                  aaaCGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCg                                                                                                         >  7:375417/1‑77 (MQ=255)
                                                                                                                   aaCGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAAc                                 <  5:180156/148‑1 (MQ=255)
                                                                                                                    aCGCTCAGGTTTTACCTTCACGCTGCGCTGCGGCTATCAGCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAgc                                                       <  4:326878/125‑1 (MQ=255)
                                                                                                                    aCGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACATGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGACCCAAgc                                                       >  3:326878/1‑125 (MQ=255)
                                                                                                                     cGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAAccag                              >  2:344142/1‑147 (MQ=255)
                                                                                                                      gCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGt                             <  5:369731/149‑1 (MQ=255)
                                                                                                                       cacaGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTa                            <  1:31800/149‑1 (MQ=255)
                                                                                                                        acaGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCgg                              >  7:125900/1‑146 (MQ=255)
                                                                                                                        acaGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCgg                              >  3:348807/1‑146 (MQ=255)
                                                                                                                         caGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAg                                      >  7:117315/1‑137 (MQ=255)
                                                                                                                         caGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAg                                      <  8:117315/137‑1 (MQ=255)
                                                                                                                          aGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTagag                         <  5:50016/149‑1 (MQ=255)
                                                                                                                             ttttACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGaaa                                       >  6:260679/1‑132 (MQ=255)
                                                                                                                             ttttACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGaaa                                       <  5:260679/132‑1 (MQ=255)
                                                                                                                                  ccTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGGAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTc                 >  1:326112/1‑149 (MQ=255)
                                                                                                                                   cTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTgg                                                                                                                           >  8:437970/1‑42 (MQ=255)
                                                                                                                                   cTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTgg                                                                                                                           <  7:437970/42‑1 (MQ=255)
                                                                                                                                   cTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTc                 >  7:106565/1‑148 (MQ=255)
                                                                                                                                    ttCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCca               >  2:145134/1‑149 (MQ=255)
                                                                                                                                       acgctgacgctgCAGGTGATCGCCGCAGCGCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGAt                       <  7:411014/138‑1 (MQ=37)
                                                                                                                                       acgctgacgctgCAGGTGATCGCCGCAGCGCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGAt                       >  8:411014/1‑138 (MQ=37)
                                                                                                                                        cgctgcgctgGGGCGATCACCTGCTGGGTGAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTTGAGATGGTTTCCACAt            <  4:361485/148‑1 (MQ=255)
                                                                                                                                        cgctgcgctgCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATc           <  1:368833/149‑1 (MQ=255)
                                                                                                                                           tgcgctgcgGCGATCACCTGCAGCGTCAGGGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCAt        >  4:401094/1‑149 (MQ=255)
                                                                                                                                           tgcgctgcgGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCAt        >  7:154816/1‑149 (MQ=255)
                                                                                                                                           tgcgctgcgGCGATCACCTGCAGCGTCAGCCTCGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCAGCCCCAAGCGCCATATCCAGTGAAAGCGCCCCGGTAGAGATGGTTTCCACATCCAt        >  2:458380/1‑149 (MQ=255)
                                                                                                                                               gctgcgGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAAc                                 >  1:207146/1‑121 (MQ=255)
                                                                                                                                               gctgcgGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAAc                                 <  2:207146/121‑1 (MQ=255)
                                                                                                                                                ctgcgGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGACCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGaa    >  5:167201/1‑149 (MQ=255)
                                                                                                                                                   cggcgATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACg  <  6:123322/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACG  >  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: