Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,852,438 C→T 11.2% A395T (GCG→ACG)  ygcL ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,852,4380CT11.2% 225.8 / 16.7 80A395T (GCG→ACG) ygcLhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/41);  new base T (7/2);  total (37/43)
Fisher's exact test for biased strand distribution p-value = 7.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

ATGAACAGAGACTCCGGCCCCTTTGAAGTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTG  >  minE/1852373‑1852507
                                                                 |                                                                     
aTGAACAGAGACTCCGGCCCCTTTGAAGTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGTCTTGc                                                                  >  1:1362241/1‑71 (MQ=255)
aTGAACAGAGACTCCGGCCCCTTTGAAGTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGCCtt                                                                    <  1:2380552/69‑1 (MQ=255)
aTGAACAGAGACTCCGGCCCCTTTGAAGTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGCCtt                                                                    <  1:2955707/69‑1 (MQ=255)
 tGAACAGAGACTCCGGCCCCTTTGAAGTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGCCTTGCg                                                                 >  1:1553786/1‑71 (MQ=255)
     atgagaCTGCGGCCCCTTTGAAGTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTaa                                                              <  1:2121405/68‑1 (MQ=255)
      agagACTCCGGCCCCTTTGAAGTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAgg                                                            >  1:3109816/1‑71 (MQ=255)
              cGGCCCCTTTGAAGTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTta                                                    >  1:944444/1‑71 (MQ=255)
              cGGCCCCTTTGAAGTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTta                                                    >  1:3465085/1‑71 (MQ=255)
               ggCCCCTTTGAAGTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTtat                                                   <  1:3703216/71‑1 (MQ=255)
                       tGAAGTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCa                                             <  1:22768/69‑1 (MQ=255)
                         aaGTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAAcc                                         <  1:2141484/71‑1 (MQ=255)
                          aGTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGTCTTGCGTAAGGCTGTTTTATATCCCAAAcc                                         >  1:129930/1‑70 (MQ=255)
                           gTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATccc                                              >  1:2503752/1‑64 (MQ=255)
                              tttATTTTTAAACCCTTCTGCAAAGGTATATAACGCCt                                                                     <  1:1469858/38‑1 (MQ=255)
                              tttATTTTTAAACCCTTCTGCAAAGGTATATAACGCCt                                                                     <  1:2363283/38‑1 (MQ=255)
                               ttATTTTTAAACCCTTCTGCAAAGGTATATAACGTCTTGCGTAAGGCTGttt                                                      <  1:3357218/52‑1 (MQ=255)
                               ttATTTTTAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGt                                                        <  1:1666180/50‑1 (MQ=255)
                               ttATTTTTAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGt                                   >  1:3106002/1‑71 (MQ=255)
                                tATTTTTAAACCCTTCTGCAAAGGTATATAACGTCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGt                                   >  1:3202885/1‑70 (MQ=255)
                                tATTTTTAAACCCTTCTGCAAAGGTATATAACGTCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGt                                   >  1:3198532/1‑70 (MQ=255)
                                tATTTTTAAACCCTTCTGCAAAGGTATATAACGTCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGt                                   >  1:2876493/1‑70 (MQ=255)
                                tATTTTTAAACCCTTCTGCAAAGGTATATAACGCCTTGCGt                                                                <  1:1544725/41‑1 (MQ=255)
                                   ttttAAACCCTTCTGCAAAGGTATATAACGCTTTGCGTa                                                               <  1:126806/39‑1 (MQ=255)
                                   ttttAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTa                                                               <  1:1512655/39‑1 (MQ=255)
                                   ttttAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTa                                                               <  1:2668566/39‑1 (MQ=255)
                                   ttttAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCAc                                >  1:1206660/1‑70 (MQ=255)
                                     ttAAACCCTTCTGCAAAGGTATATAACGCCTTGCg                                                                 >  1:495315/1‑35 (MQ=255)
                                     ttAAACCCTTCTGCAAAGGTATATAACGCCTTGCGt                                                                >  1:2063841/1‑36 (MQ=255)
                                     ttAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAg                                                             >  1:262998/1‑39 (MQ=255)
                                     ttAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGt                                                        <  1:2483062/44‑1 (MQ=255)
                                     ttAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTAt                             >  1:2132274/1‑71 (MQ=255)
                                      tAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGtt                                                       >  1:1963821/1‑44 (MQ=255)
                                      tAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGtt                                                       >  1:812763/1‑44 (MQ=255)
                                        aaCCCTTCTGCAAAGGTATATAACGTCTTGCGTAAGGCTGttt                                                      <  1:1381108/43‑1 (MQ=255)
                                        aaCCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTAttt                           <  1:2315403/70‑1 (MQ=255)
                                        aaCCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTAttt                           <  1:630576/70‑1 (MQ=255)
                                         aCCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTtata                                                  <  1:3719977/46‑1 (MQ=255)
                                         aCCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCAc                                >  1:3502953/1‑64 (MQ=255)
                                         aCCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCAc                                >  1:25436/1‑64 (MQ=255)
                                         aCCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCAc                                >  1:522037/1‑64 (MQ=255)
                                         aCCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCg                         <  1:1521520/71‑1 (MQ=255)
                                         aCCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCg                         <  1:2932711/71‑1 (MQ=255)
                                           ccTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACa                                     <  1:3600599/57‑1 (MQ=255)
                                           ccTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACa                                     <  1:1651112/57‑1 (MQ=255)
                                           ccTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACa                                     <  1:503362/57‑1 (MQ=255)
                                             ttCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTa                     >  1:594727/1‑71 (MQ=255)
                                              tCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTAt                    >  1:3611593/1‑71 (MQ=255)
                                              tCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTAt                    >  1:2303735/1‑71 (MQ=255)
                                               cTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAAc                                          >  1:903790/1‑48 (MQ=255)
                                               cTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAAc                                          >  1:112288/1‑48 (MQ=255)
                                               cTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTa                     <  1:1729611/69‑1 (MQ=255)
                                               cTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTa                     <  1:174648/69‑1 (MQ=255)
                                               cTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTAt                    >  1:732453/1‑70 (MQ=255)
                                                 gCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGt                                   <  1:2821360/53‑1 (MQ=255)
                                                 gCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGt                                   <  1:3310911/53‑1 (MQ=255)
                                                  cAAAGGTATATAACGCCTTGCGTAAGGCTGTTTta                                                    >  1:2684863/1‑35 (MQ=255)
                                                  cAAAGGTATATAACGCCTTGCGTAAGGCTGTTTta                                                    >  1:2042082/1‑35 (MQ=255)
                                                  cAAAGGTATATAACGCCTTGCGTAAGGCTGTTTta                                                    >  1:3004223/1‑35 (MQ=255)
                                                    aaGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATc                   >  1:324889/1‑66 (MQ=255)
                                                    aaGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACAtt              >  1:2656474/1‑71 (MQ=255)
                                                    aaGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACAtt              >  1:59405/1‑71 (MQ=255)
                                                      ggTATATAACGCCTTGCGTAAGGCTGTTTTATATcc                                               <  1:3689407/36‑1 (MQ=255)
                                                        tatataACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCg          <  1:1865561/71‑1 (MQ=255)
                                                        tatataACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCg          <  1:2245038/71‑1 (MQ=255)
                                                          tataACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCaa                                       <  1:3147199/40‑1 (MQ=255)
                                                          tataACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCaa                                       <  1:924606/40‑1 (MQ=255)
                                                          tataACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGt         <  1:1185052/70‑1 (MQ=255)
                                                           ataACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTAt                    <  1:726151/58‑1 (MQ=255)
                                                           ataACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTAt                    <  1:2052848/58‑1 (MQ=255)
                                                            taACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTa        <  1:1860424/69‑1 (MQ=255)
                                                            taACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTAtt      >  1:1634236/1‑71 (MQ=255)
                                                            taACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTAtt      <  1:2205134/71‑1 (MQ=255)
                                                             aaCGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGt                        <  1:1866761/52‑1 (MQ=255)
                                                             aaCGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGt                        <  1:2406548/52‑1 (MQ=255)
                                                             aaCGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTa                     <  1:3085277/55‑1 (MQ=255)
                                                              aCGTCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTg             >  1:3582379/1‑62 (MQ=255)
                                                              aCGTCTAGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTAttgt    >  1:2481270/1‑71 (MQ=255)
                                                              aCGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTa                     <  1:2218562/54‑1 (MQ=255)
                                                              aCGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTAttg     <  1:2939102/70‑1 (MQ=255)
                                                                gCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTAttgttg  >  1:2814364/1‑71 (MQ=255)
                                                                 |                                                                     
ATGAACAGAGACTCCGGCCCCTTTGAAGTCTTTATTTTTAAACCCTTCTGCAAAGGTATATAACGCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTG  >  minE/1852373‑1852507

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: