Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 99,507 T→C 39.2% T320A (ACC→GCC)  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,5070TC39.2% 38.1 / 57.1 84T320A (ACC→GCC) hofCassembly protein in type IV pilin biogenesis, transmembrane protein
Reads supporting (aligned to +/- strand):  ref base T (34/17);  new base C (19/14);  total (53/31)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

CCAGAGAGCCGGATGCCTCTCCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGAA  >  minE/99444‑99569
                                                               |                                                              
ccAGAGAGCCGGATGCCTCTCCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTAttt                                                            >  1:850760/1‑68 (MQ=255)
             tGCCTCTCCTGTTCTCACTAATTGCAAGCAGAGCGGGCTAAACTCACCGGTATTTTTTAGCgccagcca                                              <  1:2801608/69‑1 (MQ=255)
             tGCCTCTCCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCgccagcca                                              <  1:2147993/69‑1 (MQ=255)
             tGCCTCTCCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCgccagcca                                              <  1:1851431/69‑1 (MQ=255)
             tGCCTCTCCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCgccagcca                                              <  1:2144100/69‑1 (MQ=255)
             tGCCTCTCCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCgccagcca                                              <  1:1368005/69‑1 (MQ=255)
             tGCCTCTCCTGTGCTCACTAATTGCAAACAGAGCTGGCTAAACTCACCGGTATTTTTTAGCgccagcca                                              <  1:1634739/69‑1 (MQ=255)
              gCCTCTCCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGccagccag                                             >  1:226513/1‑69 (MQ=255)
              gCCTCTCCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGccagccag                                             >  1:2778613/1‑69 (MQ=255)
              gCCTCTCCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGccagccag                                             >  1:174849/1‑69 (MQ=255)
              gCCTCTCCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGccagccag                                             >  1:929528/1‑69 (MQ=255)
                  ctcCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCgccagcca                                              >  1:467390/1‑64 (MQ=255)
                  ctcCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCgccagcca                                              >  1:2564701/1‑64 (MQ=255)
                  ctcCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCgccagcca                                              >  1:574756/1‑64 (MQ=255)
                  ctcCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCgccagcca                                              >  1:2045418/1‑64 (MQ=255)
                  ctcCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCgccagcca                                              >  1:1932058/1‑64 (MQ=255)
                  ctcCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCgccagcca                                              >  1:1152115/1‑64 (MQ=255)
                  ctcCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCgccagcca                                              >  1:3412363/1‑64 (MQ=255)
                  ctcCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCgccagcca                                              >  1:2984836/1‑64 (MQ=255)
                  ctcCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATgg                                         <  1:2381645/69‑1 (MQ=255)
                  ctcCTGTTCTCACTAATCGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCgccagcca                                              >  1:2303308/1‑64 (MQ=255)
                  ctcATGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATgg                                         <  1:2960310/69‑1 (MQ=255)
                     cTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGtt                                      >  1:1998230/1‑69 (MQ=255)
                      tGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCgccagcca                                              <  1:1334743/60‑1 (MQ=255)
                          ctcACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTAt                                                              >  1:3300513/1‑40 (MQ=255)
                          ctcACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTAt                                                              >  1:1328664/1‑40 (MQ=255)
                          ctcACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTAt                                                              >  1:1534384/1‑40 (MQ=255)
                          ctcACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTAt                                                              >  1:1515827/1‑40 (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:2246421/68‑1 (MQ=255)
                            cACTAATTGCAAACAGAGCGGGCTAAACTCACCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGt                                <  1:1317550/68‑1 (MQ=255)
                            cACTAATTGCAAACAGAGCGGGCTAAACTCACCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGt                                <  1:3223972/68‑1 (MQ=255)
                            cACTAATTGCAAACAGAGCGGGCTAAACTCACCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGt                                <  1:2221395/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:1635903/69‑1 (MQ=255)
                                   tGCAAACAGAGCGGGCTAAACTCACCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGata                        <  1:454023/69‑1 (MQ=255)
                                   tGCAAACAGAGCGGGCTAAACTCACCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGata                        <  1:1916012/69‑1 (MQ=255)
                                   tGCAAACAGAGCGGGCTAAACTCACCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGata                        <  1:1201037/69‑1 (MQ=255)
                                    gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTa                                                        >  1:2087706/1‑36 (MQ=255)
                                    gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGatat                       >  1:806236/1‑69 (MQ=255)
                                    gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGatat                       >  1:2351385/1‑69 (MQ=255)
                                    gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGatat                       >  1:3322153/1‑69 (MQ=255)
                                    gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGatat                       >  1:2418175/1‑69 (MQ=255)
                                    gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGatat                       >  1:895724/1‑69 (MQ=255)
                                    gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGatat                       >  1:1959870/1‑69 (MQ=255)
                                    gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGatat                       >  1:876783/1‑69 (MQ=255)
                                    gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGatat                       >  1:1183991/1‑69 (MQ=255)
                                    gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGatat                       >  1:188407/1‑69 (MQ=255)
                                    gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGatat                       >  1:2981305/1‑69 (MQ=255)
                                    gCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGata                        >  1:3321966/1‑68 (MQ=255)
                                    gCAAACAGAGCGGGCTAAACTCACCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGata                        <  1:1706735/68‑1 (MQ=255)
                                      aaaCAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATgg                                         <  1:804753/49‑1 (MQ=255)
                                                ggCTAAACTCACCGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTtg          <  1:2791519/69‑1 (MQ=255)
                                                 gCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTtg          <  1:2244418/69‑1 (MQ=255)
                                                 gCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTtg          <  1:1486001/69‑1 (MQ=255)
                                                 gCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTtg          <  1:767939/69‑1 (MQ=255)
                                                 gCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTtg          <  1:196537/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:1094611/1‑68 (MQ=255)
                                                          aCCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa  >  1:970722/1‑68 (MQ=255)
                                                          aCCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa  >  1:3075169/1‑68 (MQ=255)
                                                          aCCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa  >  1:2790698/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: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:1153097/1‑68 (MQ=255)
                                                          aCCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa  >  1:1104205/1‑68 (MQ=255)
                                                          aCCGGCATTTTTTAGCGCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGa   >  1:3187159/1‑67 (MQ=255)
                                                          aCCGGCATTTTTTAGCCCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa  >  1:930179/1‑68 (MQ=255)
                                                          aCCGGCATTTTTTAGCCCCAGCCAGATGGGGTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGaa  >  1:934087/1‑68 (MQ=255)
                                                               |                                                              
CCAGAGAGCCGGATGCCTCTCCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGAA  >  minE/99444‑99569

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: