Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 912,872 T→A 7.7% A455A (GCA→GCT hcp ← hybrid‑cluster [4Fe‑2S‑2O] protein in anaerobic terminal reductases

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb912,8720TA 7.7% 138.1 / 4.0 52A455A (GCA→GCThcphybrid‑cluster [4Fe‑2S‑2O] protein in anaerobic terminal reductases
Reads supporting (aligned to +/- strand):  ref base T (30/18);  new base A (2/2);  total (32/20)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

CAGCCCAGTTTCTCTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTAC  >  W3110S.gb/912807‑912936
                                                                 |                                                                
cAGCCCAGTTTCTCTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATc                                                             <  1:700022/71‑1 (MQ=255)
   cccAGTTTCTCTGCCAGCGTTACCGCCAGAATAATTGCTGAGTAAGCATCGTTACACTGACCAGCATCTAc                                                          <  1:3818234/71‑1 (MQ=255)
   cccAGTTTCTCTGCCAGCGTTACCGCCAGAATAATTGCTGAGTAAGCATCGTTACACTGACCAGCATCTAc                                                          <  1:1907604/71‑1 (MQ=255)
   cccAGTTTCTCTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTAc                                                          <  1:642188/71‑1 (MQ=255)
    ccAGTTTCTCTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTAc                                                          <  1:1073011/70‑1 (MQ=255)
     cAGTTTCTCTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTAc                                                          >  1:2965863/1‑69 (MQ=255)
      aGTTTCTCTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTAc                                                          >  1:2490347/1‑68 (MQ=255)
         ttCTCTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGAcg                                                    >  1:660504/1‑71 (MQ=255)
          tctcTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGAcgc                                                   >  1:1441331/1‑71 (MQ=255)
          tctcTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGAcgc                                                   >  1:1810274/1‑71 (MQ=255)
           ctctGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGAcgcg                                                  >  1:3484196/1‑71 (MQ=255)
           ctctGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGAcgcg                                                  >  1:2581759/1‑71 (MQ=255)
              tGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCgg                                                 >  1:2364219/1‑69 (MQ=255)
                ccAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAg                                              <  1:795755/70‑1 (MQ=255)
                 cAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGAc                                            >  1:1896070/1‑71 (MQ=255)
                 cAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGAc                                            >  1:1346845/1‑71 (MQ=255)
                      tGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTc                                        >  1:142846/1‑70 (MQ=255)
                      tGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTc                                        >  1:3546639/1‑70 (MQ=255)
                      tGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTc                                        <  1:2787180/70‑1 (MQ=255)
                       gACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGCAGACCTTCGa                                      >  1:554944/1‑70 (MQ=255)
                        aCAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGa                                      <  1:924185/70‑1 (MQ=255)
                                 ataatCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCa                                               >  1:1046936/1‑52 (MQ=255)
                                 ataatCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCa                                               >  1:538957/1‑52 (MQ=255)
                                 ataatCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCaaa                            >  1:605747/1‑71 (MQ=255)
                                  taatCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAAc                           >  1:800709/1‑71 (MQ=255)
                                    atCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTc                                        <  1:401006/56‑1 (MQ=255)
                                    atCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACt                          >  1:470968/1‑70 (MQ=255)
                                    atCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACt                          >  1:2948252/1‑70 (MQ=255)
                                    atCGCTGAGTAAGCATCGTTACATTGACATGAATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACt                          >  1:3138721/1‑70 (MQ=255)
                                       gCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAg                      >  1:2474305/1‑71 (MQ=255)
                                       gCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAg                      >  1:2346518/1‑71 (MQ=255)
                                          gAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGttt                   >  1:975018/1‑71 (MQ=255)
                                           aGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGttt                   <  1:3575968/70‑1 (MQ=255)
                                           aGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGttt                   <  1:2599789/70‑1 (MQ=255)
                                            gTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAg                      <  1:2612176/66‑1 (MQ=255)
                                            gTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAg                      <  1:2581119/66‑1 (MQ=255)
                                             tAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGtt                <  1:1572786/71‑1 (MQ=255)
                                              aaGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAATTTTGtt                <  1:495843/70‑1 (MQ=255)
                                                  aTCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGc           >  1:1578236/1‑71 (MQ=255)
                                                   tCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGa                                      >  1:2905760/1‑43 (MQ=255)
                                                   tCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCa                        <  1:1416761/57‑1 (MQ=255)
                                                    cGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGa         >  1:3544907/1‑71 (MQ=255)
                                                     gTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCaaa                            >  1:3083705/1‑51 (MQ=255)
                                                     gTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCaaa                            >  1:2851406/1‑51 (MQ=255)
                                                     gTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCaaa                            >  1:2284924/1‑51 (MQ=255)
                                                     gTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGat        <  1:2502936/71‑1 (MQ=255)
                                                        aCATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGatat      <  1:910689/70‑1 (MQ=255)
                                                        aCATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATAtt     >  1:1043422/1‑71 (MQ=255)
                                                         cATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATAtt     <  1:3291917/70‑1 (MQ=255)
                                                             gACCAGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGttt                   >  1:1836102/1‑52 (MQ=255)
                                                             gACCAGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGttt                   >  1:2313983/1‑52 (MQ=255)
                                                                cTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATAtt     <  1:2820737/63‑1 (MQ=255)
                                                                 tGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTAc  <  1:3213124/65‑1 (MQ=255)
                                                                 |                                                                
CAGCCCAGTTTCTCTGCCAGAGTGACAGCCAGAATAATCGCTGAGTAAGCATCGTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTAC  >  W3110S.gb/912807‑912936

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