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.0% 110.7 / ‑3.0 33coding (127/471 nt)bcpperoxiredoxin, thiol peroxidase, thioredoxin‑dependent
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (10/22);  total (10/23)
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.34e-01

TGGAGTTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑CCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAAGAG  >  NC_000913/2600466‑2600732
                                                                                                                                           |                                                                                                                                 
tGGAGTTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGa                                                                                                                                    >  1:265164/1‑139 (MQ=255)
        aGTAATGAATCCACTGAAAGCCGGTGATATCGCGCCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑ccc‑cccgg                                                                                                                          >  1:54694/1‑136 (MQ=255)
        aGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑ccc‑cccgg                                                                                                                          <  1:53430/139‑4 (MQ=255)
             tGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTa                                                                                                                      <  1:200909/139‑1 (MQ=255)
                    aCTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACa                                                                                                               >  1:770341/1‑139 (MQ=255)
                     cTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑ccc‑cccgg                                                                                                                          >  2:77095/1‑123 (MQ=255)
                     cTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑ccc‑cccgg                                                                                                                          <  1:77095/126‑4 (MQ=255)
                               gTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCt                                                                                                    <  1:39231/139‑1 (MQ=255)
                               gTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCt                                                                                                    <  1:773742/139‑1 (MQ=255)
                                tGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCtt                                                                                                   >  1:778581/1‑139 (MQ=255)
                                 gATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTa                                                                                                  <  2:578114/139‑1 (MQ=255)
                                 gATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTa                                                                                                  <  1:798351/139‑1 (MQ=255)
                                      cGCACCGAAATTTAGCTTGCTGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCATGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGa                                                                                             <  2:99823/139‑1 (MQ=255)
                                      cGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGa                                                                                             <  2:265164/139‑1 (MQ=255)
                                             aaaTTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATg                                                                                      >  2:530982/1‑139 (MQ=255)
                                              aaTTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATgg                                                                                     <  1:673825/139‑1 (MQ=255)
                                                tttAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑CCCTCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGAt                                                                                   <  2:785301/139‑1 (MQ=255)
                                                                 aGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCggg                                                                  <  1:530982/139‑1 (MQ=255)
                                                                      gagaACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTg                                                             <  2:54694/139‑1 (MQ=255)
                                                                       agaATAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGTAAGCCATGA‑ccc‑cccgg                                                                                                                          <  2:636520/76‑4 (MQ=255)
                                                                       agaACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑ccc‑cccgg                                                                                                                          >  1:636520/1‑73 (MQ=255)
                                                                        gaACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGAt                                                           <  2:928107/139‑1 (MQ=255)
                                                                         aaCAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGAtg                                                          <  1:910880/139‑1 (MQ=255)
                                                                         aaCAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGAtg                                                          <  2:952257/139‑1 (MQ=255)
                                                                              gTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTg                                                     <  2:440805/139‑1 (MQ=255)
                                                                                   tttGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTAt                                                <  1:996003/139‑1 (MQ=255)
                                                                                      gACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAg                                             <  2:778581/139‑1 (MQ=255)
                                                                                      gACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAg                                             >  1:928450/1‑139 (MQ=255)
                                                                                         cGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTg                                                             >  1:247434/1‑120 (MQ=255)
                                                                                         cGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTg                                                             <  2:247434/120‑1 (MQ=255)
                                                                                               ccAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATaa                                    <  1:941770/139‑1 (MQ=255)
                                                                                                         cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTa                                                                                                                      >  2:107241/1‑47 (MQ=255)
                                                                                                         cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTa                                                                                                                      <  1:107241/47‑1 (MQ=255)
                                                                                                         cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGAtg                                                          >  2:289520/1‑107 (MQ=255)
                                                                                                         cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGAtg                                                          <  1:289520/107‑1 (MQ=255)
                                                                                                                            actcGAAAGCCATGACCCC‑CCGTCTGTACCTTACAGGCCTGGGGCTTACGCGATAACATGGATGAGTTGAACAAAGCGGGCGTTGCTGTGGTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGaaa       >  2:1069936/4‑139 (MQ=255)
                                                                                                                            aCCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGaaa       >  2:133111/1‑139 (MQ=255)
                                                                                                                            aCCCGAAAGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGaaa       >  1:574992/1‑139 (MQ=255)
                                                                                                                                 aaaGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAagag  <  2:917747/139‑1 (MQ=255)
                                                                                                                                 aaaGCCATGACCCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAagag  <  1:133111/139‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                 
TGGAGTTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGA‑CCC‑CCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAAGAG  >  NC_000913/2600466‑2600732

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