Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,081,858 T→C 27.3% N104S (AAC→AGC)  ydhY ← predicted 4Fe‑4S ferridoxin‑type protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,081,8580TC27.3% 65.3 / 30.8 44N104S (AAC→AGC) ydhYpredicted 4Fe‑4S ferridoxin‑type protein
Reads supporting (aligned to +/- strand):  ref base T (13/19);  new base C (4/8);  total (17/27)
Fisher's exact test for biased strand distribution p-value = 7.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

TGCATTGACGGCAGGTGTCCGCGGTATAGTTGAGATCGCCATACAGGCCGCCGCCAGAGCCAACCCCGTTGTCGCCAAAGAAATAATTGCGATGGATTTTGATACGGGAGAAGAATGTCCCTACTGAGCCATCGTTGA  >  minE/1081790‑1081927
                                                                    |                                                                     
tGCATTGACGGCAGGTGTCCGCGGTATAGTTGAGATCGCCATACAGGCCGCCGCCAGAGCCAACCCCGTTg                                                                     >  1:1321332/1‑71 (MQ=255)
tGCATTGACGGCAGGTGTCCGCGGTATAGTTGAGATCGCCATACAGGCCGCCGCCAGAGCCAACCCCGCTg                                                                     >  1:49483/1‑71 (MQ=255)
 gCATTGACGGCAGGTGTCCGCGGTATAGTTGAGATCGCCATACAGGCCGCCGCCAGAGCCAACCCCGTTg                                                                     <  1:670073/70‑1 (MQ=255)
             ggTGTCCGCGGTATAGTTGAGATCGCCATACAGGCCGCCGCCAGAGCCAACCCCGTTGTCGCCAAAGAaat                                                        <  1:1137488/71‑1 (MQ=255)
             ggTGTCCGCGGTATAGTTGAGATCGCCATACAGGCCGCCGCCAGAGCCAACCCCGCTGTCGCCAAAGAaat                                                        <  1:1016005/71‑1 (MQ=255)
                   cgcgGTATAGTTGAGATCGCCATACAGGCCGCCGCCAGAGCCAACCCCGTTGTCGCCAAAGAAATAATTGc                                                  <  1:1611726/71‑1 (MQ=255)
                   cgcgGTATAGTTGAGATCGCCATACAGGCCGCCGCCAGAGCCAACCCCGCTGTCGCCAAAGAAATAATTGc                                                  <  1:911584/71‑1 (MQ=255)
                    gcgGTATAGTTGAGATCGCCATACAGGCCGCCGCCAGAGCCAACCCCGCTGTCGCCAAAGAaat                                                        <  1:809556/64‑1 (MQ=255)
                    gcgGTATAGTTGAGATCGCCATACAGGCCGCCGCCAGAGCCAACCCCGCTGTCGCCAAAGAaat                                                        <  1:617688/64‑1 (MQ=255)
                       gTATAGTTGAGATCGCCATACAGGCCGCCGCCAGAGCCAACCCCGCTGTCGCCa                                                               <  1:1052196/54‑1 (MQ=255)
                        tataGTTGAGATCGCCATACAGGCCGCCGCCAGAGCCAACCCCGTTGTCGCCAAAGAAATAATTGCGATgg                                             <  1:1410366/71‑1 (MQ=255)
                        tataGTTGAGATCGCCATACAGGCCGCCGCCAGAGCCAACCCCGTTGTCGCCAAAGAAATAATTGCGATg                                              >  1:1602271/1‑70 (MQ=255)
                        tataGTTGAGATCGCCATACAGGCCGCCGCCAGAGCCAACCCCGCTGTCGCCAAAGAaat                                                        <  1:792552/60‑1 (MQ=255)
                        tataGTTGAGATCGCCATACAGGCCGCCGCCAGAGCCAACCCCGCTGTCGCCAAAGAAATAATTGCGATg                                              <  1:1450989/70‑1 (MQ=255)
                          taGTTGAGATCGCCATACAGGCCGCCGCCAGAGCCAACCCCGCTGTCGCCAAAGAaataa                                                      >  1:903254/1‑60 (MQ=255)
                              tGAGATCGCCATACAGGCCGCCGCCAGAGCCAACCCCGTTGTCGCCAAAGAaataa                                                      >  1:1015650/1‑56 (MQ=255)
                               gagaTCGCCATACAGGCCGCCGCCAGAGCCAACCCCGTTGTCGCCAAAGaaa                                                         >  1:306690/1‑52 (MQ=255)
                               gagaTCGCCATACAGGCCGCCGCCAGAGCCAACCCCGTTGTCGCCAAAGAaata                                                       >  1:207000/1‑54 (MQ=255)
                               gagaTCGCCATACAGGCCGCCGCCAGAGCCAACCCCGTTGTCGCCAAAGAAATAATTGCGATGGATTTTGa                                      >  1:639075/1‑71 (MQ=255)
                               gagaTAGCCATACAGGCCGCCGCCAGAGCCAACCCCGTTGTCGCCAAAGAAATAATTGCGATGGATTTTGa                                      >  1:1222516/1‑71 (MQ=255)
                                agaTCGCCATACAGGCCGCCGCCAGAGCCAACCCCGTTGTCGCCAAAGAaataa                                                      <  1:1374717/54‑1 (MQ=255)
                                agaTCGCCATACAGGCCGCCGCCAGAGCCAACCCCGTTGTCGCCAAAGAaataa                                                      <  1:2084098/54‑1 (MQ=255)
                                  aTCGCCATACAGGCCGCCGCCAGAGCCAACCCCGTTGTCGCCAAAGAAATAATTGCGATGGATTTTGATAc                                   <  1:660571/71‑1 (MQ=255)
                                  aTCGCCATACAGGCCGCCGCCAGAGCCAACCCCGTTGTCGCCAAAGAAATAATTGCGATGGATTTTGATAc                                   <  1:2204201/71‑1 (MQ=255)
                                       cATACAGGCCGCCGCCAGAGCCAACCCCGTTGTCGCCAAAGAAATAATTGCGATGGATTTTGATACGGga                               <  1:140056/70‑1 (MQ=255)
                                        aTACAGGCCGCCGCCAGAGCCAACCCCGTTGTCGCCAAAGAAATAATTGCGATGGATTTTGa                                      <  1:1217686/62‑1 (MQ=255)
                                           cAGGCCGCCGCCAGAGCCAACCCCGTTGTCGCCAAAg                                                            >  1:1919031/1‑37 (MQ=255)
                                           cAGGCCGCCGCCAGAGCCAACCCCGTTGTCGCCAAAGAAATAATTGCGATGGATTTTGATACGGGagaag                           >  1:219430/1‑70 (MQ=255)
                                           cAGGCCGCCGCCAGAGCCAACCCCGTTGTCGCCAAAGAAATAATTGCGATGGATTTTGATACGGGagaag                           <  1:906491/70‑1 (MQ=255)
                                           cAGGCCGCCGCCAGAGCCAACCCCGTTGTCGCCAAAGAAATAATTGCGATGGATTTTGATACGGGagaag                           <  1:1112744/70‑1 (MQ=255)
                                               ccgccgccAGAGCCAACCCCGTTGTCGCCAAAGAAATAATTGCGATGGATTTTGATACGGGAGAAGAATGt                      <  1:561937/71‑1 (MQ=255)
                                                 tccgccAGAGCCAACCCCGTTGTCGCCAAAGAAATAATTGCGATGGATTTTGa                                      <  1:1569827/52‑1 (MQ=255)
                                                   cgccAGAGCCAACCGCGTTGTCGCCAAAGAAATAATTGCGATGGATTTTGATACGGGAGAAGAATGTccc                   >  1:408890/1‑70 (MQ=255)
                                                   cgccAGAGCCAACCCCGCTGTCGCCAAAGAAATAATTGCGATGGATTTTGATACGGGAGAAGAATGTccc                   <  1:1023141/70‑1 (MQ=255)
                                                         agCCAACCCCGTTGTCGCCAAAGAAATAATTGCGATGGATTTTGATACGGGAGAAGAATGTCCCTACTGAg            >  1:195838/1‑71 (MQ=255)
                                                          gCCAACCCCGTTGTCGCCAAAGAAATAATTGCGATGGATTTTGATACGGGAGAAGAATGTCCCTACTGAGc           >  1:1441972/1‑71 (MQ=255)
                                                          gCCAACCCCGTTGTCGCCAAAGAAATAATTGCGATGGATTTTGATACGGGAGAAGAATGTCCCTACTGAGc           >  1:1990404/1‑71 (MQ=255)
                                                          gCCAACCCCGTTGTCGCCAAAGAAATAATTGCGATGGATTTTGATACGGGAGAAGAATGTCCCTAATGAGc           <  1:1314607/71‑1 (MQ=255)
                                                               ccccGTTGTCGCCAAAGAAATAATTGCGATGGAtttt                                        <  1:2205584/37‑1 (MQ=255)
                                                               ccccGTTGTCGCCAAAGAAATAATTGCGATGGAtttt                                        <  1:709408/37‑1 (MQ=255)
                                                               ccccGTTGTCGCCAAAGAAATAATTGCGATGGATTTTGa                                      <  1:1904312/39‑1 (MQ=255)
                                                                cccGCTGTCGCCAAAGAAATAATTGCGATGGATTTTGATa                                    >  1:1563508/1‑40 (MQ=255)
                                                                 ccGCTGTCGCCAAAGAAATAATTGCGATGGATTTTGATACGGGAGAAGAATGTCCCTACTGAGCCAt        >  1:31428/1‑67 (MQ=255)
                                                                   gTTGTCGCCAAAGAAATAATTGCGATGGATTTTGATACGGGAGAAGAATGTCCCTACTGAGCCATCGTTGa  <  1:1100445/71‑1 (MQ=255)
                                                                    |                                                                     
TGCATTGACGGCAGGTGTCCGCGGTATAGTTGAGATCGCCATACAGGCCGCCGCCAGAGCCAACCCCGTTGTCGCCAAAGAAATAATTGCGATGGATTTTGATACGGGAGAAGAATGTCCCTACTGAGCCATCGTTGA  >  minE/1081790‑1081927

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