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.C91.2% 109.3 / 3.1 34coding (127/471 nt)bcpperoxiredoxin, thiol peroxidase, thioredoxin‑dependent
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (12/19);  total (13/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑A‑CCCCC‑GGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAAG  >  NC_000913/2600471‑2600730
                                                                                                                                       |                                                                                                                               
ttAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑A‑ccccc                                                                                                                            <  2:121875/139‑1 (MQ=255)
ttAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑A‑ccccc                                                                                                                            <  2:715018/139‑1 (MQ=255)
   aGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑A‑ccccccgg                                                                                                                         >  2:431949/1‑136 (MQ=255)
     tAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑A‑ccccccggct                                                                                                                       <  2:482533/139‑6 (MQ=255)
       aTGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGt                                                                                                                     >  1:755136/1‑139 (MQ=255)
        tGAATCCACTGAAAGCCGGTGATATCGCACCGACATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGTGACAGCGTGTTCTGGTTTATTTCTCCCCTAAATCCATG‑A‑CCCCCAGGCTGTa                                                                                                                    <  1:186380/139‑1 (MQ=255)
        tGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTa                                                                                                                    <  2:133641/139‑1 (MQ=255)
        tGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTa                                                                                                                    <  2:193960/139‑1 (MQ=255)
        tGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTa                                                                                                                    <  2:471149/139‑1 (MQ=255)
        tGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTa                                                                                                                    >  1:471149/1‑139 (MQ=255)
        tGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTa                                                                                                                    <  1:431949/139‑1 (MQ=255)
        tGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTa                                                                                                                    <  1:200511/139‑1 (MQ=255)
            tCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTa                                                                                                                    <  1:453416/135‑1 (MQ=255)
            tCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTa                                                                                                                    >  2:453416/1‑135 (MQ=255)
            tCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTACCGt                                                                                                                >  2:461137/1‑139 (MQ=255)
             ccACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTACCGTa                                                                                                               <  1:189767/139‑1 (MQ=255)
              cACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTACCGTAc                                                                                                              >  2:600616/1‑139 (MQ=255)
                        cGGTGATATCGCAACGAAATTTAGCTTGCCGGATCAAGACGAAGAACAAGCTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAATGCCATG‑A‑CCCCCATCCTGTACCGTACAGGCCTGCgt                                                                                                    >  1:740194/1‑138 (MQ=255)
                          gTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTACCGTACAGGCCTGCGGCt                                                                                                  <  2:755136/139‑1 (MQ=255)
                                    aCCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTACCGTACAGGCCTGCGGCTTACGCGATaa                                                                                        <  1:600616/139‑1 (MQ=255)
                                     ccGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTACCGTACAGGCCTGCGGCTTACGCGATAAc                                                                                       >  2:584573/1‑139 (MQ=255)
                                           tttAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGAt                                                                                 <  2:467425/139‑1 (MQ=255)
                                                    ccGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGaaa                                                                        <  2:182884/139‑1 (MQ=255)
                                                              aCGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCg                                                              >  2:328972/1‑139 (MQ=255)
                                                                       aaGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGc                                                     <  1:486874/139‑1 (MQ=255)
                                                                                   ccGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCa                                         <  2:244883/139‑1 (MQ=255)
                                                                                    cGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTa                                                                                                                    >  2:130123/1‑63 (MQ=255)
                                                                                    cGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTa                                                                                                                    <  1:130123/63‑1 (MQ=255)
                                                                                        ttCCAGGGACAGTGAGTTCAGCATTATTTCTACCCGAAAGCCATGTC‑CCCCC‑GGCTGTTCCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGAt                                    <  1:681866/139‑1 (MQ=255)
                                                                                                 cAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTg                                                                                                                      >  2:300785/1‑48 (MQ=255)
                                                                                                 cAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTg                                                                                                                      <  1:300785/48‑1 (MQ=255)
                                                                                                    cGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTg                                                   <  2:362615/112‑1 (MQ=255)
                                                                                                    cGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑ACCCCCC‑GGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTg                                                   >  1:362615/1‑112 (MQ=255)
                                                                                                                  tttCTACCCGAAAGCTATG‑ACCCCCC‑GGCTGTACCGTACAGGCCTGCGGCTTAcgcg                                                                                            >  1:536901/1‑57 (MQ=255)
                                                                                                                  tttCTACCCGAAAGCCATG‑ACCCCCC‑GCCTGTACCGTACAGGCCTGCGGCTTAcgcg                                                                                            <  2:536901/57‑1 (MQ=255)
                                                                                                                          cGAAAGCCATG‑ACCCCCC‑GGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAag  <  2:460192/139‑1 (MQ=255)
                                                                                                                            aaaGCCATG‑ACCCCCC‑GGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGAtgtg                                                      <  1:236539/85‑1 (MQ=255)
                                                                                                                            aaaGCCATG‑ACCCCCC‑GGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGAtgtg                                                      >  2:236539/1‑85 (MQ=255)
                                                                                                                                       |                                                                                                                               
TTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATG‑A‑CCCCC‑GGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAAG  >  NC_000913/2600471‑2600730

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