Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,084,279 A→C 34.4% N73T (AAC→ACC)  sodB_2 → Superoxide dismutase [Mn/Fe]

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,084,2790AC34.4% 52.2 / 5.9 48N73T (AAC→ACC) sodB_2Superoxide dismutase [Mn/Fe]
Reads supporting (aligned to +/- strand):  ref base A (20/10);  new base C (2/16);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 2.42e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCATCACCAGACTTACGTGAACGGCCTCAATGCCGCCCTCGAGGGTACCGAGTGGGCCGAATGGCCGGTGGAAAAACTGGTCGCTGCCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCGTGACCAACCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTGGCGAGCCTGTGGGCCAGCTGGCGCAAGCCGTCACTGCACAACTTGGCGGCTTCGAGGCTTTCAAGGAGG  >  ACB122/1084150‑1084410
                                                                                                                                 |                                                                                                                                   
aGCATCACCAGACTTACGTCACCGGCCTCAAGGCCGCCCTCGAGGGTACCGAGTGGGCCGAATGGCCGGTGGAAAAACTGGTCGCTGCCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCGTGACCAACCAcg                                                                                                                                <  7:55862/135‑1 (MQ=255)
   aTCACCAGACTTACGTGAACGGCCTCAATGCCGCCCTCGAGGGTACCGAGTGGGCCGAATGGCCGGTGGAAAAACTGGTCGCTTCCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCGTGACCAACCAccgcg                                                                                                                             >  6:230269/1‑135 (MQ=255)
   aTCACCAGACTTACGTGAACGGCCTCAATGCCGCCCTCGAGGGTACCGAGTGGGCCGAATGGCCGGTGGAAAAACTGGTCGCTGCCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCGTGACCAACCAcggcg                                                                                                                             >  6:124859/1‑135 (MQ=255)
         aGACTTACGTGAACGGCCTCAATGCCGCCCTCGAGGGTACCGAGTGGGCCGAATGGCCGGTGGAAAAACTGGTCGCTGCCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCGTGACCAACCACGGCGGTGGTc                                                                                                                       >  6:23165/1‑135 (MQ=255)
           aCTTACGTGACCGGCCTCAAGCCCGCCCGCGAGGGCACCGAGTGGGCCGAATGGCCGGTGTAAAAACTGGGCGCTGCCGTCAAGCAACTGCCGCAAAACCTGCGGGGCGCCGTCACCACCCACGGCGGTGGTCAc                                                                                                                     <  5:348152/135‑1 (MQ=255)
              tACGTGAACGGCCTCAATGCCGCCCTCGAGGGTACCGAGTGGGCCGAATGGCCGGTGGAAAAACTGGTCGCTGCCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCGTGACCACCCACGGCGGTGGTCACGc                                                                                                                   >  6:380945/1‑134 (MQ=255)
               aCGTGAACGGCCTCAATGCCGCCCTCGAGGGTACCGAGTGGGCCGAATGGCCGGTGGAAAAACTGGTCGCTGCCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCGTGACCAACCACGGCGGTGGTCACGCCa                                                                                                                 >  7:46060/1‑135 (MQ=255)
                cGTGACCGCCCTCAAGGCCGCCCTGGAGGGTACCGAGGGGGCCGAAGGGCCGGTGTAAAAACTGGTCGCTGCCGGCAAGCAACTGCCGGAAAACCTGCGGGGCGCCGTGACCACCCACGGCGGTGGTCACGCCa                                                                                                                 <  4:265354/134‑1 (MQ=255)
                  tGAACGGCCTCAATGCCGCCCTCGAGGGTACCGAGTGGGCCGAATGGCCGGTGGAAAAACTGGTCGCTGCCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCGTGACCAACCACGGCGGTGGTCACGCCAAcc                                                                                                              >  5:189927/1‑135 (MQ=255)
                    cccGGCCTCAAGGCCGCCCTCGAGGGTCCCGAGTGGGCCTAAGGGCCGGTGAAAAACCTGGGCGCGGCCGACACGCAACTCCCGTAAACCCGGCGGGGCGCCGTCACCAACCACGGCGGTGGTCACGCCAACcac                                                                                                            <  6:268445/133‑1 (MQ=255)
                    caCGGCCGAAATGCCGCCGTGGTGGGTACGGAGTGGTCCGACTGGCGGGTGTAAAAACTGGTCGGGGCCGGCAACCAACTGCGGCAACCCCTGCGGGGCGCCGTTACCACCCACGGCGGTGGTCACGCCAACcac                                                                                                            <  6:122188/134‑1 (MQ=255)
                    aaCGGCCTCAATGCCGCCCTCGAGGGTACCTAGTGGGCCTAATGGCCGGTGTAAAAACTGGTCGCTGCCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCGTGACCAACCACGGCGGTGGTCACGCCAACcac                                                                                                            <  7:275729/135‑1 (MQ=255)
                     ccGGCCTCAACGCCGCCCTCGCGGGACCCGAGGGGGCCGAAGGGCCGGTGCAAAACCGGGGCGCTGCCGACAAGCAACTGCCGTAAAACCTGCGTGGCGCCGTCACCACCCACGGCGGTGGTCACGCCAACcac                                                                                                            <  5:23165/133‑1 (MQ=255)
                              aTGCCGCCCTCTAGGTTCCCGAGTGGGCCTAATGGCCGGTGTAAAAACTGGGCGCTGCCTTCAAGCAACTGCCTAAAAACCTGCGGGGCGCCGTGACCAACCACGGCGGTGGTCACGCCAACCACAGCCTGTTCt                                                                                                  <  5:176275/135‑1 (MQ=255)
                              aTGCCGCCCTCGAGGGATCGGCGGGGGCCGAATGGCCGTTGTACAAACTGGTCGCTGCCTTCAAGCAACTGCCGTAAAACCTGCGTGGCGCCGTGACCAACCACGGCGGTGGTCACGCCAACCACAGCCTGTTCt                                                                                                  <  2:121948/135‑1 (MQ=255)
                                      cTCGAGGGTACCGAGTGGGCCGAATGGCCGGTGGAAAAACTGGTCGCTGCCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCGTGACCAACCACGGCGGTGGTCACGCCAACCACAGCCTCGTCTGGACTGTc                                                                                          >  6:352630/1‑135 (MQ=255)
                                      cTCGAGGGTACCGAGGGGGCCGAATGGCCGGTGGAAAACCTGGGCGCTCGCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCGTGACCAACCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTc                                                                                          <  8:75307/135‑1 (MQ=255)
                                       tCGCGGGTACCGAGTGGGCCGACTGGCCGGTGTAAAAACTGGTCGCTGCCGTCAACCAACTGCCGCAAAACCGGCGGGGCGCCGTGACCACCCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCa                                                                                         <  3:158322/135‑1 (MQ=255)
                                       tCGAGGGTACCGAGTGGGCCGAATGGCCGGTGGACAAACTGGTCCCTGCCGTCAAGCAACAGCCGGAAAACATGCGAGACGCCGTGACCAACCCCCGCGGTGGTCCCGCCACCCACCGCCTGTGCGGGAGTGTcc                                                                                         >  4:18962/1‑134 (MQ=255)
                                       tCGAGGGTACCGAGTGGGCCGAATGGCCGGTGGAAAAACTGGTCGCTGCCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCGTGACCAACCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCa                                                                                         >  7:130470/1‑135 (MQ=255)
                                            ggcaCCGAGTGGGCCGACGGGCCGGTGGAAAAACTGGGCGCCGCCGCCAACAAACTGCCGCAAACCCTGCGGGGCGCCGTCACCACCCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCg                                                                                    <  5:68464/132‑1 (MQ=255)
                                                 cGAGTGGGCCGAATGGCCGGTGGAAAAACTGGTCGCTGCCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCGTGACCAACCACGGCGGTGGTCACGCCAACCCCAGCCGGGTCTGGACTGTCAAGTCGCCTcc                                                                               >  5:367361/1‑134 (MQ=255)
                                                    gTGGGCCGAATGGCCGGTGGAAAAACTGGTCGCTGCCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCCTGACCAAACACCGCCGTGGTCACGCCCACCACACCCGGTTCAGGACAGTCAAGTAGCCTCAcgg                                                                            >  4:208422/1‑135 (MQ=255)
                                                      gggCCGAATGGCCGGTGGAAAAACTGGTCGCTGCCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCGTGACCAACCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCg                                                                          >  7:154747/1‑135 (MQ=255)
                                                        gCCGAATGGCCGGTGGAAAAACTGGTCGCTGCCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCGTGACCAACCACGGCGGTGGTCACGCCAACCACAGCCTGGTCTGGCCTGTCATGTCGCCTCAGGGCGGt                                                                        >  1:160897/1‑135 (MQ=255)
                                                        gCCGAAGGGCGGGTGGAAAACGTGGGCGCTGCCGTCAACCAACTGCCGTAAACCCTGCGGGGCGCCGTGACCAACCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGt                                                                        <  8:130470/135‑1 (MQ=255)
                                                        gCCGAAGGGCGGGTGGAAAAACTGGTCGCTGCCGGCAAGCAACGGCCGGAAAACCTGCGTGGCGCCGTGACCAACCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGt                                                                        <  7:78700/135‑1 (MQ=255)
                                                         ccGAGTGGCCGGTGGAAAACCTGGTCGCTGCCGTCAAGCAACTGCCGGAAACCCTGCGTGGCGCCGTGACCACCCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTg                                                                       <  8:5690/135‑1 (MQ=255)
                                                          cGAATGGCCGGTGGAAAAACTGGTCGCTGCCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCGTGACCAACCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTgg                                                                      >  7:178399/1‑135 (MQ=255)
                                                                gCCGGTGGAAAAACTGGTCGCTGCCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCGTGACCAACCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTGGCGAGcc                                                                >  1:293256/1‑135 (MQ=255)
                                                                    gTGGAAAAACTGGTCGCTGCCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCGTGACCAACCACGGCGGTGGTCCCGCCCCCCACCGCCTGTTCTGGACGGGCATGTCGCCTCAGGGCGGGGGCGCGCCtgtg                                                            >  3:4122/1‑135 (MQ=255)
                                                                    gTGGAAAAACTGGTCGCTGCCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCGTGACCAACCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTGGCGAGccagag                                                            >  6:98326/1‑131 (MQ=255)
                                                                    gTGGAAAAACTGGTCGCTGCCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCGTGACCAACCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTGGCGAGCCtgtg                                                            >  2:16742/1‑135 (MQ=255)
                                                                             cTGGTCGCTGCCGTCAAGCACCTGCCGGAAAACCGGCGTGGCGCGGTGACCACCCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTGGCGAGCCTGTGGGCCAGCTg                                                   <  2:204096/135‑1 (MQ=255)
                                                                             cTGGGCGCTGCCTTCAAGCAACTGCCGTAAAACCTGCGTGGCGCCGTGACCACCCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTGGCGAGCCTGTGGGCCAGCTg                                                   <  2:122247/135‑1 (MQ=255)
                                                                                  cgcgGCCGGCAAGCACCGGCCGCAAACCCGGCGGGGCGCCGTGACCACCCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTGGCGAGCCTGTGGGCCAGCTGGCGCa                                              <  6:215107/135‑1 (MQ=255)
                                                                                  cGCTGCCGTCAAGCAACTGCCGGAAACCCTGCGGGGCGCCGTCACCAACCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTGGCGAGCCTGTGGGCCAGCTGGCGCa                                              <  4:242668/135‑1 (MQ=255)
                                                                                  cGCTGCCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCGTGACCAACCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTGGCGAGCCTGTGGGCCAGCTGgcg                                                >  1:72652/1‑133 (MQ=255)
                                                                                  cGCTGCCGTCAACCAACTGCCGGAAAACCTGCGTGGCGCGGTGACCACCCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTGGCGAGCCTGTGGGCCAGCTGGCGCa                                              <  2:115032/135‑1 (MQ=255)
                                                                                  cGCTGCCGGCAAGCAACTGCCGGAAACCCGGCGGGGCGCCGTGACCACCCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTGGCGAGCCTGTGGGCCAGCTGGCGCa                                              <  3:76522/135‑1 (MQ=255)
                                                                                          tCACGCAACTGGCGGACAACCGGCGTGGCGCCGGGACCACCCACGTCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTGGCGAGCCTGTGGGCCAGCTGGCGCAAGCCGTCa                                      <  7:33903/135‑1 (MQ=255)
                                                                                                     ccGGAAAACCTGCGTGGCGCCGTGACCAACCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTGGCGAGCCTGTGGGCCAGCTGGCGCAAGCCGTCACTGCACAACtt                           >  1:85512/1‑135 (MQ=255)
                                                                                                       ggAAAACCTGCGTGGCGCCGTGACCAACCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTGGCGAGCCTGTGGGCCAGCTGGCGCAAGCCGTCACTGCACAACTTgg                         >  8:57038/1‑135 (MQ=255)
                                                                                                             ccTGCGTGGCGCCGTGACCAACCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTGGCGAGCCTGTGGGCCAGCTGGCGCAAGCCGTCACTGCACAACTTGGCGGCtt                   <  2:79266/135‑1 (MQ=255)
                                                                                                               ggCGTGGCGCCTTGACCACCCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTGGCGAGCCTGTGGGCCAGCTGGCGCAAGCCGTCACTGCACAACTTGGCGGCTTCg                 <  7:81358/134‑1 (MQ=255)
                                                                                                                 cggggCGCGGGGCCCACCCACGGCGGTGGTCACGCCAAACACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTGGCGAGCCTGTGGGCCAGCTGGCGCAAGCCGTCACTGCACAACTTGGCGGCTTCGAg               <  1:262899/132‑1 (MQ=255)
                                                                                                                   tGGCGCCGTGCCCACCCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTGGCGAGCCTGTGGGCCAGCTGGCGCAAGCCGTCACTGCACAACTTGGCGGCTTCGAGGc             >  4:177259/1‑135 (MQ=255)
                                                                                                                       gcCGTCACCACCCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTGGCGAGCCTGTGGGCCAGCTGGCGCAAGCCGTCACTGCACAACTTGGCGGCTTCGAGGCTTTc         <  3:208422/135‑1 (MQ=255)
                                                                                                                             accacccaCGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTGGCGAGCCTGTGGGCCAGCTGGCGCAAGCCGTCACTGCACAACTTGGCGGCTTCGAGGCTTTCAaggag   <  6:360518/135‑1 (MQ=255)
                                                                                                                             accaaccaCGGCGGTGGTCACGCCAACCTCAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTGGCGAGCCTGTGGGCCAGCTGGCGCAAGCCGTCACTGCACAACTTGGCGGCTTCGAGGCTTTCAaggag   <  5:230269/135‑1 (MQ=255)
                                                                                                                              ccaaccaCCCCGGTGGTCACGCCACCCCCAGCCTGTTCTGGACGGTCATGTCGCCTCCGGGCGGTGGCGCGCCTGTGGGCCAGCTGGCGCAAGCCGTCACTGCACAACTTGGCGGCGGCGAGGCTTTCAaggagg  >  8:174800/1‑135 (MQ=255)
                                                                                                                                 |                                                                                                                                   
AGCATCACCAGACTTACGTGAACGGCCTCAATGCCGCCCTCGAGGGTACCGAGTGGGCCGAATGGCCGGTGGAAAAACTGGTCGCTGCCGTCAAGCAACTGCCGGAAAACCTGCGTGGCGCCGTGACCAACCACGGCGGTGGTCACGCCAACCACAGCCTGTTCTGGACTGTCATGTCGCCTCAGGGCGGTGGCGAGCCTGTGGGCCAGCTGGCGCAAGCCGTCACTGCACAACTTGGCGGCTTCGAGGCTTTCAAGGAGG  >  ACB122/1084150‑1084410

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