Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,699,171 C→T 6.3% T16T (ACC→ACT yrfC → predicted fimbrial assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,699,1710CT 6.3% 195.0 / 2.6 63T16T (ACC→ACTyrfCpredicted fimbrial assembly protein
Reads supporting (aligned to +/- strand):  ref base C (29/30);  new base T (2/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

CGCTTGGGAAGGCGTACTGATGAACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGTCGGGATAACGCTAATA  >  minE/2699105‑2699240
                                                                  |                                                                     
cgcTTGGGAAGGCGTACTGATGAACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCtt                                                                   >  1:2176815/1‑71 (MQ=255)
 gcTTGGGAAGGCGTACTGATGAACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCttt                                                                  >  1:3205063/1‑71 (MQ=255)
 gcTTGGGAAGGCGTACTGATGAACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCtt                                                                   <  1:294072/70‑1 (MQ=255)
    tGGGAAGGCGTACTGATGAACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCt                                                               >  1:4365533/1‑71 (MQ=255)
     gggAAGGCGTACTGATGAACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTg                                                              >  1:4636450/1‑71 (MQ=255)
      ggAAGGCGTACTGATGAACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCt                                                               <  1:775736/69‑1 (MQ=255)
       gAAGGCGTACTGATGAACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCg                                                            >  1:1761992/1‑71 (MQ=255)
         aGGCGGACTGATGAACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCttt                                                          <  1:2925496/71‑1 (MQ=255)
          ggCGTACTGATGAACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGttt                                                         >  1:2359874/1‑71 (MQ=255)
           gCGTACTGATGAACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCt                                                               <  1:514854/64‑1 (MQ=255)
              tACTGATGAACCCGCCAATTAATTTTTTGCCCTGGCTACAGCAACGCCGGACCGCtt                                                                   <  1:3798977/57‑1 (MQ=255)
              tACTGATGAACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTgg                                                     <  1:2306804/71‑1 (MQ=255)
                cTGATGAACCCGCGAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGtt                                                   >  1:1631137/1‑71 (MQ=255)
                cTGATGAACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTc                                                                <  1:2453437/58‑1 (MQ=255)
                cTGATGAACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGtt                                                   >  1:2496444/1‑71 (MQ=255)
                cTGATGAACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGtt                                                   >  1:2451394/1‑71 (MQ=255)
                  gatgaACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCTGGACTGCTTTTCTTCGttt                                                         >  1:2816422/1‑63 (MQ=255)
                  gatgaACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACTGCTTTTCTTCGttt                                                         >  1:2139866/1‑63 (MQ=255)
                   atgaACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCg                                                                      >  1:2407948/1‑49 (MQ=255)
                   atgaACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGTTg                                                  <  1:4592762/69‑1 (MQ=255)
                   atgaACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCt                                                <  1:2136395/71‑1 (MQ=255)
                   atgaACACGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCt                                                >  1:280425/1‑71 (MQ=255)
                    tgaACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTg                                               <  1:1602654/71‑1 (MQ=255)
                    tgaACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTg                                               <  1:1325563/71‑1 (MQ=255)
                     gaACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCt                                                <  1:3550827/69‑1 (MQ=255)
                      aaCCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCtttt                                                                 >  1:1968420/1‑51 (MQ=255)
                      aaCCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGa                                              <  1:3300654/70‑1 (MQ=255)
                       aCCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATg                                            <  1:1531804/71‑1 (MQ=255)
                          cGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTc                                         <  1:1612067/71‑1 (MQ=255)
                           gCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCg                                        <  1:926482/71‑1 (MQ=255)
                             cAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGa                                              >  1:3149100/1‑63 (MQ=255)
                               attaattTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGt                                                    <  1:2896231/55‑1 (MQ=255)
                               attaattTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTgc                                    <  1:3161564/71‑1 (MQ=255)
                               attaattTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTgc                                    <  1:2547736/71‑1 (MQ=255)
                                 taattTTTTGCCCTGGCAACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTgcgc                                  >  1:4337274/1‑71 (MQ=255)
                                  aattTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTg                                     >  1:3132480/1‑67 (MQ=255)
                                    ttttttGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTg                                                      <  1:4389614/48‑1 (MQ=255)
                                         tGCCCTGGCGACAGCAACGCCGGACTGCTTTTCTTCGTTTCTGGTTGCTGATGTTCGTTGCGCCTctgct                           <  1:3663990/70‑1 (MQ=255)
                                               ggCGACAGCAACGCCGGACTGCTTTTCTTCGTTTCt                                                       <  1:2154746/36‑1 (MQ=38)
                                                gCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCg                   >  1:4155724/1‑71 (MQ=255)
                                                 cGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGt                  >  1:239/1‑71 (MQ=255)
                                                     aGCAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGTCgg               >  1:312980/1‑70 (MQ=255)
                                                     aGCAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGTCgg               >  1:1219729/1‑70 (MQ=255)
                                                       cAACGCCGGACCGCTTTTCTGGGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGTCGGGAt            <  1:478031/71‑1 (MQ=255)
                                                       cAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGTCGGGa             >  1:3810207/1‑70 (MQ=255)
                                                       cAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGTCGGGa             >  1:3815601/1‑70 (MQ=255)
                                                       cAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGTCGGGAt            <  1:2763628/71‑1 (MQ=255)
                                                       cAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGTCGGGAt            <  1:745501/71‑1 (MQ=255)
                                                        aaCGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGTCGGGAt            >  1:3181253/1‑70 (MQ=255)
                                                        aaCGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGTCGGGAt            >  1:2837487/1‑70 (MQ=255)
                                                         aCGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGt                  <  1:931360/63‑1 (MQ=255)
                                                          cGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTctgctgc                         <  1:3896442/55‑1 (MQ=255)
                                                          cGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGTCGGGATAAc         <  1:395812/71‑1 (MQ=255)
                                                          cGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGTCGGGATAAc         >  1:2098191/1‑71 (MQ=255)
                                                           gCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGTCGGGATAACg        >  1:2600560/1‑71 (MQ=255)
                                                             cGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGTCGGGATaa          >  1:2282436/1‑67 (MQ=255)
                                                             cGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGTCGGGATAACGCt      <  1:4034423/71‑1 (MQ=255)
                                                             cGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGTCGGGATAACGCt      <  1:4245627/71‑1 (MQ=255)
                                                             cGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGTCGGGATAACGCt      <  1:767364/71‑1 (MQ=255)
                                                              ggACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGATGGCCGTCGGGATAACGCt      >  1:4294051/1‑70 (MQ=255)
                                                               gACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTctgctgc                         >  1:3962830/1‑50 (MQ=255)
                                                               gACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTctgctgc                         >  1:2313065/1‑50 (MQ=255)
                                                               gACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGTCGGGATAACGCTaa    <  1:1887521/71‑1 (MQ=255)
                                                                 ccGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGTCGGGAt            >  1:2523840/1‑61 (MQ=255)
                                                                 ccGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGTCGGGATAACGCt      <  1:1228106/67‑1 (MQ=255)
                                                                 ccGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGTCGGGATAACGCTAATa  >  1:2502228/1‑71 (MQ=255)
                                                                  |                                                                     
CGCTTGGGAAGGCGTACTGATGAACCCGCCAATTAATTTTTTGCCCTGGCGACAGCAACGCCGGACCGCTTTTCTGCGTTTCTGGTTGCTGATGTTCGTTGCGCCTCTGCTGCTGGCCGTCGGGATAACGCTAATA  >  minE/2699105‑2699240

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: