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.C94.4% 122.3 / ‑0.7 36coding (127/471 nt)bcpperoxiredoxin, thiol peroxidase, thioredoxin‑dependent
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (18/16);  total (19/17)
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

TGGAGTTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAAGAGCTGC  >  NC_000913/2600466‑2600736
                                                                                                                                            |                                                                                                                                     
tGGAGTTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGa                                                                                                                                        >  1:1143312/1‑139 (MQ=255)
   aGTTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccc                                                                                                                                    <  2:597672/139‑1 (MQ=255)
    gTTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAcccc                                                                                                                                   <  1:422178/139‑1 (MQ=255)
        aGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGACATTTAGCTTGCCGGGTCAAGACGGAGAACAAGTTAATATGACCGACTTCCAGCGACAGCGCGTTCTGGTTTATTTCTACCCGAAAGCCATGAccccccgg                                                                                                                              >  2:939511/1‑136 (MQ=255)
        aGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccccccgg                                                                                                                              >  1:494344/1‑136 (MQ=255)
        aGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccccccgg                                                                                                                              <  2:495477/139‑4 (MQ=255)
        aGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGAccccccgg                                                                                                                              <  2:582469/139‑4 (MQ=255)
           aaTGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTg                                                                                                                            <  1:252764/139‑1 (MQ=255)
             tGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCAGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTa                                                                                                                          <  1:124984/139‑1 (MQ=255)
                               gTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCt                                                                                                        <  2:1143312/139‑1 (MQ=255)
                                tGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCtt                                                                                                       <  2:669897/139‑1 (MQ=255)
                                tGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCtt                                                                                                       >  1:669897/1‑139 (MQ=255)
                                      cGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGa                                                                                                 <  2:33958/139‑1 (MQ=255)
                                        cACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTa                                                                                                                          >  1:358895/1‑112 (MQ=255)
                                        cACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTa                                                                                                                          <  2:358895/112‑1 (MQ=255)
                                         aCCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGa                                                                                                 <  2:913439/136‑1 (MQ=255)
                                         aCCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGa                                                                                                 >  1:913439/1‑136 (MQ=255)
                                                       tGCCGGATCAAGACGGATAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTATACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGa                                                                                >  1:821147/1‑139 (MQ=255)
                                                        gCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGa                                                                                                 >  2:1122746/1‑121 (MQ=255)
                                                        gCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGa                                                                                                 <  1:1122746/121‑1 (MQ=255)
                                                             aTCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAg                                                                          >  2:49006/1‑139 (MQ=255)
                                                                   aCGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCg                                                                    <  2:427743/139‑1 (MQ=255)
                                                                        gaACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGac                                                               <  2:821147/139‑2 (MQ=255)
                                                                            aaGTTAATTTGACCGACTTCCAGGGGCAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGc                                                           >  1:573784/1‑139 (MQ=255)
                                                                            aaGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCACGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGc                                                           <  1:939511/139‑1 (MQ=255)
                                                                              gTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGa                                                                                                 <  1:750922/99‑1 (MQ=255)
                                                                              gTTAATTTGACCGACTTACAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGa                                                                                                 >  2:750922/1‑99 (MQ=255)
                                                                                        ctgatttCCAGTGAAAGCGTGTTCGGGTTTATTCCTACCCGCAAGCTATGAATCCCCTGCTGAACCGTACACGGATGCGGCTTACGCCATAATATGGATGACTTTCATAAAGCGGGCCTTGATGTGCTGGGTATCAGCa                                               >  2:136029/6‑139 (MQ=255)
                                                                                                      cAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGa                                                                                        <  2:538459/84‑1 (MQ=255)
                                                                                                      cAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGa                                                                                        >  1:538459/1‑84 (MQ=255)
                                                                                                         cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGc                                                           <  1:1016078/110‑1 (MQ=255)
                                                                                                         cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGc                                                           >  2:1016078/1‑110 (MQ=255)
                                                                                                         cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGaaaaa                              >  1:551234/1‑139 (MQ=255)
                                                                                                         cGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGaaaaa                              >  1:1170104/1‑139 (MQ=255)
                                                                                                             ttCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTg                                                                 >  2:539992/1‑100 (MQ=255)
                                                                                                             ttCTGGTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTg                                                                 <  1:539992/100‑1 (MQ=255)
                                                                                                                  gTTTATTTCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGt                     >  2:650417/1‑139 (MQ=255)
                                                                                                                       tttCTACCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGc                >  2:783217/1‑139 (MQ=255)
                                                                                                                            aCCCGAAAGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGaaa           >  2:423349/1‑139 (MQ=255)
                                                                                                                                 aaaGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAagag      <  2:1170104/139‑1 (MQ=255)
                                                                                                                                 aaaGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAagag      >  2:700406/1‑139 (MQ=255)
                                                                                                                                  aaGCCATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGCTAAAAAAAGCGGGCGTTGATGTGCTGGGTACCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAAGAgc     >  1:1168546/1‑139 (MQ=255)
                                                                                                                                     ccATGACCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAATATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAAGAgctgc  <  1:783217/139‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                     
TGGAGTTAAGTAATGAATCCACTGAAAGCCGGTGATATCGCACCGAAATTTAGCTTGCCGGATCAAGACGGAGAACAAGTTAATTTGACCGACTTCCAGGGACAGCGTGTTCTGGTTTATTTCTACCCGAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTCTCCCGTTTTGCGGAAAAAGAGCTGC  >  NC_000913/2600466‑2600736

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: