Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,600,609 (C)5→6 coding (132/471 nt) bcp → peroxiredoxin, thiol peroxidase, thioredoxin‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,600,6041.C100.0% 121.1 / NA 35T43P (ACC→CCC) bcpperoxiredoxin, thiol peroxidase, thioredoxin‑dependent
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (16/19);  total (16/19)

AGTTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAAGAG  >  NC_000913/2600469‑2600732
                                                                                                                                        |                                                                                                                                
aGTTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccc                                                                                                                               <  1:1569413/139‑1 (MQ=255)
aGTTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccc                                                                                                                               >  2:1590992/1‑139 (MQ=255)
     aGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccccccgg                                                                                                                          <  1:1490605/139‑4 (MQ=255)
          tGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTa                                                                                                                      <  1:731920/139‑1 (MQ=255)
          tGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTa                                                                                                                      <  1:368828/139‑1 (MQ=255)
                 aCTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACa                                                                                                               >  2:534060/1‑139 (MQ=255)
                   tGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAgg                                                                                                             >  2:524727/1‑139 (MQ=255)
                            gTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCt                                                                                                    >  1:504718/1‑139 (MQ=255)
                              gATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTa                                                                                                  <  1:1590992/139‑1 (MQ=255)
                                        cGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACa                                                                                        <  1:517379/139‑1 (MQ=255)
                                           aaTTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATgg                                                                                     <  1:754831/139‑1 (MQ=255)
                                                     gCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGGTGAGTTGaa                                                                           <  2:504718/139‑1 (MQ=255)
                                                           tCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGc                                                                     >  2:1782916/1‑139 (MQ=255)
                                                               gACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGc                                                                                                     >  1:711193/1‑103 (MQ=255)
                                                               gACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGc                                                                                                     <  2:711193/103‑1 (MQ=255)
                                                               gACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccc                                                                                                                               <  1:1226631/76‑1 (MQ=255)
                                                               gACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGATCTGGTTTATTTCTACCCGAAAGCCATGAccc                                                                                                                               >  2:1226631/1‑76 (MQ=255)
                                                                    atcaCAAGTTGATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGa                                                            <  1:524727/136‑1 (MQ=255)
                                                                                tttGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTAt                                                <  1:1782916/139‑1 (MQ=255)
                                                                                tttGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTAt                                                >  1:983932/1‑139 (MQ=255)
                                                                                    aCCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTa                                                                                                                      >  2:560623/1‑65 (MQ=255)
                                                                                    aCCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTa                                                                                                                      <  1:560623/65‑1 (MQ=255)
                                                                                      cGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTAt                                                <  2:1771807/133‑1 (MQ=255)
                                                                                      cGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTAt                                                >  1:1771807/1‑133 (MQ=255)
                                                                                          ttCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGAt                                      <  1:534060/139‑1 (MQ=255)
                                                                                                    aGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGaaa                            >  2:1374334/1‑139 (MQ=255)
                                                                                                    aGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTACGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGc                                                       >  1:1672171/1‑112 (MQ=255)
                                                                                                    aGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTACGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGc                                                       <  2:1672171/112‑1 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTa                                                                                                  <  1:1204987/67‑1 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTa                                                                                                  >  2:1204987/1‑67 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGAt                                                                                   <  1:244158/82‑1 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGAt                                                                                   >  2:244158/1‑82 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGAtg                                                          >  2:770123/1‑107 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGAtg                                                          <  1:770123/107‑1 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGaaaaa                          >  1:1760245/1‑139 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccccccgg                                                                                                                          >  2:516551/1‑39 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccccccgg                                                                                                                          <  1:516551/42‑4 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccc                                                                                                                               >  1:925813/1‑37 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccc                                                                                                                               <  2:925813/37‑1 (MQ=255)
                                                                                                                tttATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGtt                >  2:1041998/1‑139 (MQ=255)
                                                                                                                              aaaGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTAcgcg                                                                                              >  2:1751958/1‑47 (MQ=255)
                                                                                                                              aaaGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTAcgcg                                                                                              <  1:1751958/47‑1 (MQ=255)
                                                                                                                              aaaGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAagag  <  2:1697504/139‑1 (MQ=255)
                                                                                                                              aaaGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAagag  <  1:1374334/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                
AGTTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAAGAG  >  NC_000913/2600469‑2600732

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: