Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 99,532 T→G 51.7% Q311H (CAA→CAC hofC ← assembly protein in type IV pilin biogenesis, transmembrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE99,5320TG51.7% 13.2 / 62.9 62Q311H (CAA→CAChofCassembly protein in type IV pilin biogenesis, transmembrane protein
Reads supporting (aligned to +/- strand):  ref base T (16/14);  new base G (18/14);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGAAGTTGCACCCAGTACGGGCAGCG  >  minE/99465‑99591
                                                                   |                                                           
cTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGtt                                                            >  1:1998230/1‑69 (MQ=255)
     ctcACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGAcc                                                        >  1:1864966/1‑68 (MQ=255)
       cACTAATTTCAAACAGAGCGGGCTAAACTCACCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGt                                                      <  1:373994/68‑1 (MQ=255)
       cACTAATTGCAAACAGAGCGGGCTAAACTCACCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGt                                                      <  1:1959136/68‑1 (MQ=255)
       cACTAATTGCAAACAGAGCGGGCTAAACTCACCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGt                                                      <  1:873958/68‑1 (MQ=255)
       cACTAATTGCAAACAGAGCGGGCTAAACTCACCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGt                                                      <  1:3223972/68‑1 (MQ=255)
       cACTAATTGCAAACAGAGCGGGCTAAACTCACCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGt                                                      <  1:2246421/68‑1 (MQ=255)
       cACTAATTGCAAACAGAGCGGGCTAAACTCACCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGt                                                      <  1:2221395/68‑1 (MQ=255)
       cACTAATTGCAAACAGAGCGGGCTAAACTCACCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGt                                                      <  1:1317550/68‑1 (MQ=255)
       cACTAATTGCAAACAGAGCGGGCTAAACTCACCGGCATTTTTTAGCGCCAGCCAGATGGGGCGACCGt                                                      <  1:4119/68‑1 (MQ=255)
        aCTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTa                                                    >  1:3176997/1‑69 (MQ=255)
              tGCAAACAGAGCGGGCTAAACTCACCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGata                                              <  1:377610/69‑1 (MQ=255)
              tGCAAACAGAGCGGGCTAAACTCACCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGata                                              <  1:1916012/69‑1 (MQ=255)
              tGCAAACAGAGCGGGCTAAACTCACCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGata                                              <  1:1635903/69‑1 (MQ=255)
              tGCAAACAGAGCGGGCTAAACTCACCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGata                                              <  1:454023/69‑1 (MQ=255)
              tGCAAACAGAGCGGGCTAAACTCACCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGata                                              <  1:1201037/69‑1 (MQ=255)
               gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGatat                                             >  1:3322153/1‑69 (MQ=255)
               gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGatat                                             >  1:2351385/1‑69 (MQ=255)
               gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGatat                                             >  1:2418175/1‑69 (MQ=255)
               gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGatat                                             >  1:2981305/1‑69 (MQ=255)
               gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGatat                                             >  1:895724/1‑69 (MQ=255)
               gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGatat                                             >  1:1959870/1‑69 (MQ=255)
               gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGatat                                             >  1:188407/1‑69 (MQ=255)
               gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGatat                                             >  1:876783/1‑69 (MQ=255)
               gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGatat                                             >  1:806236/1‑69 (MQ=255)
               gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGatat                                             >  1:1183991/1‑69 (MQ=255)
               gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGata                                              >  1:3321966/1‑68 (MQ=255)
               gCAAACAGAGCGGGCTAAACTCACCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGata                                              <  1:1706735/68‑1 (MQ=255)
                           ggCTAAACTCACCGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTtg                                <  1:2791519/69‑1 (MQ=255)
                            gCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTtg                                <  1:2244418/69‑1 (MQ=255)
                            gCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTtg                                <  1:767939/69‑1 (MQ=255)
                            gCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTtg                                <  1:196537/69‑1 (MQ=255)
                            gCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTtg                                <  1:1486001/69‑1 (MQ=255)
                            gCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTtg                                <  1:1982867/69‑1 (MQ=255)
                               aaaCTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTTgtgt                             <  1:2262858/69‑1 (MQ=255)
                                    cACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa                        >  1:1265049/1‑69 (MQ=255)
                                     aCCGGCATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa                        >  1:2149089/1‑68 (MQ=255)
                                     aCCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa                        >  1:497461/1‑68 (MQ=255)
                                     aCCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa                        >  1:3075169/1‑68 (MQ=255)
                                     aCCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa                        >  1:970722/1‑68 (MQ=255)
                                     aCCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa                        >  1:2631769/1‑68 (MQ=255)
                                     aCCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa                        >  1:2607989/1‑68 (MQ=255)
                                     aCCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa                        >  1:2464275/1‑68 (MQ=255)
                                     aCCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa                        >  1:2790698/1‑68 (MQ=255)
                                     aCCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa                        >  1:2348154/1‑68 (MQ=255)
                                     aCCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa                        >  1:211963/1‑68 (MQ=255)
                                     aCCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa                        >  1:2057736/1‑68 (MQ=255)
                                     aCCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa                        >  1:200294/1‑68 (MQ=255)
                                     aCCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa                        >  1:1988884/1‑68 (MQ=255)
                                     aCCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa                        >  1:1094611/1‑68 (MQ=255)
                                     aCCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa                        >  1:1153097/1‑68 (MQ=255)
                                     aCCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa                        >  1:1104205/1‑68 (MQ=255)
                                     aCCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGa                         >  1:3187159/1‑67 (MQ=255)
                                     aCCGGCATTTTTTAGCCCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa                        >  1:934087/1‑68 (MQ=255)
                                     aCCGGCATTTTTTAGCCCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa                        >  1:930179/1‑68 (MQ=255)
                                            ttttttAGCGCCAGCCAGATGGGTTGACCGTTACTGat                                               >  1:1013949/1‑38 (MQ=255)
                                                ttAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGAAGTTGCACCCAg             <  1:1643750/68‑1 (MQ=255)
                                                          ccagATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGAAGTTGCACCCAGTACGGGCAgcg  <  1:792565/69‑1 (MQ=255)
                                                          ccagATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGAAGTTGCACCCAGTACGGGCAgcg  <  1:1166313/69‑1 (MQ=255)
                                                          ccagATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGAAGTTGCACCCAGTACGGGCAgcg  <  1:1072335/69‑1 (MQ=255)
                                                          ccagATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGAAGTTGCACCCAGTACGGGCAgcg  <  1:102737/69‑1 (MQ=255)
                                                          ccagATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGAAGTTGCACCCAGTACGGGCAgcg  <  1:1529144/69‑1 (MQ=255)
                                                          ccagATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGAAGTTGCACCCAGTACGGGCAgcg  <  1:2401596/69‑1 (MQ=255)
                                                                   |                                                           
CTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGAAGTTGCACCCAGTACGGGCAGCG  >  minE/99465‑99591

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: