Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 53,789 C→T 26.2% K333K (AAG→AAA imp ← exported protein required for envelope biosynthesis and integrity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE53,7890CT26.2% 137.2 / 43.5 114K333K (AAG→AAAimpexported protein required for envelope biosynthesis and integrity
Reads supporting (aligned to +/- strand):  ref base C (52/32);  new base T (19/11);  total (71/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

Only 100 of 146 total aligned reads displayed.

TGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGTAGA  >  minE/53723‑53850
                                                                  |                                                             
tgtgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                                                             <  1:2647594/66‑1 (MQ=255)
tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                                                             <  1:3234362/69‑1 (MQ=255)
tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                                                             <  1:2817122/69‑1 (MQ=255)
tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                                                             <  1:2977069/69‑1 (MQ=255)
tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                                                             <  1:3144948/69‑1 (MQ=255)
tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                                                             <  1:2313475/69‑1 (MQ=255)
tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                                                             <  1:2311756/69‑1 (MQ=255)
tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                                                             <  1:2270104/69‑1 (MQ=255)
tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                                                             <  1:3203727/69‑1 (MQ=255)
tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                                                             <  1:192655/69‑1 (MQ=255)
tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                                                             <  1:3223025/69‑1 (MQ=255)
tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                                                             <  1:3380199/69‑1 (MQ=255)
tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                                                             <  1:145483/69‑1 (MQ=255)
tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                                                             <  1:65872/69‑1 (MQ=255)
tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                                                             <  1:1017530/69‑1 (MQ=255)
tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                                                             <  1:697937/69‑1 (MQ=255)
tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                                                             <  1:869080/69‑1 (MQ=255)
 gcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                                                             <  1:631914/68‑1 (MQ=255)
 gcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                                                             <  1:102086/68‑1 (MQ=255)
 gcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                                                             <  1:1171011/68‑1 (MQ=255)
 gcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                                                             <  1:3190794/68‑1 (MQ=255)
 gcgtAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTg                                                            >  1:3060091/1‑69 (MQ=255)
 gcgtAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTg                                                            >  1:220889/1‑69 (MQ=255)
            aGTACTGGAACCGTACTTGTTATCGAAATCTTTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCgc                                                  >  1:620848/1‑67 (MQ=255)
            aGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGAc                                                 >  1:2014728/1‑69 (MQ=255)
            aGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGAc                                                 >  1:3128899/1‑69 (MQ=255)
            aGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGAc                                                 >  1:1216823/1‑69 (MQ=255)
            aGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGAc                                                 >  1:2165881/1‑69 (MQ=255)
            aGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGAc                                                 >  1:319390/1‑69 (MQ=255)
            aGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGAc                                                 >  1:1930901/1‑69 (MQ=255)
            aGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGAc                                                 >  1:696079/1‑69 (MQ=255)
            aGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGAc                                                 >  1:1515753/1‑69 (MQ=255)
            aGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGAc                                                 >  1:1868940/1‑69 (MQ=255)
              tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt                                               <  1:1193483/69‑1 (MQ=255)
              tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt                                               <  1:1008402/69‑1 (MQ=255)
              tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt                                               <  1:2594279/69‑1 (MQ=255)
              tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt                                               <  1:529168/69‑1 (MQ=255)
              tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt                                               <  1:2523425/69‑1 (MQ=255)
              tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt                                               <  1:1607940/69‑1 (MQ=255)
              tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt                                               <  1:38647/69‑1 (MQ=255)
              tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt                                               <  1:1757760/69‑1 (MQ=255)
              tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt                                               <  1:3378098/69‑1 (MQ=255)
              tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt                                               <  1:1245649/69‑1 (MQ=255)
              tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt                                               <  1:200293/69‑1 (MQ=255)
              tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTACTCGACGt                                               <  1:2571563/69‑1 (MQ=255)
               aCTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt                                               <  1:1387042/68‑1 (MQ=255)
                cTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGtt                                              >  1:1199632/1‑68 (MQ=255)
                cTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGtt                                              >  1:899636/1‑68 (MQ=255)
                cTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGtt                                              >  1:1965963/1‑68 (MQ=255)
                  ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                                                            >  1:841725/1‑52 (MQ=255)
                  ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                                                            >  1:387248/1‑52 (MQ=255)
                  ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                                                            >  1:2765329/1‑52 (MQ=255)
                  ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                                                            >  1:2820069/1‑52 (MQ=255)
                  ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                                                            >  1:2022502/1‑52 (MQ=255)
                  ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                                                            >  1:191660/1‑52 (MQ=255)
                  ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                                                            >  1:1770557/1‑52 (MQ=255)
                  ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                                                            >  1:1337653/1‑52 (MQ=255)
                  ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                                                            >  1:1569004/1‑52 (MQ=255)
                  ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                                                            >  1:1507414/1‑52 (MQ=255)
                  ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                                                            >  1:1452030/1‑52 (MQ=255)
                  ggAACCGTATTTGTTATCGAAGTCATTCAAGTAGCTAGGATCGCTGACTTTg                                                            >  1:2935780/1‑52 (MQ=255)
                  ggAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGa                                            <  1:2525945/68‑1 (MQ=255)
                  ggAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGa                                            <  1:563198/68‑1 (MQ=255)
                  ggAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGa                                            <  1:2698557/68‑1 (MQ=255)
                      ccGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGc                                       <  1:2248658/69‑1 (MQ=255)
                           cTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt                                               >  1:806029/1‑56 (MQ=255)
                           cTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCacac                                  >  1:1378874/1‑69 (MQ=255)
                           cTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCacac                                  >  1:1316477/1‑69 (MQ=255)
                           cTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCacac                                  >  1:1606981/1‑69 (MQ=255)
                           cTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCacac                                  >  1:2069906/1‑69 (MQ=255)
                                               aGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAg               <  1:1163199/68‑1 (MQ=255)
                                                         aTCGCTGACTTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGt     <  1:3179292/68‑1 (MQ=255)
                                                         aTCGCTGACTTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGt     <  1:3103472/68‑1 (MQ=255)
                                                         aTCGCTGACTTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGt     <  1:2837596/68‑1 (MQ=255)
                                                         aTCGCTGACTTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGt     <  1:2783583/68‑1 (MQ=255)
                                                         aTCGCTGACTTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGt     <  1:1289935/68‑1 (MQ=255)
                                                         aTCGCTGACTTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGt     <  1:795842/68‑1 (MQ=255)
                                                         aTCGCTGACTTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGt     <  1:443751/68‑1 (MQ=255)
                                                         aTCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAg      <  1:1596064/67‑1 (MQ=255)
                                                          tCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTcc        >  1:329070/1‑64 (MQ=255)
                                                          tCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTcc        >  1:2405686/1‑64 (MQ=255)
                                                          tCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGTa    >  1:53536/1‑68 (MQ=255)
                                                          tCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGTAg   <  1:1712932/69‑1 (MQ=255)
                                                          tCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGTAg   >  1:1822056/1‑69 (MQ=255)
                                                          tCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGTAg   >  1:926737/1‑69 (MQ=255)
                                                          tCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGTAg   <  1:293371/69‑1 (MQ=255)
                                                          tCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGTAg   >  1:2860360/1‑69 (MQ=255)
                                                          tCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGTAg   >  1:2834852/1‑69 (MQ=255)
                                                          tCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGTAg   >  1:796733/1‑69 (MQ=255)
                                                           cGCTGAGCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGTAGa  >  1:439484/1‑69 (MQ=255)
                                                           cGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGTAGa  >  1:342300/1‑69 (MQ=255)
                                                           cGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGTAGa  >  1:983007/1‑69 (MQ=255)
                                                           cGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGTAGa  >  1:986585/1‑69 (MQ=255)
                                                           cGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGTAGa  >  1:821875/1‑69 (MQ=255)
                                                           cGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGTAGa  >  1:575099/1‑69 (MQ=255)
                                                           cGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGTAGa  >  1:1702744/1‑69 (MQ=255)
                                                           cGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGTAGa  >  1:2396501/1‑69 (MQ=255)
                                                           cGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGTAGa  >  1:313851/1‑69 (MQ=255)
                                                           cGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGTAGa  >  1:3006243/1‑69 (MQ=255)
                                                           cGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGTAGa  >  1:2776526/1‑69 (MQ=255)
                                                                  |                                                             
TGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACACCTGATCCATGACCCCGGAGTGGTTCCAGTAGA  >  minE/53723‑53850

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: