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.C97.9% 169.0 / ‑3.4 48coding (127/471 nt)bcpperoxiredoxin, thiol peroxidase, thioredoxin‑dependent
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (23/24);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

AGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAAGA  >  NC_000913/2600474‑2600731
                                                                                                                                   |                                                                                                                                
aGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccccccgg                                                                                                                         <  1:982245/139‑4 (MQ=255)
aGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccccccgg                                                                                                                         <  1:380457/139‑4 (MQ=255)
     tGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTa                                                                                                                     <  1:717427/139‑1 (MQ=255)
         tCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGt                                                                                                                 >  1:73272/1‑139 (MQ=255)
         tCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGt                                                                                                                 >  1:139517/1‑139 (MQ=255)
            aCTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACa                                                                                                              >  1:1687289/1‑139 (MQ=255)
                   gCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTa                                                                                                                     >  1:639745/1‑125 (MQ=255)
                   gCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTa                                                                                                                     <  2:639745/125‑1 (MQ=255)
                   gCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGc                                                                                                       >  1:1505877/1‑139 (MQ=255)
                          atatCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTa                                                                                                                     >  1:871410/1‑118 (MQ=255)
                          atatCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTa                                                                                                                     <  2:871410/118‑1 (MQ=255)
                                cACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATa                                                                                          >  1:1377355/1‑139 (MQ=255)
                                     aaaTTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATg                                                                                     >  2:1364390/1‑139 (MQ=255)
                                        tttAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGAt                                                                                  <  1:446519/139‑1 (MQ=255)
                                          tAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGa                                                                                <  2:139517/139‑1 (MQ=255)
                                          tAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGa                                                                                <  2:73272/139‑1 (MQ=255)
                                           aGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTAcgcg                                                                                             <  2:1348402/125‑1 (MQ=255)
                                           aGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTAcgcg                                                                                             >  1:1348402/1‑125 (MQ=255)
                                                      tCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGc                                                                    <  2:1687289/139‑1 (MQ=255)
                                                      tCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGc                                                                    >  2:115142/1‑139 (MQ=255)
                                                         aGACGGAGAACAAGTTAATTTGACCGACTCCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCACCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCggg                                                                 <  2:1119084/139‑1 (MQ=255)
                                                           aCGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCg                                                               <  2:1377355/139‑1 (MQ=255)
                                                                 aaCAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGAtg                                                         >  1:1541933/1‑139 (MQ=255)
                                                                    aaGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGc                                                      <  2:1541933/139‑1 (MQ=255)
                                                                      gTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGAtgtg                                                       >  1:759058/1‑136 (MQ=255)
                                                                      gTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGAtgtg                                                       <  2:759058/136‑1 (MQ=255)
                                                                           tttGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTAt                                               <  2:1505877/139‑1 (MQ=255)
                                                                           tttGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTAt                                               <  2:838650/139‑1 (MQ=255)
                                                                           tttGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTAt                                               <  1:115142/139‑1 (MQ=255)
                                                                                         aGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAAc                                 <  2:1342677/139‑1 (MQ=255)
                                                                                              cAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACa                                                                                                              >  2:1685665/1‑57 (MQ=255)
                                                                                              cAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACa                                                                                                              <  1:1685665/57‑1 (MQ=255)
                                                                                               aGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATgg                                                                                    >  1:808217/1‑82 (MQ=255)
                                                                                               aGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATgg                                                                                    <  2:808217/82‑1 (MQ=255)
                                                                                               aGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATg                                                                                 <  2:1601476/85‑1 (MQ=255)
                                                                                               aGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATg                                                                                 >  1:1601476/1‑85 (MQ=255)
                                                                                               aGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTg                                                            <  2:1022364/106‑1 (MQ=255)
                                                                                               aGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTg                                                            >  1:1022364/1‑106 (MQ=255)
                                                                                               aGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGaaa                           >  2:1542770/1‑139 (MQ=255)
                                                                                                 cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGAtgtg                                                       <  1:72997/109‑1 (MQ=255)
                                                                                                 cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGAtgtg                                                       >  2:72997/1‑109 (MQ=255)
                                                                                                 cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccccc                                                                                                                            <  2:1358402/39‑1 (MQ=255)
                                                                                                 cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccccc                                                                                                                            >  1:1358402/1‑39 (MQ=255)
                                                                                                     ttCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAActct                     >  1:1247280/1‑139 (MQ=255)
                                                                                                          gTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTg                                                            >  2:1419204/1‑95 (MQ=255)
                                                                                                          gTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTg                                                            <  1:1419204/95‑1 (MQ=255)
                                                                                                           tttATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGAt                                                                                  >  2:1665485/1‑72 (MQ=255)
                                                                                                           tttATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGAt                                                                                  <  1:1665485/72‑1 (MQ=255)
                                                                                                                    aCCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGaaa      >  1:880958/1‑139 (MQ=255)
                                                                                                                     cccGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGaaaa     >  2:896403/1‑139 (MQ=255)
                                                                                                                       cGAAAGCCATTACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAag   <  1:1542770/139‑1 (MQ=255)
                                                                                                                        gAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAaga  <  2:880958/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                
AGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAAGA  >  NC_000913/2600474‑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: