Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,708,543 C→T 7.7% P222P (CCG→CCA)  cysI ← sulfite reductase hemoprotein subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,708,5430CT 7.7% 150.6 / 2.2 52P222P (CCG→CCA) cysIsulfite reductase hemoprotein subunit beta
Reads supporting (aligned to +/- strand):  ref base C (26/22);  new base T (2/2);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.12e-01

CCATTTCGGCCTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTCCTGC  >  NC_002947/2708395‑2708690
                                                                                                                                                     |                                                                                                                                                   
ccATTTCGGCCTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCAc                                                                                                                                                     >  4:262370/1‑149 (MQ=255)
         ccTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAAc                                                                                                                                            >  6:217030/1‑149 (MQ=255)
          cTCGACCTTTTCGGCGAACTCTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCGTGGCAGGACCt                                                                                                                                           <  4:257867/149‑1 (MQ=255)
             gACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTc                                                                                                                                          <  3:262370/147‑1 (MQ=255)
              aCCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGtt                                                                                                                                       <  1:81416/149‑1 (MQ=255)
               ccctttCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTtg                                                                                                                                      <  8:212466/146‑1 (MQ=255)
                    tCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCGGGTCCTGCCACGGCAGGAACTCGTtgatga                                                                                                                                  <  3:366753/148‑1 (MQ=255)
                      gtcGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAAc                                                                                                                                <  4:234573/146‑1 (MQ=255)
                                 ccGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCccacca                                                                                                                      <  4:250136/147‑1 (MQ=255)
                                  cGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACTGCAGGAACTCGTTCATCAACGAgttc                                                                                                                          <  5:146586/142‑4 (MQ=255)
                                  cGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCg                                                                                                                   >  4:317904/1‑149 (MQ=255)
                                  cGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTCATCAACGAgttc                                                                                                                          >  6:146586/1‑139 (MQ=255)
                                      gTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCtggg                                                                                                                >  4:20383/1‑148 (MQ=255)
                                            gCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGcc                                                                                                         <  3:317904/149‑1 (MQ=255)
                                               ttGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAg                                                                                                      >  3:160676/1‑149 (MQ=255)
                                                    cAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCGTGGCAGGacc                                                                                                                                            <  2:42262/106‑3 (MQ=255)
                                                    cAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCGTGGCAGGacc                                                                                                                                            >  1:42262/1‑104 (MQ=255)
                                                         tCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCaccacc                                                                                             >  5:365543/1‑148 (MQ=255)
                                                               tCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGc                                                                                      <  3:181253/149‑1 (MQ=255)
                                                                  gcgcCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGcaca                                                                                   >  8:23839/1‑149 (MQ=255)
                                                                               gttgtCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCcca                                                                                                                         >  6:25684/1‑98 (MQ=255)
                                                                               gttgtCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCcca                                                                                                                         <  5:25684/98‑1 (MQ=255)
                                                                                   tCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGcc                                                                  >  5:328724/1‑149 (MQ=255)
                                                                                   tCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGcc                                                                  >  3:158357/1‑149 (MQ=255)
                                                                                   tCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGcc                                                                  >  3:180332/1‑149 (MQ=255)
                                                                                     aCGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCgg                                                                >  3:175833/1‑149 (MQ=255)
                                                                                          gACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGCACGAGCCCACCACCGGGGTACGGCCCAGGCAACCACATACCAGCACACGGAAGCCCAGTTCGCCGGCAGCg                                                           >  7:342552/1‑149 (MQ=255)
                                                                                            ccGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCgct                                                         <  4:160676/149‑3 (MQ=255)
                                                                                             cGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTg                                                        >  6:160590/1‑149 (MQ=255)
                                                                                                 aCGGTTGTACACGCGCAGGATGGCTGTCGAGGTAGCTGATCAGGACCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTgcg                                                      >  7:234063/1‑148 (MQ=255)
                                                                                                         acacGCGGAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTATGGCCCAGGCCACCACCTACCAGCACACGGaa                                                                              <  5:320372/115‑1 (MQ=255)
                                                                                                         acacGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTATGGCCCAGGCCACCACCTACCAGCACACGGaa                                                                              >  6:320372/1‑115 (MQ=255)
                                                                                                         acacGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTc                                             <  3:20383/148‑1 (MQ=255)
                                                                                                          cacGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCa                                                                                                       >  3:192758/1‑89 (MQ=255)
                                                                                                          cacGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCa                                                                                                       <  4:192758/89‑1 (MQ=255)
                                                                                                           acGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAg                                          >  5:20861/1‑149 (MQ=255)
                                                                                                             gcgcAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGc                                        <  4:180332/149‑1 (MQ=255)
                                                                                                                 aGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAAt                                    >  3:29021/1‑149 (MQ=255)
                                                                                                                   gATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGt                                  >  1:13945/1‑149 (MQ=255)
                                                                                                                    aTGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTc                                 <  6:53956/149‑1 (MQ=255)
                                                                                                                    aTGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTc                                 <  4:175833/149‑1 (MQ=255)
                                                                                                                                   gCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATg                    <  6:363791/148‑1 (MQ=255)
                                                                                                                                    cTGATCAGGACCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATggc                  <  8:234063/149‑1 (MQ=255)
                                                                                                                                     tGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCAGCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGGAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATggcg                 >  1:265108/1‑149 (MQ=255)
                                                                                                                                     tGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATggcg                 >  2:174799/1‑149 (MQ=255)
                                                                                                                                       acgagttCCTGCCGTGACAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACGTCGATggcggc               >  4:125093/7‑149 (MQ=255)
                                                                                                                                        tCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTCCGGCCCAGGCCTCCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGgcggcg              >  5:325475/1‑149 (MQ=255)
                                                                                                                                          aGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGgcgc             <  5:217030/148‑1 (MQ=255)
                                                                                                                                           ggTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGTGgtgctg           <  4:29021/149‑3 (MQ=255)
                                                                                                                                                tGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATgg                   <  4:338261/136‑1 (MQ=255)
                                                                                                                                                tGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATgg                   >  3:338261/1‑136 (MQ=255)
                                                                                                                                                  ccACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACAGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTCCt    >  7:126270/1‑149 (MQ=255)
                                                                                                                                                    aCGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACGCGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTAGTGCACTTCGATGGCGGCGCGGGCTTCCTGc  >  8:140956/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCATTTCGGCCTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGC‑GTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTCCTGC  >  NC_002947/2708395‑2708690

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