Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 5,856,018 C→A 7.7% M202I (ATG→ATT)  cysA ← Sulfate/thiosulfate import ATP‑binding proteinCysA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,856,0180CA 7.7% 157.4 / 4.5 52M202I (ATG→ATT) cysASulfate/thiosulfate import ATP‑binding proteinCysA
Reads supporting (aligned to +/- strand):  ref base C (27/21);  new base A (2/2);  total (29/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 = 4.86e-01

CAGCCAGTTGCCCTGGCCGACGAACTCGGCCACGAACGGCGTGGCCGGGTGGTCATAGATTTCCTGCGGGGTGGCGTACTGCTCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCAGTTCGGTGCACAAGTGCTCGCGCACCCGGGCGTCCAGCGCCGACATCGGCTCATCGAGCAGCAGCA  >  KT2440/5855892‑5856151
                                                                                                                              |                                                                                                                                     
cAGCCAGTTGCCCTGGCCGACGAACTCGGCCACGAACGGCGTGGCCGGGTGGTCATAGATTTCCTGCGGGGTGGCGTACTGCTCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTc                                                                                                                               <  1:251658/135‑1 (MQ=255)
      gTTGCCCTGGCCGACGAACTCTGCCACGAACGGCGTGGCCGGGTGGTCATAGATTTCCTGCGGGGTGGCGTACTGCTCGAGCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTgg                                                                                                                        <  8:229713/136‑1 (MQ=255)
           ccTGGCCGACGAACTCGGCCACGAACGGCGTGGCCGGGTGGTCATAGATTTCCTGCGGGGTGGCGTACTGCTCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTtgtg                                                                                                                   <  7:203005/136‑1 (MQ=255)
            cTGGCCGACGAACTCGGCCACGAACGGCGTGGCCGGGTGGTCATAGATTTCCTGCGGGGTGGCGTACTGCTCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTtgtg                                                                                                                   <  4:108918/135‑1 (MQ=255)
                ccGACGAACTCGGCCACGAACGGCGTGGCCGGGTGGTCATAGATTTCCTGCGGGGTGGCGTACTGCTCGACCTGGCCGGTGGTCATCACCGCGATGCGGTCGGCCATCCACATGGCCTCGTCCTGGTTGGGGGtt                                                                                                               >  2:107033/1‑134 (MQ=255)
                  gacgaACTCGGCCACGAACGGCGTGGCCGGGTGGTCATAGATTTCCTGCGGGGTGGCGTACTGCTCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAAc                                                                                                             >  8:232249/1‑135 (MQ=255)
                       aCTCGGCCACGAACGGCGTGGCCGGGTGGTCATAGATTTCCTGCGGGGTGGCGTACTGCTCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCa                                                                                                        >  1:177639/1‑135 (MQ=255)
                       aCTCGGCCACGAACGGCGTGGCCGGGTGGTCATAGATTTCCTGCGGGGTGGCGTACTGCTCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCa                                                                                                        >  3:79245/1‑135 (MQ=255)
                                     ggCGTGGCCGGGTGGTCATAGATTTCCTGCGGGGTGGCGTACTGCTCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGccc                                                                                          <  3:193960/135‑1 (MQ=255)
                                        gTGGCCGGGTGGTCATAGATTTCCTGCGGGGTGGCGTACTGCTCGCCCTGGCCGTTGTTGATCACCGCGATGCGGTCGGCCATCACCCTTGCCGCGGCCGGGTGGTGGGGTACCATCCACGGGGGGATGGCCAGt                                                                                       >  6:157911/1‑135 (MQ=255)
                                          ggCCGGGTGGTCATAGATTTCCTGCGGGGTGGCGTACTGCTCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTgg                                                                                    >  7:88311/1‑136 (MQ=255)
                                               ggtggtCATAGATTTCCTGCGGGGTGGCGTTCTGCTCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTg                                                                               <  7:112512/136‑1 (MQ=255)
                                               ggtggtCATAGATTTCCTGCGGGGTGGCGTACTGCTCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATAAACATGGCCTCGTCCTGGTTGTTGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCt                                                                                >  2:78469/1‑135 (MQ=255)
                                                  ggtCATAGATTTCCTGCGGGGTGGCGTACTGCACGACCTGGCCCTTGTTCATAACCGCGATGGGGTCGGCCCACAACAAGGCCTCGGCGTGGGTGTGGGTTACCCAAAAAGTTGTGAAGCCCCCTTGgcccagcc                                                                             >  2:259290/1‑129 (MQ=255)
                                                           tttCCTGCGGGGTGGCGTACTGCTCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCa                                                                    <  5:15623/135‑1 (MQ=255)
                                                           tttCCTGCGGGGTGGCGTACTGCTCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCa                                                                    <  4:267803/135‑1 (MQ=255)
                                                           tttCCTGCGGGGTGGCGTACTGCTCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCa                                                                    <  2:177639/135‑1 (MQ=255)
                                                              ccTGCGGGGTGGCGTACTGCTCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCAGt                                                                  <  2:88108/134‑1 (MQ=255)
                                                                tGCGGGGTGGCGTACTGCTCGACCTGGACGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGCCGCAGTTCg                                                               >  4:168083/1‑135 (MQ=255)
                                                                 gCGGGGTGGCGTACTGCTCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTtgtg                                                                                                                   <  1:165698/82‑1 (MQ=255)
                                                                 gCGGGGTGGCGTACTGCTCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTtgtg                                                                                                                   >  2:165698/1‑82 (MQ=255)
                                                                     ggTGGCGTACTGCTCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCAGTTCGGTGc                                                           >  3:245525/1‑134 (MQ=255)
                                                                               tGCTCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCa                                                                    >  2:227198/1‑115 (MQ=255)
                                                                               tGCTCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCa                                                                    <  1:227198/115‑1 (MQ=255)
                                                                                 cTCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCAGTTCGGTGCACAAGTGCTcgc                                               <  7:27074/134‑1 (MQ=255)
                                                                                  tCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCTGCCATCAACATGGCCTCGGCCTGGTTGTGGGTAACCATCAAAGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCAGTTCGGTGCACAAGTGCTcgc                                               <  4:174059/133‑1 (MQ=255)
                                                                                    gACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCAGTTCGGTGCACAAGTGCTCGCGCAc                                           >  3:159749/1‑135 (MQ=255)
                                                                                    gACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGCAACCATCa                                                                                                        <  8:175958/74‑1 (MQ=255)
                                                                                    gACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGCAACCATCa                                                                                                        >  7:175958/1‑74 (MQ=255)
                                                                                      ccTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCAGTTCGGTGCACAAGTGCTCGCGCAccc                                         >  3:192787/1‑135 (MQ=255)
                                                                                       cTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCAGTTCGGTGCACAAGTGCTCGCGCACCCgg                                       >  8:180182/1‑136 (MQ=255)
                                                                                         ggCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCAGTTCGGTGCACAAGTGCTCGCGCACCCggg                                      >  4:136542/1‑135 (MQ=255)
                                                                                            cGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCAGTTCGGTGCACAAGTGCTCGCGCACCCGGGCGTc                                  >  7:221680/1‑136 (MQ=255)
                                                                                            cGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCAGGTCGGTGCACAAGTGCTCGCGCACCCGGGCGt                                   >  3:156052/1‑135 (MQ=255)
                                                                                            cGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCGCGTCCTGGTTGTGGGTTACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCAGGCGCAGTTCGGTGCACACGTGCTCGCGCCCCCGGGCGt                                   >  3:190562/1‑135 (MQ=255)
                                                                                            cGTTGTTCATCACCGCGATGCGGTCGGCCATCAAAATgg                                                                                                                                   >  8:197531/1‑39 (MQ=255)
                                                                                            cGTTGTTCATCACCGCGATGCGGTCGGCCATCAAAATgg                                                                                                                                   >  8:193998/1‑39 (MQ=255)
                                                                                            cGTTGTTCATCACCGCGATGCGGTCGGCCATCAAAATgg                                                                                                                                   <  7:193998/39‑1 (MQ=255)
                                                                                            cGTTGTTCATCACCGCGATGCGGTCGGCCATCAAAATgg                                                                                                                                   <  7:197531/39‑1 (MQ=255)
                                                                                            cGTTGTTCATCACCCCGATGCGGTCGGCCATCAACATGGCCACGTCCATGTTGTGGGTAACAAACAACGGGGTGCAGCCCACGTGGCGCTGCCGCTGACGCAGTTCGGTGCACAACTGCTCGCGCACCCGGGCGt                                   >  5:30424/1‑135 (MQ=255)
                                                                                             gttgttCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCAGTTCGGTGCACAAGTGCTCGCGCACCCGGGCGTc                                  <  1:259290/135‑1 (MQ=255)
                                                                                               tgttCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCAGTTCGGTGCACAAGTGCTCGCGCACCCGGGCGTCCa                                <  4:100803/135‑1 (MQ=255)
                                                                                                   catcaCCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGCTGCCCAGTTGGCGCTGCAGCTGCCGCAGTTCGGTGCACCACTGCTCGCGCACCCGGGCGTCCAccgc                            >  6:71114/1‑135 (MQ=255)
                                                                                                           cgATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCAGTTCGGTGCACAAGTGCTCGCGCACCCGGGCGTCCAGCGCCGACATCg                    <  4:79245/135‑1 (MQ=255)
                                                                                                            gATGCGGTCGGCCATCAACATGGCCTCGTCCTGGGTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCAGTTCGGTGCACAAGTGCTCGCGCACCCGGGCGTCCAGCGCCGACATCGGc                  <  8:192941/136‑1 (MQ=255)
                                                                                                              tGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGCTGCCCAGTTGGCGCTGCCGCTGACGCAGTTCGGTGCACAAGTGCTCGCGCACCCGGGCGTCCCGCGCCGCCATCGGcg                 >  5:247396/1‑134 (MQ=255)
                                                                                                              tGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCAGTTCGGTGCACAAGTGCTCGCGCACCCGGGCGTCCAGCGCCGACATCGGc                  >  6:89590/1‑134 (MQ=255)
                                                                                                              tGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCAGTTCGGTGCACAAGTGCTCGCGCACCCGGGCGTCCAGCGCCGACATCGGCt                 >  7:15480/1‑135 (MQ=255)
                                                                                                              tGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCAGTTCGGTGCACAAGTGCTCGCGCACCCGGGCGTCCAGCGCCGACATCGGCt                 <  5:71187/135‑1 (MQ=255)
                                                                                                                 ggTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCAGTTCGGTGCACAAGTGCTCGCGCACCCGGGCGTCCAGCGCCGACATCGGCTCATc             <  8:222134/136‑1 (MQ=255)
                                                                                                                     ggCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCAGTTCGGTGCACAAGTGCTCGCGCACCCGGGCGTCCAGCGCCGACCTCGGCTCATCGa           >  4:231268/1‑134 (MQ=255)
                                                                                                                             aCATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCAGTTCGGTGCACAAGTGCTCGCGCACCCGGGCGTCCAGCGCCGCCATCGGCTCATAGagcagcagca  >  3:61590/1‑135 (MQ=255)
                                                                                                                              |                                                                                                                                     
CAGCCAGTTGCCCTGGCCGACGAACTCGGCCACGAACGGCGTGGCCGGGTGGTCATAGATTTCCTGCGGGGTGGCGTACTGCTCGACCTGGCCGTTGTTCATCACCGCGATGCGGTCGGCCATCAACATGGCCTCGTCCTGGTTGTGGGTAACCATCAACGTGGTGATGCCCAGTTGGCGCTGCAGCTGACGCAGTTCGGTGCACAAGTGCTCGCGCACCCGGGCGTCCAGCGCCGACATCGGCTCATCGAGCAGCAGCA  >  KT2440/5855892‑5856151

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