Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,699,330 G→A 100% K69K (AAG→AAA yrfC → predicted fimbrial assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,699,3300GA78.7% 34.2 / 4.2 28K69K (AAG→AAAyrfCpredicted fimbrial assembly protein
Reads supporting (aligned to +/- strand):  ref base G (1/5);  new base A (0/22);  total (1/27)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAGCCACGTTTGCTGGAGCAACAGCAGTTACGCGAACAGCGTTCGCAACGGCAGCGCCAGCGACAA  >  minE/2699267‑2699393
                                                               |                                                               
acgcgTAAACGCCGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAACCACg                                                            <  1:52725/68‑1 (MQ=255)
acgcgTAAACGCCGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAACCACg                                                            <  1:401013/68‑1 (MQ=255)
acgcgTAAACGCCGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAACCACg                                                            <  1:1363984/68‑1 (MQ=255)
acgcgTAAACGCCGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAACCACg                                                            <  1:150539/68‑1 (MQ=255)
acgcgTAAACGCCGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAACCACg                                                            <  1:2974977/68‑1 (MQ=255)
acgcgTAAACGCCGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAACCACg                                                            <  1:1888591/68‑1 (MQ=255)
acgcgTAAACGCCGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAACCACg                                                            <  1:1996006/68‑1 (MQ=255)
acgcgTAAACGCCGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAACCACg                                                            <  1:2006716/68‑1 (MQ=255)
acgcgTAAACGCCGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAACCACg                                                            <  1:2411036/68‑1 (MQ=255)
  gcATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAGCCACGtt                                                          <  1:1161276/69‑1 (MQ=255)
  gcATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAGCCACGtt                                                          <  1:3083764/69‑1 (MQ=255)
         cGCCGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAAGGAAACCACGTTTGCTgg                                                    <  1:2924982/68‑1 (MQ=255)
         cGCCGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAACCACGTTTGCTgg                                                    <  1:2466982/68‑1 (MQ=255)
         cGCCGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAACCACGTTTGCTgg                                                    <  1:16240/68‑1 (MQ=255)
         cGCCGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAACCACGTTTGCTgg                                                    <  1:1348465/68‑1 (MQ=255)
         cGCCGTTTTGCTTAAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAACCACGTTTGCTgg                                                    <  1:3073090/68‑1 (MQ=255)
         cGCCGTTTTGCTTAAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAACCACGTTTGCTgg                                                    <  1:3360324/68‑1 (MQ=255)
         cGCCGTTTTGCTTAAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAACCACGTTTGCTgg                                                    <  1:3368341/68‑1 (MQ=255)
         cGCCGTTTTGCTTAAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAACCACGTTTGCTgg                                                    <  1:1219149/68‑1 (MQ=255)
         cGCCGTTTTGCTTAAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAACCACGTTTGCTgg                                                    <  1:677923/68‑1 (MQ=255)
         cGCCGTTTTGCTTAAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAACCACGTTTGCTgg                                                    <  1:730466/68‑1 (MQ=255)
           ccgTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAACCACGTTTGCTgg                                                    <  1:3128123/64‑1 (MQ=255)
               tttGCTTAAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAACCACGTTTGCTgg                                                    <  1:1064212/62‑1 (MQ=255)
                ttGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAGCCACGTTTGCTGGAGCAAc                                              >  1:3273990/1‑67 (MQ=255)
                     tCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAGCCACGTTTGCTGGAGCAACAGCAGt                                        <  1:2439133/68‑1 (MQ=255)
                     tCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAGCCACGTTTGCTGGAGCAACAGCAGt                                        <  1:147608/68‑1 (MQ=255)
                           ggAACAACAACTCGCCCGCAGCTTACAGATAACGAAACCACGTTTGCTgg                                                    <  1:334380/50‑1 (MQ=255)
                                                          aCGAAGCCACGTTTGCTGGAGCAACAGCAGTTACGCGAACATCGTTCGCAACGGCAGCGCCAGCGACaa  <  1:2295984/69‑1 (MQ=255)
                                                               |                                                               
GCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAGCCACGTTTGCTGGAGCAACAGCAGTTACGCGAACAGCGTTCGCAACGGCAGCGCCAGCGACAA  >  minE/2699267‑2699393

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