Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,600,609 (C)5→7 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% 152.3 / NA 44T43P (ACC→CCC) bcpperoxiredoxin, thiol peroxidase, thioredoxin‑dependent
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (23/21);  total (23/21)
*NC_0009132,600,6042.C100.0% 154.1 / NA 44T43P (ACC→CCC) bcpperoxiredoxin, thiol peroxidase, thioredoxin‑dependent
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (23/21);  total (23/21)

AGTTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑‑CCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAAGA  >  NC_000913/2600469‑2600731
                                                                                                                                        ||                                                                                                                               
aGTTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑‑ccc                                                                                                                              >  2:175494/1‑139 (MQ=255)
     aGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGTGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑‑cccccccg                                                                                                                         >  1:965694/1‑136 (MQ=255)
     aGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑‑cccccccg                                                                                                                         >  1:108658/1‑136 (MQ=255)
     aGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑‑cccccccg                                                                                                                         >  1:1329363/1‑136 (MQ=255)
      gTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑‑cccccccgg                                                                                                                        <  2:1329363/139‑5 (MQ=255)
      gTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑‑cccccccgg                                                                                                                        <  2:108658/139‑5 (MQ=255)
          tGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGGCGTAGAACAAGTTAATTTGACCGACTGCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑‑cccccccggctgt                                                                                                                    <  1:1065834/139‑9 (MQ=255)
            aaTCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTAc                                                                                                                    >  2:363639/1‑139 (MQ=255)
                cACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTa                                                                                                                <  2:135660/139‑1 (MQ=255)
                             tGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCt                                                                                                   <  1:363639/139‑1 (MQ=255)
                               atatCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑‑cccccccgg                                                                                                                        <  2:529618/114‑5 (MQ=255)
                               atatCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑‑cccccccgg                                                                                                                        >  1:529618/1‑110 (MQ=255)
                                   cGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTAcgcg                                                                                             >  2:1368129/1‑139 (MQ=255)
                                    gCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGa                                                                                            <  1:58662/139‑1 (MQ=255)
                                    gCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGa                                                                                            <  1:175494/139‑1 (MQ=255)
                                    gCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGa                                                                                            <  1:1368129/139‑1 (MQ=255)
                                       ccGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTAcgcg                                                                                             <  2:780197/135‑1 (MQ=255)
                                       ccGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTAcgcg                                                                                             >  1:780197/1‑135 (MQ=255)
                                                               gACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑‑cccccccgg                                                                                                                        <  2:427428/82‑5 (MQ=255)
                                                               gACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑‑cccccccgg                                                                                                                        >  1:427428/1‑78 (MQ=255)
                                                                 cGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑‑c                                                                                                                                <  2:9367/72‑1 (MQ=255)
                                                                 cGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑‑c                                                                                                                                >  1:9367/1‑72 (MQ=255)
                                                                    agaACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTg                                                            >  1:802156/1‑139 (MQ=255)
                                                                       aCAAGTTAATTTGACCGACTTCCAGGGATAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGAtg                                                         <  2:965694/139‑1 (MQ=255)
                                                                         aaGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑‑cccc                                                                                                                             <  1:1413412/67‑1 (MQ=255)
                                                                         aaGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑‑cccc                                                                                                                             >  2:1413412/1‑67 (MQ=255)
                                                                          aGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGc                                                      >  1:1558957/1‑139 (MQ=255)
                                                                           gTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCt                                                     >  1:1237594/1‑139 (MQ=255)
                                                                              aaTTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTggg                                                  >  1:224467/1‑139 (MQ=255)
                                                                              aaTTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑‑cccccccgg                                                                                                                        <  1:1404579/67‑5 (MQ=255)
                                                                              aaTTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑‑cccccccgg                                                                                                                        >  2:1404579/1‑63 (MQ=255)
                                                                                        aCTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCt                                                                                                   <  1:1389885/80‑1 (MQ=255)
                                                                                        aCTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCt                                                                                                   >  2:1389885/1‑80 (MQ=255)
                                                                                             cAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATaa                                   <  2:224467/139‑1 (MQ=255)
                                                                                                    aGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCt                                                                                                   >  1:446229/1‑68 (MQ=255)
                                                                                                    aGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCt                                                                                                   <  2:446229/68‑1 (MQ=255)
                                                                                                    aGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑‑ccccc                                                                                                                            <  1:1382519/41‑1 (MQ=255)
                                                                                                    aGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑‑ccccc                                                                                                                            >  2:1382519/1‑41 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCg                                                                                                      <  2:1198333/63‑1 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCg                                                                                                      >  1:1198333/1‑63 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGGGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATaa                                   <  2:220825/130‑1 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGGGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATaa                                   >  1:220825/1‑130 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGc                                                                    <  2:1113751/97‑1 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGc                                                                    >  1:1113751/1‑97 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGAtgtg                                                       >  2:351312/1‑110 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGAtgtg                                                       <  1:351312/110‑1 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGaaaa                          >  2:981766/1‑139 (MQ=255)
                                                                                                      cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGaaaa                          >  1:997823/1‑139 (MQ=255)
                                                                                                       gtgtTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGaaaaa                         <  1:981766/139‑1 (MQ=255)
                                                                                                               gTTTATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCg                 >  2:681411/1‑139 (MQ=255)
                                                                                                                tttATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTa                                                                                                                <  2:1106298/43‑1 (MQ=255)
                                                                                                                tttATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTa                                                                                                                >  1:1106298/1‑43 (MQ=255)
                                                                                                                tttATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTg                                                    <  2:606024/103‑1 (MQ=255)
                                                                                                                tttATTTCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTg                                                    >  1:606024/1‑103 (MQ=255)
                                                                                                                     ttCTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGc           <  2:802156/139‑1 (MQ=255)
                                                                                                                       cTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAAc                                 >  2:1463894/1‑115 (MQ=255)
                                                                                                                       cTACCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAAc                                 <  1:1463894/115‑1 (MQ=255)
                                                                                                                         aCCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTg                                                    >  2:1262996/1‑94 (MQ=255)
                                                                                                                         aCCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTg                                                    <  1:1262996/94‑1 (MQ=255)
                                                                                                                         aCCCGAAAGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGaa       >  2:58813/1‑139 (MQ=255)
                                                                                                                              aaaGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGaa       <  1:908952/134‑1 (MQ=255)
                                                                                                                              aaaGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGaa       >  2:908952/1‑134 (MQ=255)
                                                                                                                              aaaGCCATGACCCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAaga  >  2:259064/1‑139 (MQ=255)
                                                                                                                                        ||                                                                                                                               
AGTTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑‑CCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAAGA  >  NC_000913/2600469‑2600731

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: