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% 135.3 / NA 39T43P (ACC→CCC) bcpperoxiredoxin, thiol peroxidase, thioredoxin‑dependent
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (17/22);  total (17/22)

AGTTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAAGAGCTGCT  >  NC_000913/2600469‑2600737
                                                                                                                                        |                                                                                                                                     
aGTTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccc                                                                                                                                    <  2:119268/139‑1 (MQ=255)
aGTTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccc                                                                                                                                    <  2:1001395/139‑1 (MQ=255)
     aGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccccccgg                                                                                                                               <  2:985242/139‑4 (MQ=255)
     aGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccccccgg                                                                                                                               <  2:910158/139‑4 (MQ=255)
     aGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccccccgg                                                                                                                               <  2:379172/139‑4 (MQ=255)
        aaTGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTg                                                                                                                             <  2:659926/139‑1 (MQ=255)
         aTGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGt                                                                                                                            >  2:88208/1‑139 (MQ=255)
          tGAATCCACTGAAAGCGGGTAATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACATCCAGGGACAGCGTGTTCTGGTTTATTTCTGCCCGAGTGCCTTGACCCCCCGGCTGTa                                                                                                                           <  1:730522/139‑1 (MQ=255)
                     aaaGCCGGTGATATCGCACCGAAATTTCGCTTGCCGGATCAAGACGGAGGACAAGTTAATTTGACCGACTTCCAGGGATAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGcc                                                                                                                >  2:314155/1‑139 (MQ=255)
                     aaaGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGcc                                                                                                                <  1:752433/139‑1 (MQ=255)
                                 atCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTAcgc                                                                                                    <  2:17182/139‑1 (MQ=255)
                                    gCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGAt                                                                                                 <  2:69412/139‑1 (MQ=255)
                                           aaTTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTg                                                                                                                             <  1:618816/104‑1 (MQ=255)
                                           aaTTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTg                                                                                                                             >  2:618816/1‑104 (MQ=255)
                                                      ccGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGaaa                                                                               >  1:332116/1‑139 (MQ=255)
                                                              aGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCggg                                                                       <  1:495670/139‑1 (MQ=255)
                                                                aCGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCg                                                                     <  1:380335/139‑1 (MQ=255)
                                                                    agaACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCggg                                                                       <  1:37259/133‑1 (MQ=255)
                                                                    agaACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCggg                                                                       >  2:817844/1‑133 (MQ=255)
                                                                    agaACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCggg                                                                       <  1:817844/133‑1 (MQ=255)
                                                                    agaACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCggg                                                                       >  2:37259/1‑133 (MQ=255)
                                                                      aaCAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGAtg                                                               >  2:832793/1‑139 (MQ=255)
                                                                           gTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTg                                                          <  1:610107/139‑1 (MQ=255)
                                                                                        aCTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccc                                                                                                                                    <  1:1057637/51‑1 (MQ=255)
                                                                                        aCTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccc                                                                                                                                    >  2:1057637/1‑51 (MQ=255)
                                                                                                    aGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCa                                                >  1:1002653/1‑124 (MQ=255)
                                                                                                    aGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCa                                                <  2:1002653/124‑1 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTa                                                                                                                           >  1:157963/1‑47 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTa                                                                                                                           >  2:850817/1‑47 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTa                                                                                                                           <  1:850817/47‑1 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTa                                                                                                                           <  2:157963/47‑1 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGaaaaa                               >  1:143588/1‑139 (MQ=255)
                                                                                                            cTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTa                                                                                                                           >  1:898470/1‑41 (MQ=255)
                                                                                                            cTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTa                                                                                                                           <  2:898470/41‑1 (MQ=255)
                                                                                                                    tttCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTa                                                                                                                      <  2:814526/38‑1 (MQ=255)
                                                                                                                    tttCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTa                                                                                                                      >  1:814526/1‑38 (MQ=255)
                                                                                                                    tttCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGa                                                                                                  >  2:36494/1‑58 (MQ=255)
                                                                                                                    tttCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGa                                                                                                  <  1:36494/58‑1 (MQ=255)
                                                                                                                         aCCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGttt                    <  2:507638/131‑1 (MQ=255)
                                                                                                                         aCCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGttt                    >  1:507638/1‑131 (MQ=255)
                                                                                                                             gAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAaga        <  1:319172/139‑1 (MQ=255)
                                                                                                                              aaaGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTg                                                                  >  2:53109/1‑80 (MQ=255)
                                                                                                                              aaaGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTg                                                                  <  1:53109/80‑1 (MQ=255)
                                                                                                                              aaaGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAagag       >  1:857419/1‑139 (MQ=255)
                                                                                                                              aaaGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAagag       <  1:314155/139‑1 (MQ=255)
                                                                                                                                   cATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAAGAgctgct  <  1:88208/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                     
AGTTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAAGAGCTGCT  >  NC_000913/2600469‑2600737

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